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1.0  Sea  Ice  Mechanics  Research 

The  work  conducted  under  this  Contract  has  covered  Arctic  field  measurements,  sea  ice 
model  development,  and  the  organization  of  a  workshop  on  Sea  Ice  Mechanics  and  Arctic 
Modeling.  The  research  is  best  summarized  by  the  published  papers.  They  are  presented  in  the 
Appendix,  and  the  titles  are  listed  here.  A  discussion  of  the  workshop  follows  this  section. 

2.0  Published  papers  (full  text  with  Appendix) 

“Pack  Ice  Anisotropic  Constitutive  Model,’’  M.D.  Coon,  D.C.  Echert,  and  G.S.  Knoke  (1992), 
IAHR  Ice  Symposium,  Banff,  Alberta,  pp  1 188-1205. 

“Microwave  Response  to  Stress  in  Sea  Ice  -  A  Report  on  Initial  Field  Tests,”  M.D.  Coon,  R.G. 
Onstott,  D.C.  Echert  (1993),  Cold  Regions  Science  and  Technology,  pp  79-83. 

“Phenomenological  Constitutive  Model  for  Columnar  Ice,”  M.D.  Coon  and  G.S.  Knoke  (1993), 
ISOPE-93  Singapore:  3rd  International  Offshore  and  Polar  Engineering  Conference, 
Singapore. 

“Contemporaneous  Field  Measurements  of  Pack  Ice  Stress  and  Ice  Strain  Measurements  from 
SAR  Imagery,”  M.D.  Coon,  G.S.  Knoke,  D.C.  Echert,  and  H.L.  Stem  (1993), 
OCEANS'93,  Oct.  18-21,  Victoria,  British  Columbia. 

“The  Sea  Ice  Mechanics  Initiative,”  M.D.  Coon,  G.S.  Knoke,  and  D.C.  Echert  (1994),  26th 
Annual  Offshore  Technology  Conference  (OTC),  Houston,  TX,  pp  499-506. 

“Sea  Ice  Mechanics  Research,”  D.M.  Coon,  D.C.  Echert,  and  G.S.  Knoke  (1995),  Sea  Ice 
Mechanics  and  Arctic  Modeling  Workshop,  Anchorage,  AK,  pp  150-159. 

“Panel  Report  on  100  km  Scale  (Pack  Ice),”  M.D.  Coon,  R.  Preller,  A.  Proshutinsky  (1995),  Sea 
Ice  Mechanics  and  Arctic  Modeling  Workshop,  Anchorage,  AK,  pp  10-15. 

“Sea  Ice  Deformation  and  Stress,  A  Comparison  across  Space  Scales,”  M.D.  Coon,  D.C.  Echert, 
G.S.  Knoke,  J.E.  Overland,  S.  Salo,  R.S.  Pritchard,  D.A.  Rothrock,  and  H.L.  Stern 
(1995),  Sea  Ice  Mechanics  and  Modeling  Workshop,  Anchorage,  AK,  pp  216-229. 

“Acoustic  Signatures  of  Falling  Ice  Blocks,”  R.S.  Pritchard,  Y.Xie,  M.D.  Farmer,  M.D.  Coon, 
G.S.  Knoke,  and  D.C.  Echert  (1995),  AMD-Ice  Mechanics,  pp  219-230. 

“Force-Displacement  Measurements  of  a  First-Year  Pressure  Ridge  Keel,”  D.M.  Coon,  D.C. 
Echert,  and  G.S.  Knoke  (1995),  AMD-Ice  Mechanics,  pp  239-254. 

“An  Oriented  Thickness  Distribution  for  Sea  Ice,”  D.M.  Coon,  D.C.  Echert,  and  G.S.  Knoke  and 
R.S.  Pritchard  (undated),  submitted  to  JGR,  1997. 


3.0  The  Sea  Ice  Mechanics  and  Arctic  Modeling  Workshop 

Northwest  Research  Associates,  Inc.  (NWRA)  organized  the  Sea  Ice  Mechanics  and 
Arctic  Modeling  Workshop  and  published  the  2-volume  proceedings. 

The  Sea  Ice  Mechanics  and  Arctic  Workshop  was  held  in  Anchorage,  Alaska, 
April  25-28,  1995.  The  workshop  focused  on  current  state  of  the  practice  and  future  research 
needs  relative  to  offshore  oil  and  gas  facilities,  and  also  reviewed  the  results  of  the  U.S.  Navy’s 
Office  of  Naval  Research  (ONR)  Sea  Ice  Mechanics  Initiative  (SIMI).  SIMI  was  an  Office  of 
Naval  Research  Accelerated  Research  Initiative  (ARI),  spanning  the  years  1991  to  1996.  Its 
goals  were  to: 

•  understand  sea  ice  constitutive  laws  and  fracture  mechanics  over  the  full  range  of 
geophysical  scales;  and 

•  determine  the  scaled  response  to  applied  external  forces  and  develop  physically  based 
constitutive  and  fracture  models. 

The  goals  of  SIMI  were  established  at  a  planning  meeting  held  at  Airlie,  Virginia,  in 
November  1990.  Oil  industry  representatives  participated  in  this  planning  process.  The  report 
from  that  meeting  is  included  as  Appendix  A  in  volume  1  of  these  Proceedings.  The  main  field 
experiments  were  defined  at  a  second  SIMI  planning  meeting  held  at  Sidney,  British  Columbia, 
in  August  1993  and  are  described  in  the  Summary  Plan  included  as  Appendix  B  in  Volume  1. 
The  main  field  experiments  were  conducted  in  the  Beaufort  Sea  during  1993  and  1994.  The 
Proceedings  of  the  workshop  are  presented  in  Volume  2.  The  Table  of  Contents  of  Vol.  1&2  are 
presented  below. 

The  workshop  received  financial  support  from  a  number  of  government  agencies  and  oil 
companies,  including: 

U.S.  Navy  Office  of  Naval  Research 
U.S.  Minerals  Management  Service 
Canadian  National  Energy  Board 
Amoco,  Arco,  Chevron,  Exxon,  and  Mobil 

The  workshop  sponsored  travel  and  expenses  for  five  Russian  scientists  and  provided 
partial  support  for  a  sixth  Russian  scientist  to  attend  the  workshop.  Funding  for  these  individuals 
was  provided  by  Minerals  Management  Service,  Arco  and  Mobil.  These  scientists  were: 

•  Prof.  Robert  V.  Goldstein  of  the  Institute  for  Problems  in  Mechanics,  Moscow, 
Russia; 

•  Dr.  Anatoly  M.  Polomoshnov  of  the  Sakhalin  Oil  and  Gas  Institute,  Sakhalin  Island, 
Russia; 

•  Dr.  Vladimir  E.  Ryabinin  of  the  Hydropmetcentre  of  Russia,  Moscow,  Russia; 

•  Dr.  Victor  N.  Smirnov  of  the  Arctic  and  Antarctic  Research  Institute,  St.  Petersburg, 
Russia; 
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•  Dr.  Pavel  A.  Truskov  of  the  Sakhalin  Oil  and  Gas  Institute,  Sakhalin  Island,  Russia; 
and 

•  Dr.  Alexander  Bekker  of  the  Far-Eastern  State  Technical  University,  Vladivostok, 
Russia. 

The  workshop  began  with  a  day  of  providing  information  to,  and  among,  the  attendees. 
On  Tuesday  morning,  a  series  of  speakers  provided  background  information  to  all  participants 
That  afternoon  was  a  poster  session  in  which  43  posters  were  presented  and  discussed.  The 
remainder  of  the  week  was  spent  in  working  panel  sessions  and  in  presenting  results  to,  and 
discussing  them  with,  the  other  attendees.  On  Wednesday  morning  (April  26),  five  panel 
sessions  were  held,  each  focusing  on  processes  at  different  physical  scales.  On  Thursday 
morning,  four  panels  were  held  on  topics  of  current  importance.  The  afternoons  of  Wednesday 
and  Thursday  were  spent  in  presenting  and  discussing  the  results  of  the  morning  sessions  to  the 
other  attendees.  Friday  morning  was  a  plenary  wrap-up  and  summary  session  of  the  workshop 
results. 

A  summary  of  workshop  results,  “Recent  Accomplishments  and  Research  Needed,”  is 
presented  after  the  Table  of  Contents  of  Vol.  1  and  2. 
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3.1  Proceedings  of  the  Sea  Ice  Mechanics  and  Arctic  Modeling  Workshop 


Sea  Ice  Mechanics  and  Arctic  Modeling  Workshop 

Volume  1 

Table  of  Contents 


1.  Introduction 

2.  Acoustics 

Farmer,  D.,  and  Y.  Xie,  "Acoustic  and  Seismic  Studies  of  Ice  Mechanics" 
Mikhalevsky,  P„  A.  B.  Baggeroer,  H.  Schmidt,  K.  von  der  Heydt,  E.  K. 

Sheer,  and  A.  Gavrilov,  "Transarctic  Acoustic  Propagation" 
Pritchard,  R.  S.,  "Sea  Ice  Failure  Mechanisms" 

Rajan,  S.  D.,  "Sea  Ice  Mechanics  Research:  Tomographic  Imaging  of 
Wave  Speeds  and  Acoustic  Emission  Event  Localization" 

Schmidt,  H„  A.  B.  Baggeroer,  I.  Dyer,  K.  von  der  Heydt,  and  E.  K. 

Scheer,  "Seismo- Acoustic  Remote  Sensing  of  Ice-Mechanical 
Processes  in  the  Arctic" 

Stein,  P.  J.,  D.  W.  Andersen,  A.  Bahlavouni,  S.  E.  Euerle,  G.  M.  Santos, 
and  R.  K.  Menoche,  "SIMI  Winter-Over  Geophone/Hydrophone 
System" 

3.  Fracture  Mechanics 

Bazant,  Z.  P.,  Y.-N.  Li,  M.  Jirdsek,  Z.  Li,  and  J.-J.  Kim,  "Effect  of  Size  on 
Distributed  Damage  and  Fracture  of  Sea  Ice" 

Dempsey,  J.,  and  R.  Adamson,  "Scale  Effects  on  the  Fracture  and 
Constitutive  Behavior  of  Sea  Ice" 

4.  Ice  Properties 

Cole,  D.  M.,  "Field  and  Laboratory  Experiments  and  Modeling  of  the 
Constitutive  Behavior  of  Sea  Ice" 

Gupta,  V.,  R.  C.  Picu,  J.  Bergstrom,  and  H.  J.  Frost,  "Crack  Nucleation 
Mechanisms  in  Columnar  Ice  -  Recent  Developments" 

Schulson,  E.  M.,  "Compressive  Failure  of  Columnar  Saline  Ice  under 
Multiaxial  Loading" 

Shapiro,  L.  H.,  W.  F.  Weeks,  and  W.  D.  Harrison,  "Studies  of  the 

Influence  of  Fabric  and  Structure  on  the  Flexural  Strength  of  Sea 
Ice  and,  of  the  Consolidation  of  First-Year  Pressure  Ridges" 
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5.  Ice  Stress,  Ice  Strain,  and  Ice  Conditions 

Coon,  M.  D.,  D.  C.  Echert,  and  G.  S.  Knoke,  "Sea  Ice  Mechanics 
Research" 

O’Hara,  S.  and  J.  Ardai,  Jr.,  "SIMI  GPS  Position  and  CTD  Cast  Data" 
Overland,  J.,  S.  Salo,  S.  Li,  and  L.  McNutt,  "Regional  and  Floe-Floe 
Deformation" 

Richter-Menge,  J.  A.,  B.  C.  Elder,  W.  B.  Tucker  IE,  and  D.  K.  Perovich 
"Pack  Ice  Stresses  and  their  Relationship  to  Regional  Deformation” 

6.  Modeling 

Connor,  J.  J.,  S.  Shyam-Sunder,  A.  Elvin,  D.  Choi,  and  J.  Kim,  "Physically 
Based  Constitutive  Modeling  of  Ice" 

Hopkins,  M.  A.,  "Numerical  Simulation  of  Arctic  Pressure  Ridging" 

Lewis,  J.  K„  and  P.  J.  Stein,  "Sea  Ice  Mechanics  Related  to  Thermally- 
Induced  Stresses  and  Fracturing  of  Pack  Ice" 

Petrenko,  V.,  and  O.  Gluschenkov,  "Measurements  of  Crack  Velocity  in 
Sea  Ice  Using  Electromagnetic  Techniques" 

Rodin,  G.,  R.  Shapery,  K.  Abdel-Tawab,  and  L.  Wang,  "Constitutive 
Equations  and  Fracture  Models  for  Sea  Ice” 

Wu,  M.  S.,  J.  Niu,  Y.  Zhang,  and  H.  Zhou,  "Physically-Based  Constitutive 
Modeling  of  Ice:  Damage  and  Failure" 

Appendix  A  —  Air-Ice-Ocean  Interaction:  Lead  Dynamics,  Ice  Mechanics,  Ice  Acoustics 
(Sea  Ice  Mechanics  Workshop  at  Airlie,  VA,  1990) 

Appendix  B  --  Sea  Ice  Mechanics  Initiative  (SIMI)  Summary  Plan,  FY94-95  (Workshop 
at  Sidney,  BC,  1993) 
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3.1  Proceedings  of  the  Sea  Ice  Mechanics  and  Arctic  Modeling  Workshop 

Volume  2 


TABLE  OF  CONTENTS 

Preface  i 

Introduction  1 

1.  PART  1-RESULTS  OF  WORKSHOP  SESSIONS  5 

1.1  Summary  of  Workshop  Results  -  Recent  Accomplishments  and  Research 

Needed  5 

1.2  Results  of  Panel  Discussions  5 

1.2. 1  100  km  Scale  (pack  ice)  10 

1.2.2  10  km  Scale  (floe  cluster)  16 

1.2.3  1  km  Scale  (floe  ridge  crack)  22 

1.2.4  10  m  Scale  (macro-crack)  28 

1.2.5  1  cm  to  1  m  Scale  (materials)  45 

1.2.6  Ice  Loads  on  Fixed  and  Floating  Structures,  Including  Vessels  60 

1.2.7  Ice  Forecasting  for  Ship  Routing  and  Operations  79 

1.2.8  Ice  Gouging  of  the  Sea  Floor  84 

1.2.9  Oil  Spills  and  Oil  Transport  by  the  Ice  97 

2.  PART  II  -  INVITED  SPEAKERS  AND  POSTERS  102 

2.1  Speakers  on  Tuesday,  April  25,  Opening  Session  102 

2.1.1  “Industry-Government  Cooperation  in  Ice  Mechanics  and  Modeling,” 

Albert  Wang  (Exxon)  103 

2.1.2  “Program  for  Energy  Research  and  Development  (PERD),  Ice-Structure 
Interaction  Program,”  Ibrahim  Konuk,  (National  Energy  Board  of 

Canada)  118 

2.1.3  “MMS  Programs,”  Jeff  Walker  (MMS)  (no  text  available) 

2.1.4  “Application  of  Ice  Cover  Mechanics  in  Design  and  Operations  of  Marine 
Structures”,  Jukka  Tuhkuri,  (Arctic  Offshore  Research  Centre,  Helsinki 

Univ.  of  Tech.)  123 

2.1.5  “Ice- Structure  Interaction  Research  and  Arctic  Development,”  Ben 

Danielewicz  (Canadian  Marine  Drilling)  135 

2.1.6  “BP  Program,”  Graham  A.  N.  Thomas  (BP  International  Ltd.)  141 

2.2  Invited  Speaker  at  Banquet,  Thursday,  April  27 

“A  Review  of  Select  Exploration  and  Development  Activities  in  State 
Nearshore  Areas,”  James  Hansen,  Manager  of  Oil  and  Gas  Leasing  Program, 
Alaska  Department  of  Natural  Resources.  145 


6 


2.3  SIMI  Papers  not  included  in  Volume  1  175 

“Ice  Surface  Oscillation  Measurements  on  SIMI  using  Strain,  Heave  and  Tilt  Sensors,” 

Peter  Wadhams  and  Stephen  C.  S.  Wells,  Scott  Polar  Research  Institute,  Cambridge,  UK.  176 

“AUV  Operations  in  the  Arctic,”  J.  G.  Bellingham,  J.  J.  Leonard,  J.  Vaganay,  C.  A. 

Goudey,  D.  K.  Atwood,  T.  R.  Consi,  J.  W.  Bales,  H.  Schmidt  and  C.  Chryssostomidis, 

Mass.  Inst,  of  Tech.  190 

“Multipath  Navigation  in  the  Arctic:  A  Feasibility  Test,”  Max  Deffenbaugh,  Henrik 

Schmidt  and  James  G.  Bellingham,  Mass.  Inst,  of  Tech.  199 

“Performance  Assessment  of  the  High  Resolution  GPS  Ranging  System,”  Max 

Deffenbaugh,  Henrik  Schmidt  and  James  G.  Bellingham,  Mass.  Inst  of  Tech.  204 

2.4  Poster  Session  207 

Author(s)  and  Title 

Anderson,  Robert  M.,  and  Burton  L.  Markham,  "Submarine  Force  Support  of  Arctic  212 

Science:  The  SCICEX  Program." 

Bazant,  Z.  P.,  Y.-N.  Li,  M.  Jirasek,  Z.  Li,  and  J.-J.  Kim,  "Effect  of  Size  on  VI,  73 

Distributed  Damage  and  Fracture  of  Sea  Ice." 

Bekker,  A.  T.,  "Ice  Structure  Dynamic  Interaction:  Destruction  Ice  Model."  213 

Bekker,  A.  T.,  and  S.  G.  Gomolsky  "Determination  of  the  Ice  Strength  for  Calculation  214 

of  the  Ice  Load." 

Cole,  D.  M.,  "Field  and  Laboratory  Experiments  and  Modeling  of  the  Constitutive  VI,  101 
Behavior  of  Sea  Ice." 

Cole,  David  M.,  Lewis  H.  Shapiro,  Wilford  F.  Weeks,  Carl  Byers,  John  P.  Dempsey,  215 

Robert  M.  Adamson,  Victor  F.  Petrenko,  and  Oleg  V.  Gluschenkov,  "The 
Barrow  Experiments." 

Connor,  J.  J.,  S.  Shyam-Sunder,  A.  Elvin,  D.  Choi,  and  J.  Kim,  "Physically  Based  VI,  189 

Constitutive  Modeling  of  Ice." 

Coon,  M.  D.,  D.  C.  Echert,  and  G.  S.  Knoke,  "Sea  Ice  Mechanics  from  20  cm  to  80  VI,  151 
km." 

Coon,  M.  D.,  D.  C.  Echert,  G.  S.  Knoke,  J.  E.  Overland,  S.  Salo,  R.  S.  Pritchard,  D.  216 

A.  Rothrock,  and  H.  L.  Stern  "Sea  Ice  Deformation  and  Stress,  a  Comparison 
across  Space  Scales." 

Dempsey,  J.,  and  R.  Adamson,  "Scale  Effects  on  the  Fracture  and  Constitutive  VI,  84 

Behavior  of  Sea  Ice." 

Farmer,  D.,  and  Y.  Xie,  "Acoustic  and  Seismic  Studies  of  Ice  Mechanics."  VI,  7 
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Gupta,  V.,  R.  C.  Picu,  J.  Bergstrom,  and  H.  J.  Frost,  "Crack  Nucleation  Mechanisms  VI,  1 10 
in  Columnar  Ice  --  Recent  Developments." 

Hedstrom,  Katherine,  "A  Coupled  Ice-Ocean-Atmosphere  Model  of  the  Bering  Sea. '  230 

Hopkins,  M.  A.,  "Numerical  Simulation  of  Arctic  Pressure  Ridging."  VI,  199 

Kerman,  Bryan,  "Information  States  in  Sea  Ice  Imagery."  239 

Kristensen,  Dirk,  Vera  Alexander,  Thomas  Royer,  and  Robert  Eisner,  "The  Planned  249 

NSF/UNOLS  Arctic  Research  Vessel  for  the  United  States  Arctic  Science". 

Lai,  Ronald  J.,  and  Walter  R.  Johnson,  "Minerals  Management  Service  Environmental  262 

Studies  Program  and  Oil  Spill  Risk  Analysis  for  the  Beaufort  Sea". 

Lewis,  J.  K.,  and  P.  J.  Stein,  "Sea  Ice  Mechanics  Related  to  Thermally-Induced  VI,  209 

Stresses  and  Fracturing  of  Pack  Ice." 

Mikhalevsky,  P„  A.  B.  Baggeroer,  H.  Schmidt,  K.  von  der  Heydt,  E.  K.  Sheer,  and  A.  VI,  24 
Gavrilov,  "Transarctic  Acoustic  Propagation." 

O'Hara,  S.,  and  J.  Ardai,  Jr.,  "SIMI  GPS  Position  and  CTD  Cast  Data."  VI,  160 

Ostoja-Starzewski,  Martin,  "Micromechanically  Based  Constitutive  Laws  for  Ice  263 

Fields." 

Overland,  J.,  S.  Salo,  S.  Li,  and  L.  McNutt,  "Regional  and  Floe-Roe  Deformation."  VI,  168 

Petrenko,  V.,  and  O.  Gluschenkov,  "Measurements  of  Crack  Velocity  in  Sea  Ice  VI,  216 

Using  Electromagnetic  Techniques." 

Pham,  Thomas  Thang,  "The  Polar  Ice  Prediction  System  at  Reet  Numerical  264 

Meteorology  and  Oceanography  Center  (FLENUMMETOCCEN)." 

Pritchard,  R.  S.,  "Sea  Ice  Failure  Mechanisms."  VI,  34 

Proshutinsky,  Andrey,  and  Mark  Johnson,  "The  Arctic  Ocean  Ice  and  Water  265 

Transport  During  1946-1988." 

Rajan,  S.  D.,  "Sea  Ice  Mechanics  Research,  "Tomographic  Imaging  of  Wave  Speeds  VI,  41 
and  Acoustic  Emission  Event  Localization." 

Richter-Menge,  J.  A.,  B.  C.  Elder,  W.  B.  Tucker  HI,  and  D.  K.  Perovich,  "Pack  Ice  VI,  178 
Stresses  and  their  Relationship  to  Regional  Deformation." 

Rodin,  G.,  R.  Schapery,  "Constitutive  Equations  and  Fracture  Models  for  Sea  Ice. "  VI,  226 

Ryabinin,  Vladimir  E.,  Alexander  I.  Danilov,  Vladimir  V.  Elisov,  Alexander  V.  276 

Klepikov,  Vladimir  A.  Kuidumov,  Valery  N.  Malek,  Victor  N.  Smirnov,  Igor  N. 

Stepanov,  Oleg  Ya.  Timofeev,  "Marine  Ice  Bottom  Scouring:  Some  Mechanisms 
and  an  Approach  to  Depth  Evaluation." 
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3.2  Summary  of  Workshop  Results 

Recent  Accomplishments  and  Research  Needed 

This  section  is  a  summary  of  the  Friday  (April  28)  morning  session.  The  objective  of  this 
session  was  to  describe  the  research  needed  in  each  workshop  topic  area: 

•  100  km  scale  (pack  ice) 

•  10  km  Scale  (floe  cluster) 

•  1  km  Scale  (floe,  ridge,  crack) 

•  10  meter  Scale  (macro-crack) 

•  1  cm  to  1  m  Scale  (micro-crack,  crystal) 

•  Ice  Loads  on  Fixed  and  Floating  Structures,  including  Vessels 

•  Ice  Forecasting  for  Ship  Routing  and  Operations 

•  Ice  Gouging  of  the  Sea  Floor 

•  Oil  Spills  and  Oil  Transport  by  the  Ice 

For  most  of  the  topic  areas,  the  most  important  accomplishment  by  SIMI  and  most 
important  area  for  future  research,  as  judged  by  the  panels,  are  contained  on  the  table  below. 
The  topics  on  ice  forecasting,  ice  gouging  of  the  seafloor,  and  oil  spills  were  not  specifically 
addressed  by  SIMI.  The  main  text  for  the  ice  forecasting  topic  includes  a  description  of  the 
current  state  of  the  practice.  For  the  ice  gouging  topics,  we  include  a  summary  of  what  currently 
is  known.  The  needs  and  findings  on  any  topic  are  best  determined  by  reading  the  panel  report 
on  that  topic. 

3.3  Results  of  Panel  discussions 

The  results  of  the  discussions  are  summarized  in  Table  1. 
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Table  1  -  Summary  of  Recent  Accomplishments  and  Research  Needed 
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Improved  wind  fields  from  atmospheric  models  has  been  applied 
for  input  to  ice-ocean  models. 
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Arctic  sea  ice  may  prove  to  be  an  important  indicator  of  global  climate  change. 
However,  the  inter-annual  variability  of  ice  extent  and  thickness  clouds  the  issue.  The 
focus  of  this  work  is  to  improve  the  accuracy  of  ice  models  so  that  effects  due  to 
climate  change  can  be  distinguished  from  inter-annual  variability.  This  paper  describes 
the  development  of  a  2-D  plastic  anisotropic  constitutive  model  of  the  ice  cover  which 
includes  the  creation  of  new  leads  in  thick  ice  and  the  effects  of  leads  and  their 
orientation  on  the  strength  of  the  ice  cover.  The  model  requires  an  oriented  thickness 
distribution  to  account  for  the  growth  and  ridging  of  ice  in  refrozen  leads.  The  shear 
and  compressive  strength  of  each  lead  are  functions  of  the  lead  ice  thickness 
distribution.  Plans  to  verify  the  new  model  with  ice  stress  field  measurements,  Arctic 
buoy  ice  motion  data,  ice  vectors,  and  lead  orientation  data  provided  by  synthetic 
aperture  radar  satellite  imagery,  are  discussed. 


INTRODUCTION 

J  An  understanding  of  the  heat,  mass,  and  momentum  balance  of  the  Arctic  has 
been  known  to  be  important  to  global  climate  change  for  over  30  years  (e.g.,  National 
Academy  of  Sciences,  1967  and  Barry,  1985  and  1986).  According  to  Parkinson  (1990), 
sea  ice  thickness  and  distribution  is  a  good  indicator  of  climate  change.  It  is  highly 
sensitive  to  changes  in  climate,  and  it  can  be  routinely  measured  on  a  global  scale 
through  satellite  remote  sensing.  Parkinson  identified  a  shortcoming  of  its  use  as  its 
poorly-understood  inter-annual  variability;  the  climate  signal  may  be  lost  in  the 
year-to-year  variation.  If  the  causes  of  the  variability  were  understood,  the  utility  of  sea 
ice  as  an  indicator  would  be  greatly  enhanced.  This  understanding  can  be  obtained 
through  improved  ice  dynamics  and  thermodynamics  models. 
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In  polar  regions,  the  sea  ice  cover  acts  as  a  barrier  between  the  atmosphere  and 
ocean,  curtailing  the  transfer  of  thermal  and  mechanical  energy  as  well  as  mass  trans¬ 
port.  The  energy  exchange  over  open  leads  and  young  sea  ice  in  the  Central  Arctic  has 
pronounced  effects  on  the  atmosphere  and  ocean.  For  example,  Badgley  (1966)  reports 
that  the  turbulent  energy  input  to  the  atmosphere  over  a  freezing  lead  is  more  than  2 
orders  of  magnitude  larger  than  that  over  typical  3  m  ice.  Maykut  (1978)  has  shown  that 
heat  input  to  the  atmosphere  over  0-40  cm  thick  ice  during  the  cold  months  is  between 
1  and  2  orders  of  magnitude  larger  than  that  over  perennial  ice.  Thus,  even  though 
leads  account  for  a  small  percentage  of  the  ice  pack,  their  effects  dominate  the  large 
scale  turbulent  heat  exchange  and  alter  boundary  layer  stability.  For  these  reasons, 
proper  characterization  of  lead  formation  and  fate  in  Arctic  climate  models  is  required. 

The  first  comprehensive  ice  dynamics  and  thermodynamics  model  was  developed 
by  the  Arctic  Ice  Dynamics  Joint  Experiment  (AIDJEX)  in  the  1970’s.  This  model  is 
summarized  in  Coon,  et  al.  (1974)  and  Coon  (1980).  The  model  contains  a  constitutive 
law  consisting  of  a  yield  surface  which  is  a  prescribed  function  of  the  yield  strength  and 
internal  ice  stress.  This  model,  and  subsequent  modifications  and  improvements  to  this 
model  (e.g.  Hibler,  1986),  assumed  that  sea  ice  is  isotropic,  i.e.,  that  leads  are  randomly 
distributed  and  oriented.  These  investigators  recognized  the  shortcomings  of  this 
assumption  but,  since  they  lacked  a  suitable  alternative,  the  approximation  was 
reasonable.  This  paper  offers  a  mathematical  description  of  the  anisotropic  nature  of 
sea  ice  which  is  suitable  for  incorporation  into  ice  dynamics  and  thermodynamics 
models. 

In  the  section  below  we  describe  the  development  of  the  model.  Plans  for 
validating  the  model  and  a  summary  of  how  the  model  can  be  used  in  an  ice 
dynamics/thermodynamics  model  are  described  in  following  sections. 

MODEL  DEVELOPMENT 

Our  motivation  for  describing  the  ice  behavior  more  accurately  is  so  that  climate, 
ice  dynamics,  and  ice-generated  noise  can  be  better  understood  and  forecasted.  Existing 
models  of  these  features  are  of  the  continuum  type.  Therefore,  we  elected  to  use  a 
continuum  approach  for  our  constitutive  model  so  that  it  will  readily  adapt  to  existing 
models. 


Overview 

The  behavior  of  the  pack  ice  element  shown  in  Figure  1  is  controlled  by  the 
behavior  of  the  thick  ice  and  the  lead  ice.  The  plastic  failure  of  the  thick  ice  and  the 
active  leads  is  characterized  by  ice  strength  parameters  written  as  stress  resultants  (the 
integral  of  the  stresses  through  the  ice  thickness)  to  be  independent  of  the  local 
variations  in  ice  thickness. 

The  thick  ice  is  assumed  to  be  plastic  with  a  Coulomb  yield  criterion  and  a 
normal  flow  rule.  Coon  (1974)  used  the  Coulomb  yield  criterion  for  pack  ice  and,  in 
addition,  the  normal  flow  rule  has  been  used  in  many  sea  ice  models  (e.g.,  Coon,  1980; 
Hibler,  1979;  and  Pritchard,  1981).  More  recently,  Hopkins  and  Hibler  (1991)  found  a 
normal  flow  rule  through  numerical  simulations  of  compact  ice  floes.  The  thick  ice  is 
assumed  to  be  homogeneous  and  isotropic.  The  cohesive  strength  and  angle  of  friction 
of  the  thick  ice  are  assumed  to  be  controlled  by  the  flawed  (cracked)  ice  floes.  Another 
view  of  the  Coulomb  criterion  is  that  in  two  dimensions  the  ice  behavior  is  controlled 
by  a  compressive  strength  which  is  much  greater  than  its  tensile  strength.  Here  again, 
the  cracked  floes  produce  this  smaller  tensile  strength.  Lastly,  new  leads  are  formed  is 
this  thick  ice  along  characteristic  directions  which  admit  velocity  discontinuities 
(Pritchard,  1988). 
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Figure  1.  Element  of  pack  ice  with  a  new  lead  forming  in  the  thick  ice  and  a  refrozen 


Figure  2.  Failure  Envelope  of  a  Mohr-Coulomb  Solid  using  a  Mohr  Circle 
The  following  angles  are  used  in  physical  space: 
a-  from  principal  stress  resultant  direction  to  the  new  lead  direction 
S-  from  new  lead  direction  to  the  oblique  angle  of  the  motion-of  one  side  of  the 
new  lead  relative  to  the  other  (the  friction  angle) 

0-  from  x-axis  direction  (Figure  1)  to  principal  stress  resultant  direction 
0-  from  x-axis  direction  (Figure  1)  to  new  lead  direction 
From  the  geometry  in  Figure  2,  we  can  find  the  angle  a  from  the  angle  S. 


Thus,  the  angle  a  is  also  a  material  property.  Now,  the  angle  of  the  new  lead  can  be 
found  from  angle  of  the  principal  direction  and  the  angle  6: 


0  -  0  ± 


(2) 


The  sign  depends  upon  the  sign  of  the  shear  stress  resultant.  Also  from  the  geometry 
in  Figure  2,  the  equation  for  the  line  representing  the  Mohr-Coulomb  yield  surface  is 


sin  li  -  - - -  (2) 

C 

-  _ —  +  _ 

2  tan  6 
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^  The  axis  goes  from  the  tip  through  the  origin  and  the  square  root  in  Equation  (8) 
is  the  radius  of  the  cone.  A  new  stress  resultant  coordinate  system  is  useful  for  viewing 
the  yield  surface.  We  define  the  following  variables: 


P 


a 

X  y 


S  -  a 


*y 


(10) 


where  P  is  the  pressure  resultant,  D  is  the  difference  between  the  normal  stress 
resultants,  and  S  is  the  shear  stress  resultant.  As  illustrated  in  Figure  3,  the  cone  now 
lies  along  the  P  axis  with  D  and  S  in  the  plane  of  the  circle  normal  to  the  axis  of  the 
cone.  The  yield  function  in  P-D-S  stress  resultant  space  now  follows: 


F  -  {d^S2  +  P  sin  &  -  C  cos  6 


(ID 


From  Equation  (6)  it  can  be  seen  that  2<p  is  the  angle  of  the  stress  resultant  in  the  D-S 
plane  from  the  D-axis  towards  the  S-axis. 

The  normal  flow  rule  states  that  the  plastic  strain  rate  must  be  normal  to  the 
yield  surface  at  a  smooth  point  or,  at  a  comer,  must  lie  between  the  outward  normals 
to  the  adjacent  yield  surfaces.  Therefore,  the  mechanical  properties  of  a  lead  which 
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plane 


Figure  3.  Conical  Yield  Surface  of  Thick  Ice 


Figure  4.  Single  Lead  in  the  Element 
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Figure  5 •  Yield  Surface  for  a  Lead 


Table  1.  Lead  Strain  Rates 


Strain  Rate 

Tension 

Compression 

Shear 

** 

-  (1-cos  20) 

2 

--  (1-cos  20) 

2 

-A  sin  20 

2 

s 

A  (l+cos  29) 

-A  (l+cos  20) 

A  sin  20 

2 

sin  29 

2 

A  sin  29 

2 

A  cos  20 

2 

Table  2.  Example  Material  Properties 


Thick  Ice 


h  =  3  m 


8  -  25° 


C  =  90,000  N/m 


Lead  Ice 


h  =  0.30  m 


O'  =  45,000  N/m 
oc  =  9,000  N/m 
os  =  6,000  N/m 


Figure  6a  shows  the  lead  aligned  with  the  x-axis.  For  this  case,  9  =  0,  so  the  equations 
for  the  yield  surface  are  as  follows. 

For  the  thick  ice: 


+  P  sin  25°  -  9000  cos  25° 


(18) 
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For  shearing  the  lead  (positive  shear  has  the  far  side  of  the  lead  moving  to  the 
right): 


5  -  ±6000  (21> 

The  equation  for  thick  ice  failure  forms  a  right  circular  cone  with  its  axis  along  the  P- 
axis.  The  lead  failure  is  a  long  rectangular  prism  with  two  opposite  sides  representing 
shear  failure  and  the  other  two  sides  representing  opening  and  closing  of  the  lead. 
When  the  two  yield  surfaces  are  combined,  the  rectangular  prism  intersects  with  the 
cone  to  form  a  rectangular  prism  with  rounded  ends.  Since  the  rounded  ends  came  from 
the  cone,  these  rounded  ends  represent  the  only  portion  of  the  yield  surface  which  can 
form  new  leads. 

Figure  6  also  shows  three  views  of  the  combined  yield  surface:  Figure  6b  is  the 
plane  for  P  equal  to  zero,  Figure  6c  is  for  D  equal  to  zero,  and  Figure  6d  is  for  S  equal 
to  zero.  In  the  P=0  and  D=0  views,  the  thick  ice,  represented  by  the  dashed  lines,  is 
not  involved  in  the  ice  failure.  As  seen  in  Figure  6d,  the  rectangular  prism  associated 
with  the  lead  makes  an  angle  of  45  °  with  the  P  axis.  The  effect  of  the  lead  orientation 
angle  9  is  to  rotate  the  lead’s  prism  around  the  P-axis  while  maintaining  a  constant 
angle  of  45°  between  the  prism  axis  and  the  P-axis.  The  physical  deformations 
associated  with  different  faces  of  the  yield  surface  have  been  noted  on  the  figure.  The 
plastic  strain  rate  vectors  associated  with  each  surface  is  the  outward  normal,  except  at 
the  intersections.  For  example,  the  intersection  of  ’lead  closing’  and  ’new  lead  forming’ 
in  Figure  6d  shows  a  fan  for  the  outward  normal,  which  implies  that  both  motions  can 
act  simultaneously  and  their  proportion  is  determined  by  the  given  problem. 

It  is  interesting  to  note  that  a  yield  surface  formed  by  two  leads  excludes  the 
contribution  from  the  thick  ice,  i.e.,  the  origin  is  completely  enclosed  by  the  intersection 
of  the  two  prisms.  Consequently,  a  third  lead  cannot  be  created.  Therefore,  it  appears 
that  there  are  only  three  types  of  yield  surfaces:  thick  ice  with  no  leads,  thick  ice  with 
one  lead,  and  thick  ice  with  two  leads.  Summer  ice  conditions  with  open  water 
surrounding  all  the  floes  is  represented  by  the  yield  surface  collapsed  to  the  origin,  i.e., 
the  ice  has  no  strength. 
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Very  High  Resolution  Radiometer  (AVHRR)  satellite  imagery,  will  provide  a  view 
of  ice  conditions  in  the  study  area. 

APPLICATION  TO  ICE  DYNAMICS  MODELS 

To  incorporate  this  anisotropic  pack  ice  constitutive  model  into  other  models, 
such  as  ice  dynamics  models,  a  procedure  was  developed  to  update  the  yield  surface  as 
a  result  of  ice  deformations  and  thermodynamics.  Based  on  the  strain  rates  found  by 
the  ice  dynamics  model  at  each  time  step,  the  change  in  lead  area  is  calculated.  The 
change  in  lead  area  causes  changes  in  the  distribution  of  ice  thickness  within  the  lead. 
The  lead  ice  thickness  also  changes  due  to  thermodynamics.  The  thinnest  ice  in  the 
updated  lead  provides  the  basis  for  estimating  an  updated  set  of  strength  parameters. 

The  distribution  of  ice  thickness  must  be  calculated  for  each  lead  system  to 
determine  its  ice  strength  parameters.  It  is  anticipated  that  satellite  data  will  be  used 
to  establish  and  correct  these  oriented  lead  ice  thickness  distributions  when  available. 
There  is  also  an  ice  thickness  associated  with  the  thick  ice  to  establish  its  strength 
parameters. 

The  rate  of  area  change  for  each  lead  system  is  calculated  from  the  element 
strain  rates.  Since  all  the  ice  in  the  element  is  assumed  rigid  except  for  the  active  lead, 
all  energy  dissipation  due  to  plastic  yielding  of  the  ice  occurs  in  the  active  lead,  and  the 
strain  rate  parallel  to  the  active  lead  is  zero.  The  shearing  and  normal  velocity 
components  of  one  side  of  the  lead  relative  to  the  other,  as  shown  in  Figure  4,  are 
derived  to  dissipate  energy  at  the  same  rate  as  in  the  whole  element,  with  the  result 
that 

4.—  -  -i  cot  0  +  e  tan  0 
ab  y  (22) 


where  ab  is  the  element  area,  $U  is  the  shearing  velocity  jump  across  the  lead,  and  <SV 
is  the  normal  velocity  jump  across  the  lead  (opening  velocity). 

DISCUSSION 

It  has  been  shown  that  a  plastic  model  of  the  behavior  of  pack  ice  can  be 
developed  which  accounts  the  creation  of  new  leads  and  the  strength  and  orientation 
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Abstract 

An  initial  field  test  was  conducted  to  determine  if  the  stress  state  of  sea  ice  has  an  effect  on  its  microwave  signature.  The 
tests  were  conducted  during  March  and  April  of  1992  in  the  Beaufort  Sea.  An  inflatable  jack  was  placed  into  a  vertical  slot 
in  the  ice  to  induce  stress  in  the  ice  while  a  microwave  transmitter/receiver  monitored  the  backscatter  from  the  ice. 
Measurements  were  made  at  X-Band  and  C-Band  frequencies,  with  various  polarization  combinations,  on  1.3  m  thick 
first-year  ice  and  on  a  freshwater  melt  pond  within  a  multiyear  floe.  Changes  in  backscatter  were  seen  in  both  the  first-year 
ice  and  the  melt  pond  ice  at  both  X-Band  and  C-Band  frequencies,  but  not  at  all  polarizations  and  orientations.  In  general, 
the  largest  change  in  microwave  response  took  place  when  the  ice  was  pressurized  for  the  first  time.  The  measured  changes 
in  microwave  response  were  less  than  2  dB.  However,  after  accounting  for  the  mismatch  between  the  radar  footprint  and  the 
size  of  the  area  su*essed.  the  effect  on  the  backscatter  signal  was  13  dB.  Response  was  seen  at  jack  pressures  at  low  as  about 
70  kPa,  which  resulted  in  an  average  stress  over  the  radar  footprint  of  about  13  kPa.  In  comparison,  natural  pack  ice  stresses 
range  up  to  400  kPa. 

Keywords:  Microwave  response;  Stress  in  sea  ice:  Sea  ice;  Initial  field  tests 


1.  Introduction 

The  effect  of  sea  ice  properties  on  active  mi¬ 
crowave  imaging  of  sea  ice  has  been  the  subject  of 
numerous  studies,  covering  a  period  of  more  than 
twenty  years.  It  is  generally  accepted  that  the  param¬ 
eters  that  have  been  investigated  to  date  do  not 
explain  all  of  the  variation  in  sea  ice  microwave 
signature.  A  parameter  that  has  not  previously  been 
investigated  is  the  stress  state  of  the  ice.  From 
previous  studies  (e.g..  Coon  et  al.?  1989,  and  Lau 
and  Knoke,  1991),  it  is  known  that  the  stress  state  of 
pack  ice  can  vary  markedly  over  periods  of  an  hour 


or  less.  Increased  stress  will  cause  microcracking 
within  the  ice  which  should  have  an  effect  on  the 
microwave  signature.  Additionally,  in  thin,  warm 
ice.  stress  is  expected  to  expel  brine  onto  the  ice 
surface  which  may  also  have  a  significant  effect  on 
the  microwave  signature. 

Considerable  work  has  been  done  to  document 
stress-induced  microcracking  in  the  laboratory  (e.g., 
Sinha,  1982:  Hallam  et  al./l986;  and  Cole/ 1986). 
However,  we  are  aware  of  no  previous  work  investi¬ 
gating  the  effect  of  ice  stress  state  and/or  microc¬ 
racking  on  microwave  signature. 

This  is  a  report  of  an  initial  series  of  field  tests  to 
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investigate  the  effect  of  ice  stress  on  the  microwave 
signature  of  sea  ice.  Demonstrating  and  quantifying 
such  an  effect  may  lead  to  a  better  understanding  and 
interpretation  of  radar  imagery  of  sea  ice.  and  may 
allow  changes  in  ice  stress  to  be  sensed  remotely. 

2.  Equipment  and  methodology 

The  microwave  system  used  in  this  work  was 
developed  at  The  Environmental  Research  Institute 
of  Michigan  (ERIM)  and  has  been  used  on  numerous 
sea  ice  investigations.  It  consists  of  a  frequency- 
modulated,  continuous-wave  (FM-CW)  radar  and 
operates  at  wavelengths  of  3  cm  (X-Band)  and  6  cm 
(C-Band).  There  are  four  transmit- receive  polariza¬ 
tions  possible.  These  include  VV,  HH,  HV  and  VH 
(where  the  V  represents  vertical  and  H  represents 
horizontal).  Ground  illumination  area  was  about  1  m 
X  l.l  m  at  the  6cm  wavelength  and  0.5  m  X  0.6  m  at 
the  3  cm  wavelength.  Additional  details  of  the  sys¬ 
tem  may  be  found  in  Onstott  and  Shuchman  (1988). 

An  inflatable  jack  (a  stiffened  rubber  air  bladder) 
made  by  OTC  Power  Team,  Owatonna,  MN,  was 
used  to  pressurize  the  ice.  This  device  measured 
38  X  38  X  2.5  cm  and  could  be  pressurized  to  800 
kPa.  The  jack  was  pressurized  with  a  small  air  pump 


powered  by  a  12  VDC  battery.  A  photograph  of  the 
test  configuration  is  shown  in  Fig.  1 . 

The  testing  was  done  during  the  LEADEX  field 
program  in  March  and  April  of  1992  in  the  Southern 
Beaufort  Sea.  Measurements  were  made  in  two  ice 
types:  1.3  m  thick,  first-year  sea  ice.  and  in  a  re¬ 
frozen  melt  pond  ice  on  a  multiyear  floe.  A  core  was 
taken  of  the  upper  70  cm  of  the  first-year  ice.  The 
eight  samples  analyzed  ranged  from  3.3  ppt  to  6.6 
ppt  with  an  average  of  5.2  ppt.  There  was  no  signifi¬ 
cant  trend  with  depth.  An  ice  core  collected  from  the 
melt  pond  ice  had  salinities  ranging  from  0.0  to  1.0 
ppt  to  a  depth  of  57  cm  with  an  average  of  0.55  ppt 
for  the  seven  samples.  Here  also  there  was  no  signif¬ 
icant  trend  with  depth.  From  the  surface  to  a  depth  of 
57  cm  the  ice  was  clear.  The  ice  below  this  depth 
had  a  milky  appearance,  marking  the  bottom  of  the 
pond.  At  62  and  72  cm  depths  the  salinity  measured 
0.65  and  0.60  respectively. 

The  general  sequence  of  the  tests  was  as  follows: 

1 .  Snow  cover  was  removed  from  the  test  area. 

2.  A  chainsaw  was  used  to  cut  a  slot  in  the  sea  ice. 
The  slot  measured  approximately  40  cm  long,  36 
cm  deep,  and  3  cm  wide.  After  the  initial  test, 
slots  were  cut  in  the  ice  at  the  ends  of,  and 
normal  to,  the  3  cm  slot.  These  extended  approxi¬ 
mately  60  cm  from  the  main  slot.  These  were  cut 


Fig.  1.  The  microwave  system  and  inflatable  jack  installed  in  the  ice  during  test  2  in  first-year  ice. 
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Fig.  2.  The  test  configuration  from  the  top  and  side  showing  the 
air  jack  and  the  side  slots. 


to  simplify  the  calculation  of  stress  in  the  region 
around  the  jack.  The  ice  loading  configuration, 
with  the  side  slots  cut,  is  shown  in  Fig.  2. 

3.  An  inflatable  jack  was  inserted  into  the  slot  and 
was  connected  to  a  manifold  with  a  dial  gauge, 
which  in  turn  was  connected  to  an  air  pump.  The 
top  of  the  jack  extended  above  the  ice  surface  less 
than  2  cm. 

4.  The  microwave  transmitter/receiver  was  aimed 
at  the  area  around  the  inflatable  jack  and  was  set 
to  the  desired  polarity. 

5.  The  microwave  system  was  turned  on. 

6.  The  inflatable  jack  was  pressurized  to  the  desired 
level,  and  the  pressure  and  the  microwave  re¬ 
sponse  were  manually  recorded. 

7.  The  jack  pressure  was  changed,  and  new  readings 
were  taken.  This  process  was  repeated  until  tests 
were  done  at  all  the  desired  pressures. 


8.  Tests  were  done  with  the  jack  both  parallel  and 
perpendicular  to  the  axis  of  the  radar. 


3.  Results 

The  testing  was  done  both  at  X-Band  and  C-Band 
microwave  frequencies.  The  greatest  response  was 
seen  at  C-Band  frequencies,  and,  for  this  reason, 
only  data  from  this  frequency  are  shown  in  the 
figures.  All  of  the  tests  were  done  at  a  30°  incidence 
angle.  The  results  of  two  of  these  tests  are  summa¬ 
rized  below. 

Test  2  was  done  in  first-year  ice.  The  axis  of  the 
radar  was  parallel  with  the  inflatable  jack  slot,  and 
stress  relief  slots  were  present.  This  was  the  second 
time  this  ice  had  been  pressurized.  As  can  be  seen  in 
Fig.  3,  the  VV  polarization  response  generally  fol¬ 
lowed  the  air  pressure  in  the  jack.  The  HH  response 
was  not  as  large. 

Test  5  was  done  in  melt  pond  ice.  The  melt  pond 
ice  had  many  large  cracks,  however  the  ice  within 
the  radar  footprint  had  no  large  cracks.  The  radar 
axis  was  perpendicular  to  the  inflatable  jack  slot.  Ice 
at  this  test  location  had  not  previously  been  pressur- 


o  5  10  is  20 
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Fig.  3.  The  backscatter  response  shown  is  for  test  2.  on  first-year 
ice.  The  radar  is  oriented  with  its  axis  parallel  to  the  pressure  slot. 
The  backscatter  power  scale  is  shown  on  the  left  vertical  axis 
while  jack  pressure  is  shown  on  the  right.  Jack  pressure  was 
increased  from  0  to  690  kPa  and  then  returned  to  0  kPa.  Backscat¬ 
ter  for  both  VV  and  HH  polarizations  are  shown. 
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Fig.  4.  The  backscatrer  response  shown  is  for  test  5,  on  melt  pond 
ice.  The  radar  is  oriented  with  its  axis  perpendicular  to  the 
pressure  slot.  Relief  slots  were  present  and  ran  perpendicular  to 
the  pressure  slot.  Jack  pressure  was  increased  from  0  to  480  kPa, 
the  pressure  was  relieved  and  then  was  increased  to  690  kPa. 
Responses  at  VV,  HH,  and  VH  polarizations  are  shown. 

ized.  Fig.  4  shows  the  results  for  HH,  VV,  and  HV 
polarizations.  Here  changes  were  seen  in  ail  three 
polarizations  with  the  largest  response  in  HV, 

In  general,  it  appeared  that  the  largest  change  in 
microwave  response  took  place  when  the  ice  was 
pressurized  for  the  first  time.  Subsequent  tests  on  the 
same  ice  generally  did  not  have  as  large  a  change  in 
backscatter.  It  appears  the  response  is  influenced  by 
loading  path  and  loading  history. 

It  is  clear  from  the  data  that  something  within  the 
radar  footprint  changed  during  pressurization.  We 
considered  two  possible  causes  not  directly  related  to 
the  ice  stress  state.  The  first  was  the  size  change  of 
the  exposed  upper  part  of  the  inflatable  jack  during 
pressurization.  The  exposed  part  was  less  than  two 
cm  tall.  This  size  change  was  small,  and  after  the 
jack  contained  a  nominal  pressure,  on  the  order  of  a 
few  kPa.  virtually  no  additional  size  change  oc¬ 
curred.  For  this  reason  we  believe  this  size  change 
had  little  or  no  effect.  The  second  alternative  cause 
considered  was  that  the  slope  of  the  ice  near  the 
inflatable  jack  may  have  changed  an  imperceptible 
amount  during  pressurization.  However,  if  this  did 
occur,  we  believe  the  microwave  system  would  not 
be  sensitive  to  such  a  small  change. 

The  microwave  response  changes  were  small  to 
moderate.  For  tests  2  and  5,  the  maximum  changes 
are  1.6  dB,  and  1.7  dB.  respectively.  However,  the 


radar  footprint  and  the  size  of  the  area  that  was 
stressed  were  not  well  matched.  The  radar  footprint 
is  approximately  2.6  times  larger,  and  the  effect  on 
the  backscatter  signal  should  be  about  7.6  times 
greater  due  to  the  overall  difference  in  area  and 
volume  observed  for  the  assumption  of  total  beam 
filling.  This  calculates  out  to  effective  changes  of  12 
dB  and  13  dB  for  tests  2  and  5  at  air  jack  pressures 
of  690  kPa. 

A  calculation  of  the  stress  in  the  ice  was  made 
based  on  the  solution  of  an  elastically  loaded  90° 
wedge  as  described  by  Timoshenko  and  Goodier 
(1951).  This  configuration  is  that  shown  in  Fig.  2. 
The  results  show  that  there  is  essentially  no  stress  in 
the  ice  at  a  distance  greater  than  38  cm  from  the  air 
jack.  (For  the  configuration  without  side  slots  the 
stress  zero  point  would  be  less  than  38  cm  from  the 
jack.)  Therefore,  the  average  stress  in  the  38  cm 
square  region  would  be  approximately  half  of  the 
pressure  in  the  air  jack.  For  example,  the  maximum 
test  pressure  in  the  jack,  690  kPa,  produced  an 
average  pressure  of  350  kPa  over  the  38  cm  square 
area.  Dividing  this  pressure  by  the  ratio  of  the  area 
of  the  radar  footprint  to  the  area  stressed,  2.6,  results 
in  an  average  stress  in  the  radar  footprint  of  135  kPa. 
A  response  was  seen  at  jack  pressures  as  low  as  70 
kPa  which  results  in  an  average  stress  in  the  radar 
footprint  of  13  kPa. 

How  do  the  above  stress  values  compare  with 
those  in  natural  sea  ice?  Pack  ice  stresses  were 
measured  up  to  275  kPa  by  Coon  et  al.  (1989)  and  to 
400  kPa  by  Perovich  et  al.  (1992)  and  Richter-Menge 
et  al.  (1995).  Therefore,  the  response  seen  was  at 
stresses  well  within  the  range  of  natural  stresses  in 
sea  ice. 


4.  Discussion 

Something  within  the  radar  footprint  changed 
when  the  ice  was  stressed,  which  markedly  altered 
the  radar  signature  of  the  ice.  This  took  place  for 
both  the  first-year  ice  and  the  relatively  fresh  melt- 
pond  ice.  We  suspect  that  the  change  took  place  due 
to  microcracking  or  deformation  of  existing  microc¬ 
racks  within  the  ice  or  perhaps  through  deformation 
of  air  bubbles  within  the  ice.  For  the  first-year  ice, 
other  possible  mechanisms  include  distortion  of  the 
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brine  channels  or  brine  pockets  and  expulsion  of 
brine  onto  the  surface  of  the  ice. 

More  detailed  studies  should  be  conducted  to 
better  quantify  and  characterize  the  effect  found  in 
this  preliminary  work.  We  believe  the  most  accurate 
and  reliable  results  will  be  obtained  through  testing 
of  natural  sea  ice  under  arctic  conditions.  An  im¬ 
proved  stress  fixture  should  be  used  to  provide  a 
more  uniform  stress  field  and  to  better  match  the 
radar  footprint  with  the  stressed  area.  Variables 
should  include  radar  frequency,  polarizations,  orien¬ 
tation.  look  angle,  number  of  pressure  cycles,  and  ice 
thickness  and  ice  type.  Once  the  relationship  be¬ 
tween  ice  stress  and  microwave  response  is  better 
defined,  it  would  be  appropriate  to  conduct  analyses 
of  the  physical  changes  in  the  ice  that  cause  the 
change  in  microwave  properties.  This  may  best  be 
done  as  a  combination  of  laboratory  and  polar  field 
studies. 
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ABSTRACT 

This  paper  presents  an  Ice  behavior  model  based  on  the  study  of 
columnar  ice  deformation  phenomena  from  many  experimental 
configurations.  The  important  and  significant  feature  of  this  plastic 
constitutive  law  is  its  reliance  on  slip  planes,  analogous  to  those  in  soils  and 
within  crystals.  Thus,  deformations  which  are  solutions  to  the  boundary 
value  problems  may  result  from  continuous  or  discontinuous  deformations 
on  these  planes.  The  analysis  of  any  ice  mechanics  problem  involves  the 
stresses  and  strains  in  a  coordinate  frame  tied  to  the  ice  fabric,  e.g.,  the 
dominant  grain  direction  and,  if  the  c-axes  are  aligned,  the  c-axis  direction. 
In  this  paper,  we  will  present  the  3-D  formulation  and  consider  how  the 
plastic  constitutive  model  might  be  applied  to  such  apparently  simple  cases 
as  uniaxial  and  biaxial  compression  test  samples. 

KEY  WORDS:  ice,  columnar,  grain,  constitutive,  basal,  stress,  strain 


INTRODUCTION 

The  objective  of  this  paper  is  to  develop  a  constitutive  law  For 
columnar  sea  ice  suitable  for  solving  boundary  value  problems  such  as  sea  ice 
loads  on  structures.  Such  a  constitutive  law  must  include  recoverable  and 
non-recoverable  deformations.  To  carefully  model  the  behavior  of  the  sea  ice, 
the  law  must  also  be  anisotropic  under  multi-axial  loading  conditions. 
Assuming  that  the  ice  sheet  is  a  horizontal  sheet  with  vertical  columnar  grains, 
we  have  modeled  the  anisotropic  non-recoverable  deformations  as  four  shear 
mechanisms  on  three  types  of  slip  planes:  a  vertical  plane  with  horizontal  slip 
along  grain  boundaries,  an  angled  plane  which  shears  across  grains,  the  basal 
plane,  and  a  vertical  plane  with  vertical  slip  along  grain  boundaries.  The 
constitutive  model  is,  however,  applicable  only  to  continuous  ice  pieces  with 
no  cracks  and  for  strain  rates  in  the  ductile  range  (10*6  to  10~3  1/s). 

Columnar  sea  ice  consists  of  parallel  ice  grains;  mechanically  these 
grains  behave  as  though  they  were  a  stack  of  rods.  The  rods  slide  against 
each  other  easier  than  they  can  be  sheared;  therefore  the  ice  fails  at  a  lower 
stress  in  situations  where  the  rods  can  slide  against  each  other  than  when  the 
rods  are  forced  to  be  sheared  off.  Furthermore,  since  the  ice  has  basal  planes 
that  are  aligned  within  the  grains,  the  ice  can  fail  at  an  even  lower  stress  in 
situations  that  the  basal  planes  can  slip,  i.e.,  when  an  ice  block  has  all  its  basal 
planes  aligned. 

Columnar  ice  exhibits  four  distinct  shear  mechanisms,  as  illustrated 
in  Figure  1.  In  this  figure,  we  consider  the  x-axis  to  be  normal  to  the  grain 
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Figure  1.  Shear  Deformation  Mechanisms 


direction  and  normal  to  the  basal  planes  (parallel  to  the  c-axis  within  the 
grain),  the  y-axis  also  be  normal  to  the  grain  direction  but  in  the  basal  plane, 
and  the  z-axis  be  parallel  to  the  grain  direction  and  vertical.  Basal  plane  slip 
and  out-of-plane  grain  slip  do  not  occur  if  the  principal  stress  directions  arc 
aligned  with  these  x-y-z  axes.  The  grain  direction  in  Figure  1  is  indicated  in 
the  lower  left  corner  of  each  figure;  the  hatching  approximates  the  directions 
of  the  slip  planes  for  each  mechanism.  Basal  plane  slip  and  out-of-planc  grain 
slip  result  in  shear  strains  in  the  principal  stress  directions  as  indicated  by  the 
dashed  rectangles;  in-plane  grain  slip  and  grain  shear  do  not. 

1.  In-Flane  Grain  Slip  Shear  stress  in  the  plane  of  the  ice  sheet  causes 
the  grains  to  slide  past  each  other  in  a  direction  perpendicular  to  the  grain 
direction.  A  uniaxial  compression  test  in  the  plane  of  the  ice  sheet  would  be 
effective  in  activating  this  mechanism. 

2.  Grain  Shear  The  slip  plane  for  grain  shear  cuts  across  grain 
boundaries  and  basal  planes,  at  an  angle  to  the  columnar  grain  direction.  As 
such,  grain  shear  is  the  mechanism  with  the  highest  strength.  This  type  of 
out-of-plane  deformation  occurs  when  one  of  the  principal  stress  directions  is 
perpendicular  to  the  ice  sheet.  A  uniaxial  compression  test  normal  to  the 
plane  of  the  ice  sheet  would  be  effective  in  generating  this  mechanism. 

3.  Basal  Plane  Slip  If  the  basal  planes  in  adjacent  grains  are  parallel, 
the  c-axes  are  aligned,  then  shear  stress  on  the  basal  plane  will  cause  the  basal 
planes  to  slip  past  each  other,  A  uniaxial  compression  test  at  45°  to  the  basal 
plane  would  be  effective  in  activating  this  mechanism.  Shear  stress  causes 
shearing  motion  in  the  plane  of  the  ice  sheet  and  shear  stress  aa  causes 
shearing  motion  out  of  the  plane.  We  have  assumed  that  the  response  is  the 
same  in  cither  direction. 

4.  Out-of-Plane  Grain  Slip  Shear  stress  parallel  to  the  grain  direction 
causes  the  ice  to  slide  along  the  grain  boundaries  parallel  to  the  grains.  Since 
the  ice  grains  run  vertical  to  the  plane  of  the  ice  sheet,  this  type  of  out-of- 
planc  deformation  requires  shear  stress  out  of  the  plane,  i.e.,  a ^  or  a  .  A 
uniaxial  compression  test  at  60°  to  the  plane  of  the  ice  sheet  would  be 
effective  in  generating  this  mechanism.  Since  we  expect  the  behavior  of  this 
mechanism  to  be  significantly  different  from  in-plane  slip,  we  have  modeled 
this  as  a  distinct  mechanism,  although  they  both  involve  slip  between  grains. 


The  analysis  of  an  ice  load  problem  can  be  done  with  one  coordinate 
frame  aligned  with  the  ice  fabric  and  another  coordinate  frame  aligned  with 
the  problem.  Having  defined  the  two  coordinate  frames,  the  six  stress 
components  in  the  coordinate  frame  of  the  ice  fabric  can  be  determined  from 
the  stress  state  in  the  problem  coordinate  frame  following  Meilor,  1983.  The 
same  transformation  can  be  applied  to  strains  and  strain  rates. 

The  next  section  defines  the  structure  of  the  plastic  constitutive  law 
for  columnar  sea  ice  in  terms  of  the  four  shear  mechanisms.  In  the  third 
section,  we  provide  detail,  material  properties,  and  examples  for  the  in-plane 
grain  slip  and  grain  shear  mechanisms:  in  the  fourth  section,  the  basal  plane 
slip  and  out-of-plane  grain  slip  mechanisms  are  considered  and  an  example 
provided. 

PLASTIC  CONSTITUTIVE  LAW 

The  phenomenological  plastic  constitutive  law  developed  in  this  paper 
combines  many  observations  of  columnar  ice  behavior  and  provides  for 
multiaxial  stress  and  deformation.  The  concept  is  that  all  nonrecoverable 
deformation  of  columnar  ice  occurs  on  discrete  planes.  The  first  assumption 
made  here  is  that  the  total  strain  rate  can  be  separated  into  a  recoverable 
strain  rate  £rj-  and  a  nonrecoverable  strain  rate  expressed  by: 

i’v  -  -  *lm  oa  *  U> 


where  the  recoverable  strain  rate  can  be  related  to  the  stress  rate  through  the 
elastic  properties  A^,.  The  object  of  this  paper  is  to  relate  the 
nonrecoverable  part  of  the  strain  rate  to  the  stress  and  determine  the  yield 
surface  which  encloses  the  origin  in  stress  space.  Yield  implies  ductile  failure 
in  plasticity  theory,  although  most  laboratory  tests  measure  the  failure  stress 
of  ice  as  a  peak  stress  or  a  sample  fracture. 


It  is  assumed  that  deformations  on  any  of  the  yield  surfaces  are 
controlled  by  the  stress  st3te.  Therefore,  there  is  a  yield  function  in  stress 
space  which  serves  as  a  plastic  flow  potential.  At  this  point  it  is  valuable  to 
define  the  plastic  constitutive  law  and  boundary  value  problem  solution 
procedures.  For  ice  it  is  assumed  that  there  exists  a  function  F: 

F  i  ci<**  Ci&r  O-  C3  -  3  -o  0) 

The  function  F  depends  on  the  stress  Oj:  and  the  material  properties  C2, 
C3,  ...  which  in  turn  depend  upon  tne  nonrecoverable  strain  rate  and 
properties  of  the  ice  aa  such  as  temperature,  brine  volume,  and  grain  size. 

The  function  F  in  equation  (2)  has  been  treated  as  a  yield  surface  by 
Coon  et  al.  (1984).  In  that  view,  the  yield  surface  depends  on  the  strain  rate, 
the  form  of  F  accounts  for  the  anisotropy  of  the  ice,  and  the  yield  surface  is 
also  the  plastic  flow  potential.  We  assume  that  the  ice  is  a  stable  material, 
such  that  work  cannot  be  extracted  from  plastic  deformations  in  any  situation. 
It  follows  from  this  assumption  that  the  yield  surface  of  the  material  is  convex 
in  stress  space  and  that  any  plastic  strain  increment  vector  superposed  on  the 
stress  space  with  corresponding  axes  matching  would  lie  in  the  direction  of  the 
outward  normal  to  the  yield  surface  or  between  outward  normals  at  a  corner. 
Thus,  the  nonrecoverable  strain  rates  are  directly  related  to  the  plastic  flow 
potential  as  follows. 


A  i  +  _  x  dF 
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where  X  is  a  positive  constant  to  be  determined  as  part  of  the  boundary  value 
problem  solution.  The  six  strain  rates  are  related  to  the  three  velocities  v.  by 
linear  equations  in  equation  (3).  The  equations  to  be  solved  are  the  three 
equilibrium  equations,  the  six  velocity  flow  rule  equations,  and  the  yield 
surface;  the  variables  being  solved  for  are  the  six  stresses,  three  velocities,  and 
X.  In  a  solution  procedure,  the  yield  surface  remains  constant  for  a  time  step 
and  the  plastic  increment  of  the  deformation  is  found.  The  yield  surface  is 
then  changed,  because  it  depends  on  the  strain  rates.  The  solution  then 
proceeds  step  by  step.  The  biggest  advantage  of  such  a  formulation  is  that  all 
of  the  procedures  of  plastic  limit  analysis  can  be  applied  to  problems  to  obtain 
the  limit  load,  which,  for  many  problems,  is  of  great  importance.  Large 
computer  codes  are  being  developed  to  determine  ice-structure  interactions 
(e.g.,  Kama,  1992);  the  constitutive  law  presented  here  could  be  used  in  such 
codes.  Ice  anisotropy  and  strain  rate  dependencies  are  very  important  to  the 
failure  loads  of  ice  on  a  vertical  structure,  as  shown  by  Coon  et  al.  (1984). 

The  ductile  behavior  of  ice  is  dependent  on  the  normal  stress  on  the 
slip  plane  because  the  uniaxial  tensile  strength  is  less  than  the  uniaxial 
compressive  strength.  From  test  data,  it  is  dear  that,  when  the  normal  stress 
is  somewhat  larger  than  the  compressive  strength,  this  "pressure*  sensitivity 
disappears  (Coon  et  al.,  1984).  It  may  return  for  very  large  pressures  or  for 
very  slow  (recrystallization)  flow;  we,  however,  are  not  concerned  with  such 
problems  here. 

We  now  consider  the  yield  surfaces  for  the  first  two  mechanisms,  in- 
plane  grain  slip  and  grain  shear,  to  be  those  of  the  Coulomb  and  Tresca  yield 
criteria,  after  Coon  et  al.  (1984).  The  Coulomb  yield  criterion 
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states  that  the  shear  stress  S  on  any  surface  must  not  exceed  the  shear 
strength  of  the  surface.  The  shear  strength,  in  turn,  depends  linearly  upon  the 
normal  stress  N.  The  material  constants,  C  and  fi,  are  called  the  cohesion 
and  the  angle  of  internal  friction.  Shield  (1955)  extended  Coulomb’s  criterion 
to  isotropic  materials  subjected  to  a  three  dimensional  state  of  stress,  by  a 
literal  three  dimensional  interpretation  of  equation  (4).  If  the  material  is 
isotropic,  the  orientation  of  the  stress  state  has  no  relevance,  and  the  stress 
state  is  completely  defined  for  our  purposes  by  Us  principal  values  aj,  o2,  and 
a3.  For  definitiveness,  take  <jx  >  a2  >  ay  Then  the  possible  values  of  S  and 
N  lie  inside  the  at-a3  Mohr’s  circle  and  outside  the  at-a2  circle  and  the  c?2-a3 
circle  (see  Figure  2).  On  the  other  hand,  equation  (4)  plots  as  a  pair  of 
straight  lines  in  N-S  space,  so  the  three  dimensional  extension  of  Coulomb 
simply  states  that  the  aj-cr3  Mohr’s  circle  must  lie  inside  the  wedge  shaped 
region  defined  by  equation  (4). 


Figure  2.  Mohr's  Representation  of  Stress  and  the  Coulomb  Yield  Surface 


Yield  occurs  when  the  wedge  is  tangent  to  the  largest  Mohr's  circle, 
i.e.,  when 

cfj  -  <ij  -  2  C  cos (3  -  (a{  +  Oj)  sinfJ  (5) 

The  material  constants  C  and  J3  may  be  calculated  from  the  uniaxial  tensile 
and  compressive  strengths  aT  and  oc.  In  this  paper,  we  shall  treat  and  <jc 
as  positive  numbers  but  retain  the  usual  compression  negative  convention  for 
calculated  stresses.  Then,  application  of  equation  (5)  to  uniaxial  tension  and 
compression  tests  leads  to 


GRAIN  SHEAR  AND  IN-PLANE  GRAIN  SLIP  MECHANISMS 

To  extend  this  model  to  anisotropic  materials,  the  material  properties 
C  and  13  in  equation  (4)  must  become  functions  of  the  orientation.  There  are 
many  ways  to  do  this,  but  the  columnar  grain  model  suggests  a  particularly 
simple  one.  We  assume  that  there  are  two  sets  of  values  for  C  and  13,  for 
grain  shear  and  in-plane  grain  slip  mechanisms.  First,  we  consider  grain 
shear,  illustrated  in  Figure  lb.  Yield  for  grain  shear  is  governed  by  Cs  and  /3S. 
Consider  a  test  specimen  whose  axis  is  oriented  in  the  z  direction,  subjected 
to  a  uniaxial  state  of  stress  (Figure  lb  illustrates  this  configuration  but  with 
a  small  lateral  load).  If  and  aTs  denote  the  corresponding  compressive 
and  tensile  strengths,  then  from  equation  (6)  we  have 


Next,  consider  an  in-plane  specimen,  as  illustrated  in  Figure  la.  The 
values  Cj  and  I3i  are  the  cohesion  and  angle  of  internal  friction  for  in-plane 
grain  slip.  For  definitiveness,  let  y  be  along  the  axis  of  the  specimen.  Let  ay 
be  non-zero,  and  all  the  other  stress  components  be  zero  (Figure  la  illustrates 
this  configuration  but  with  a  small  lateral  load).  Thus  the  load  is  applied  in 
a  plane  having  less  strength  than  a  vertical  plane,  and,  denoting  the  in-plane 
strengths  as  aCi  and  equation  (6)  yields 


For  the  special  case  where  one  of  the  principal  axes  of  stress  is 
perpendicular  to  the  ice  sheet,  the  equations  of  the  yield  surface  can  be 
derived  from  equations  (5),  (7),  and  (8).  Thus,  for  the  special  case  of  aligned 
axes  of  stress  and  material,  the  yield  surface  is  defined  be  the  following  six 
equations: 


where  the  x-y  coordinates  are  the  principal  axes  of  stress  in  the  plane  of  the 
ice  sheet. 

At  some  point  along  the  negative  N  axis,  the  experimental  evidence 
supports  a  transition  to  a  yield  surface  which  is  independent  of  N,  i.e.,  a 
cylinder  with  generators  parallel  to  the  N  axis  (see  Figure  3,  after  Coon  et  a!., 
1984).  An  anisotropic  Tresca  criterion  is  consistent  with  the  Coulomb 
generalization  which  has  been  used  and  is  indeed  the  limit  of  equation  (4)  as 
tan  f3  ■+  0.  Two  parameters  are  required  to  define  the  yield  surface:  the 
critical  value  of  shear  stress  for  yield  in  the  weaker  x-y  plane  and  the  critical 
value  of  shear  stress  for  yield  in  the  other  planes.  These  critical  values  of 
shear  stress  must  be  measured  in  a  test  conducted  at  high  enough  normal 
stress  on  the  slip  plane  so  that  yield  occurs  on  the  Tresca  portion  of  the  yield 
surface.  Data  from  in-plane  or  out-of-plane  triaxial  tests  provide  the  critical 
shear  strengths  crSj  and  cr^.  When  one  of  the  principal  stress  directions  is 
perpendicular  to  the  ice  sheet,  the  Tresca  portion  of  the  yield  surface  is  giveu 
by  the  following  six  equations: 


Figure  3.  Typical  Set  of  Yield  Surfaces 


Yield  Strength  for  In-Plane  Slip  and  Grain  Shear  Mechanisms 

The  compressive  and  tensile  yield  strengths  of  sea  ice  are  functions 
of  temperature,  grain  size,  brine  volume,  and  strain  rate.  In  Coon  and 
.Associates  (1984),  least  squares  regression  analyses  were  used  to  investigate 
;  these  various  functional  dependencies  to  obtain  empirical  expressions  for  first^ 
year  columnar  sea  ice  based  on  laboratory  and  field  data.  The  form  chosen 
was 


O  -  [  p0  +  p,r  -  *  P,y77  ] 


(11) 


where 


a  is  the  strength  in  MPa, 

T  is  the  ice  temperature  in  °C, 

g  is  the  average  grain  size  of  an  ice  sample  in  mm, 

k  is  the  strain  rate  (from  the  platen  speed)  in  1/s, 

vb  is  the  brine  volume  expressed  as  a  decimal  fraction,  and 

J30,  flj,  fl7,  B3,  and  fl4  are  coefficients  found  from  regression  analysts. 

It  was  assumed  that  equation  (11)  could  be  written  as 


log  a  -  log  p  +  p4  log  k 


The  fundamental  way  of  viewing  the  dependence  of  the  yield 
properties  in  equations  (7)  and  (8)  is  to  find  the  C’s  and  fl’s  as  functions  of 
the  shear  strain  rate  on  the  active  slip  planes.  This  is  not  accomplished  by 
combining  equations  (7),  (8),  (14),  and  (15)  because  then  the  C’s  and  fl’s  will 
depend  on  the  x-y-z  strain  rates.  The  transformation  to  the  active  slip  planes 
requires  the  use  of  the  plastic  flow  potential  function,  the  flow  rule,  and  tensor 
equations  for  stress  and  strain  coordinate  transformations.  This  has  not  been 
done  in  this  paper. 

Example  with  Principal  Stresses  Aligned  with  Ice  Fabric 

Timco  and  Frederking  (1986)  report  test  data  on  columnar  ice  in 
which  the  specimens  are  loaded  in  five  types  of  test  configurations: 


equation  (11)  as 


4"  -  Po  *  f>!r  -  02 VF  *  03 


By  then  rewriting 
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(13) 

D. 

°x  -  o ; 

'  0  '  °«  - 

E. 

a,  -  0  ; 

ay  -  0  ;  (j2  -  -a 

(17) 


A  linear  regression  analysis  was  used  to  determine  the  other  four  coefficients. 
The  following  expressions  were  derived  based  on  the  data  sources  listed  in 
Table  1. 

In-Plane  Slip: 

a  a  -  (65.0  -  0.32T  -  \2Ay/g  -  56.5^  )  kQzt 

an  -  0.71  -  1.74^7.  (14) 

-  0.82  aa 


To  compare  this  test  data  with  the  present  model,  the  material  properties 
must  be  found  from  equations  (14)  and  (15).  The  test  conditions  were  as 
follows: 


k  -  2xl0‘4  s~l 
T  -  -2*C 
Salinity  -  1.65  ppt 

y FI  -  0-2 


(18) 


Grain  Shear: 


a  a  -  (40.5  -  1.027*  -  56.5^  ) 
aTs  -  1.48  -  2.18/77 
“  0.565  o0 


No  grain  size  is  reported  and  therefore  a  grain  size  of  10  mm  has  been 
assumed.  Timco  and  Frederking  (1986)  account  for  porosity  and  brine 
volume,  which  has  not  been  done  here.  Therefore,  the  yield  strengths  in  MPa 
are 

°cr  °n  asi  acr  °Ts  as»  (19) 

1.52  0.37  1.245  4.05  1.04  2.29 


Table  1.  Sources  of  Material  Properties  Dnta 


<*Ti 

*T. 

Frederking  and  Timco  (1981) 

X 

Frederking  and  Timco  ( 1983) 

X 

Hausler  (1981) 

X 

X 

Sinha  (1983) 

X 

Vaudrey  (1977) 

X 

X 

X 

X 

Wang  ( 1979) 

X 

The  only  data  for  tensile  strength  were  Vaudrey’s;  it  was  not  possible 
to  determine  any  variation  in  tensile  strength  with  grain  size  or  strain  rate 
from  the  data.  Also,  there  was  no  grain  size  dependence  found  on  a ^  The 
Tresca  shear  strengths  and  o&  were  found  only  in  the  data  of  Hausler 
(1981).  There  was  insufficient  data  to  determine  all  dependencies;  however, 
the  Hausler  data  has  in-plane  and  out-of-plane  compression  strengths  as  well 
as  test  specimens  subjected  to  multi-axial  stress  states.  Therefore,  it  was 
decided  to  determine  the  ratio  of  the  Tresca  strength  as  to  the  unconfined 
compressive  strength  oc  from  the  data  sets  and  assume  that  the  dependence 
on  temperature,  grain  size,  brine  volume,  and  strain  rate  was  the  same  for  the 
high  normal  stress  in  Tresca  as  for  Coulomb.  This  comparison  was  done  by 
plotting 

\  a  *  -  o  x  \  vs.  \  a  t  +  o  x\ 

and  (16) 

I  oy  -  ox  I  vs. I  ay  +  ax  I 

to  find  the  Tresca  shear  strengths  and,  therefore,  the  ratios  specified  in 
equations  (14)  and  (15). 


Table  2  lists  the  faces  of  the  yield  surface  based  on  equations  (9)  and  (10),  the 
flow  rules  based  on  equation  (3),  and  the  type  of  test  configuration 
intersecting  the  faces  of  the  yield  surface. 

The  Timco  and  Frederking  test  data  are  compared  to  the  yield 
surfaces  from  Table  2  in  Figures  4  and  5.  The  yield  surfaces  from  Table  2 
represent  the  theory  presented  in  this  paper  using  material  properties 
corresponding  with  the  test  conditions  shown  in  equation  (18),  but  none  of  the 
Timco  and  Frederking  test  results  were  used  in  determining  the  yield  surfaces 
in  Table  2.  In  Figure  4,  Type  A  tests  are  on  the  Coulomb  part  of  the  grain 
shear  yield  surface  (Face  1)  but  are  very  near  the  Tresca  part  of  the  yield 
surface  (Face  7).  In  Type  A  tests,  it  is  possible  to  develop  a  small  amount  of 
oz  due  to  friction  from  the  side  constraints  when  the  deformation  is  out  of  the 
plane  of  the  ice.  Figure  4b  shows  that  a  small  amount  of  compression  in  the 
az-direction  could  increase  the  ax  yield  stress  and  that  again  the  test  is  near 
the  Coulomb  and  Tresca  transition  of  Faces  1  and  7. 

Type  B  and  C  tests  are  shown  in  Figures  4a  and  5.  In  Type  B  tests, 
it  is  possible  to  predict  a  small  amount  of  az  which  has  little  effect  on  the 
yield  because  it  moves  the  stress  state  parallel  to  Face  6  of  the  yield  surface. 
It  was  noted  by  Timco  and  Frederking  that  the  confining  platens  had  to  be 
hand-tightened  after  yield  for  Type  B  tests.  This  accounts  for  the  rectangular 
area  of  data  for  Type  B  tests  in  Figure  5.  The  amount  of  az  predicted  is 
small  because  the  deformation  is  in  the  plane  and  therefore  does  not  develop 
much  friction. 

Figure  5  also  shows  the  Type  D  and  E  tests.  There  may  be  scatter 
in  Type  D  and  E  tests  because  the  loading  path  is  again  nearly  parallel  to  the 
yield  surfaces  and  because  the  Coulomb  Face  2  and  Tresca  Face  8  surfaces 
are  nearly  the  same.  In  Type  D  tests,  a  small  amount  of  ox  may  have  been 
introduced  by  friction  with  the  confining  platens  (which  had  to  be  tightened). 
This  is  shown  in  Timco  and  Frederking’s  data  and  will  produce  large  apparent 
fluctuations  in  the  c2  yield  strength. 


The  representation  of  this  data  by  the  present  model  leads  to  insights 
about  the  behavior  of  columnar  sea  ice  by  accounting  for  anisotropic  behavior 
and  "pressure"  saturation.  However,  the  model  also  provides  a  helpful  and 
succinct  way  of  perceiving  the  shear  failure  mechanisms,  ft  is  encouraging  to 
see  such  good  agreement  between  the  behavior  observed  in  these  tests  with 
predictions  based  on  prior  data. 

To  test  the  yield  surface  in  Figures  4  and  5,  it  is  important  to  control 
the  loading  paths.  This  has  been  done  by  Hausler  (1981)  and  by  Schulson  and 
Smith  (1992);  however,  the  strain  rates  in  the  second  case  were  greater  than 
those  considered  in  this  paper.  In  most  investigations,  the  strain  rates  at  yield 
have  not  been  compared  with  the  flow  rule.  These  comparisons  would 
provide  a  much  better  insight  into  the  behavior  of  columnar  sea  ice.  Another 
’  way  of  controlling  the  load  path  is  to  use  a  confining  platen  restrained  by  an 
elastic  spring  in  a  uniaxial  test  machine,  as  reported  by  Coon  et  al.  (1991). 
The  spring  stiffness  provides  the  different  load  paths;  the  test  can  be 
performed  in  a  very  simple  machine. 


Figure  5.  Timco  and  Fredcrking  (1986)  Data  for  Types  8.  C,  D.  and  E  Tests 
in  Shaded  Rectangles  with  Pertinent  Yield  Surfaces 


BASAL  PLANE  SLIP  AND  OUT-OF-PLANE  GRAIN  SLIP  MECHANISMS 


Table  2.  Yield  Surfaces  for  Timco  and  Fredcrking  (t986) 
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We  now  consider  the  other  two  shear  mechanisms:  basal  plane  slip 
and  out-of-plane  grain  slip.  It  is  assumed  that  the  yield  surfaces  of  these 
mechanisms  also  exhibit  Coulomb/Tresca  behavior. 

The  basal  plane  slip  mechanism  occurs  in  ice  with  the  c-axes  aligned 
in  the  plane  of  the  ice  sheet.  The  x-axis  is  defined  as  normal  to  the  basal 
plane  (parallel  to  the  c-axis),  and  the  normal  stress  on  the  basal  plane  to  be 
used  in  equation  (4)  is  simply  cx. 

Nb  -  a,  (20) 

The  slip  plane  is  the  y-z  plane.  The  ice  will  slip  in  the  direction  of  maximum 
shear  stress  in  the  slip  plane,  defined  by  the  two  shear  stresses  and  a ^ 
Let  9b  be  the  angle  of  the  slip  direction  from  the  y-axis  towards  the  z-axis  in 
the  slip  plane. 


The  largest  shear  stresses  on  the  basal  plane  are  given  by 


Sb  ~ 


2 


+  o 


2 


(22) 


Following  equation  (4),  the  Coulomb  plastic  flow  potential  function  for  basal 
plane  slip  is  expressed  as  a  linear  function  of  the  normal  stress  Nb  and  the 
shear  stress  Sb  on  the  slip  plane. 

Fb-Sb*Ni  tanfit,  -  q,  (23) 

where  J3b  is  the  friction  angle  and  Cj,  is  the  cohesion  for  the  basal  plane  shear 
mechanism.  The  friction  angle  and  cohesion  vary  with  temperature,  grain  size, 
brine  volume,  and  the  shear  strain  rate  on  the  basal  slip  plane.  On  the  Tresca 
yield  surface,  the  plastic  flow  potential  is  independent  of  the  normal  stress,  /3b 
is  zero,  and  the  constant  is  replaced  by  oSb,  the  Tresca  shear  strength, 
which  also  varies  with  temperature,  grain  size,  brine  volume,  and  shear  strain 
rate  on  the  basal  slip  plane. 


In  terms  of  the  x-y-z  stresses,  the  Coulomb  and  Tresca  yield  criteria 
and  plastic  flow  potential  functions  are 

Coulomb :  Fb  -  +  o\  +  ox  tan0fr  -  Ch  (24) 

Tresca :  Fh  -  +  o\  -  o  sb 

Using  equations  (3)  and  (21),  the  corresponding  shear  and  normal  strain  rates 
for  the  basal  slip  mechanism  are  therefore 


Figure  4.  Timco  and  Fredcrking  ( 1 986)  Data  for  Types  A,  B.  and  C  Tests 
in  Shaded  Rectangles  with  Pertinent  Yield  Surfaces 


(31) 
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.  X 
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*  2 


±—  COS0, 
2  b 


±-~  sin0fc 


(25) 


where  X  is  a  proportionality  constant  equal  to  the  maximum  shear  strain  rate 
on  the  basal  plane.  The  factor  Vi  is  needed  in  the  derivative  of  the  plastic 
flow  potential  with  respect  to  shear  stresses  since  aV]  and  ^  must  be  treated 
as  distinct  in  the  differentiation  (Hill,  1950).  The  other  strain  rates  are  zero 
for  this  mechanism. 

The  slip  plane  for  out-of-plane  grain  slip  contains  the  z-axis,  and  the 
grains  slide  vertically.  The  shear  stress  on  the  slip  plane  causing  plastic  flow 
is  the  vector  sum  of  the  components  of  and  normal  to  the  slip  plane. 

S.  -  oa  cos0,  +  sinQ,  (26) 

where  0Q  is  the  angle  of  the  normal  to  the  slip  plane  from  the  x-axis  towards 
the  y-axis.  The  normal  stress  on  the  slip  plane  is  found  using  a  textbook 
stress  transformation  of  ax,  a ^  and  V 

N,  -  ox  cos20,  *  oy  sin20,  +  a v  sin  20,  (27) 

The  out-of-plane  plastic  flow  potential  function  in  terms  of  the  stresses  and 
0O  is  therefore 

F  -  la^cosO,  +a^stn0j  +  [ozcos20,  +o,sin2e,  20,]  tan|S,  -  C 

(28) 


The  angle  of  rotation  of  the  slip  plane  about  the  z-axis  is  determined  from  the 
maximum  value  of  the  plastic  flow  potential  function,  which  is  equivalent  to 
solving  the  following  fourth-order  equation  for  sin  0O, 

tanP  [  sin  20,  -  2  a ^  cos  20,  ]  @9) 

±  [  Oyz  cos 0,  -  0a  sine,  ]  -  0 


The  strain  rates  have  a  shear  component  parallel  to  the  shear  stress 
and  a  normal  component  perpendicular  to  the  slip  plane.  In  the  x-y-z 
coordinate  frame,  the  strain  rates  are 


-  X  cos20,  tanP , 

-  X  sin20,  tanP, 
%  -  |  sin  29,  tanp, 

.  *  i  cos0. 

*y.  '  *  |  S“9. 
e,-0 


(30) 


Example  of  Basal  Plane  Slip 

Richter-Menge  (1991)  has  measured  the  effects  of  confinement  on  the 
compressive  strength  of  aligned,  columnar  first-year  sea  ice  using  a 
conventional  triaxial  test  configuration.  All  of  the  ice  samples  were  taken 
from  the  horizontal  plane  of  the  ice  sheet,  perpendicular  to  the  long  direction 
of  the  columnar  grains.  The  angle  between  the  applied  load  and  the  mean  c- 
axis  direction  were  estimated  before  the  tests  and  measured  after  the  tests, 
resulting  in  some  scatter  from  the  desired  angles  of  0°,  45° ,  and  90°. 

For  our  analysis  of  the  Richter-Menge  data,  we  have  assumed  that 
the  stresses  normal  to  the  applied  load  are  homogeneous  and  equal  to  the 
confining  pressure.  This  is  the  conventional  approach  although 
Haythornth waite  (1960)  points  out  that  such  solutions  are  not  necessarily 
unique.  In  our  coordinates,  a  compressive  load  a  is  applied  at  an  angle  a  to 
the  x  direction  in  the  x-y  plane  with  a  confining  pressure  P.  The  stresses  on 
the  basal  plane  are 


<j  t  -  -cr  cosja  -  P  sin2 a 
-  -(a  -  P)  sin  a  cos  a 
aa  -  0 


Combining  equations  (24)  and  (31)  and  solving  for  the  maximum  load  as  a 
function  of  angle  and  confining  pressure  yields 


Coulomb : 

Tresca: 


Cb  +  P  (since  coso  +  sin2 a  tanp  j) 

since  coso  -  cos2 a  tanp  b 
oa  ♦  P  sina  cor.ce 
since  coso 


(32) 


The  material  constants  in  equation  (32)  have  been  written  as 
functions  of  the  shear  strain  rate  on  the  basal  plane;  however,  the  normal 
strain  rate  at  the  angle  cr  was  a  controlled  test  parameter,  not  the  shear  strain 
rate.  Therefore,  the  shear  strain  rate  varies  with  the  load  angle  throughout 
this  test  series.  Given  that  the  x-z  stresses  and  shear  strain  rates  are  zero,  we 
can  write  the  controlled  strain  rate  in  terras  of  the  x-y  shear  strain  rate. 

6-2  sina  coso  *  cos2 a  tanp  b )  $3) 


Richter-Menge  reported  an  unconfined  compressive  strength  of  4.27 
MPa  at  a  platen  strain  rate  of  10*3  1/s  with  the  load  applied  at  an  angle  of 
43.3°  to  the  c-axis,  and  an  unconfined  compressive  strength  of  0.59  MPa  at  a 
platen  strain  rate  of  10"^  1/s  with  the  load  applied  at  an  angle  of  393  .  These 
findings  were  summarized  with 

o  -  83.7  404®  <34) 


Assuming  that  and  vary  with  shear  strain  rate  as  in  equation  (34),  then 
equations  (32),  (33),  and  (34)  can  be  combined  to  yield 


sin  2a  -  2  cos2a  tanP  b 

0  435 


-  43 


6 

sin  2a] 


(35) 


for  no  confining  pressure.  For  the  Coulomb  yield  surface,  a  value  for  f3b  of 
30°  was  used,  and,  for  Tresca,  the  friction  angle  is  zero.  Strain  rates  normal 
to  the  load  angle  need  to  be  measured  to  determine  the  variation  of  ^  with 
shear  strain  rate. 

Equations  (32)  are  plotted  in  Figure  6  using  the  Coulomb  and  Tresca 
coefficients  from  equation  (35),  along  with  Richter-Menge^  data  for 
unconfined  compressive  strength  as  a  function  of  the  angle  between  the  load 
and  the  c-axis  at  a  strain  rate  of  10’3  1/s.  The  in-plane  grain  slip  unconfined 
compressive  strength  from  her  0°  data  is  also  plotted  as  a  straight  line.^  It 
appears  that,  in  an  unconfined  compressive  strength  test,  any  of  three  yield 
surfaces  may  be  intercepted.  It  also  appears  that  basal  slip  may  have  some 
effect  on  lower  compressive  strength  value  measured  at  90°  as  compared  to 
0°,  considered  the  scatter  of  about  5°  in  the  measured  angle.  A  test  at  an 
angle  of  60°  or  70°  would  be  very  interesting. 

Figure  7  shows  additional  data  from  Richter-Menge  with  significant 
confining  pressures.  Here  the  difference  between  the  compressive  stress  and 
the  confining  pressure  is  plotted  versus  the  angle  between  the  applied  load 
and  the  c-axis.  Two  lines  are  used  to  describe  the  in-plane  slip  yield  surface: 
one  for  an  unconfined  Coulomb  mechanism  and  one  for  a  Tresca  mechanism. 
The  confining  pressure  moves  the  basal  plane  coulomb  yield  surface  seen  in 
Figure  6  up  and  off  the  graph.  It  appears  that  basal  slip  mechanism^may 
again  contribute  to  the  lower  values  measured  at  90°  as  compared  to  0  . 
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Figure  6.  Unconfined  Compressive  Strength  due  to  Basal  Slip  and  In-Plane  Slip 
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Abstract  -  Material  constitutive  laws  are  commonly  developed 
from  stress/strain  data  which  are  obtained  from  laboratory 
measurements;  e.g.,  a  sample  of  the  material  is  stressed  in  a 
load  machine,  and  the  resulting  deformation  is  measured. 
Similarly,  developing  or  verifying  a  constitutive  law  for  polar 
pack  ice  is  best  done  by  obtaining  contemporaneous  ice  stress 
and  strain  measurements-  A  pilot  program  to  measure  these 
parameters  was  performed  during  April  1992  in  the  southern 
Beaufort  Sea.  Directional  stress  was  measured  at  two 
locations  by  gauges  imbedded  in  the  sea  ice.  Strain  was 
obtained  from  sequential  ERS-1  satellite  Synthetic  Aperture 
Radar  (SAR)  images  (typically  three  days  apart).  Unfortu¬ 
nately,  no  significant  ice  stress  or  ice  motion  events  occurred 
during  the  brief  measurement  period.  However,  a  failure 
(yield)  envelope  for  thick  ice,  based  on  ice  stress  measure¬ 
ments  taken  in  the  eastern  Arctic  Ocean,  was  developed. 
These  pilot  studies  verified  the  techniques  to  be  used  during 
an  eight-month  measurement  program  which  will  take  place 
in  the  Beaufort  Sea  during  the  1993-1994  Ice  season. 

L  INTRODUCTION 

Arctic  sea  ice  may  prove  to  be  an  important  early 
indicator  of  global  climate  change.  This  will  require  an 
understanding  (in  the  form  of  a  model)  of  Arctic  pack  ice 
dynamics  and  thermodynamics  to  understand  the  causes  of 
inter-annual  variability.  At  the  heart  of  current  ice  dynam¬ 
ics/thermodynamics  models  is  a  pack  ice  constitutive  law 
that  assumes  pack  ice  is  isotropic.  This  has  simplified 
model  calculations  but  is  a  gross  simplification  of  pack  ice 
behavior,  since  the  ice  thickness  distribution,  which 
controls  the  heat  and  mass  transfer  from  the  ocean  to  the 
atmosphere,  is  highly  anisotropic  and  irregular.  Northwest 
Research  Associates  (NWRA),  with  support  from  NASA,  is 
developing  an  anisotropic  constitutive  law  with  an  oriented 
ice  thickness  distribution  for  pack  ice  [1].  It  is  the  first 
model  of  pack  ice  to  account  for  leads  and  ridges  explicitly. 
Therefore,  the  model  can  be  validated  by  measuring  what 
happens  during  individual  events,  not  their  ensemble  as  in 
the  case  of  an  isotropic  model.  Stress  and  deformation 
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contract  No.  NASW-470I  and  by  the  Office  of  Naval  Research  under 
contract  No.  N00014-92-C-0027.  The  Polar  Science  Center  was 
sponsored  by  NASA  under  grant  No.  NAGW-2513 


measurements  for  this  modeling  approach  can  now  be 
made,  thanks  to  routine  imaging  provided  by  the  European 
Remote  Sensing  Satellite-1  (ERS-1)  SAR,  low-powered 
Global  Positioning  System  (GPS)  navigation  modules,  and 
the  development  of  ice  stress  buoys.  The  model  will  be 
verified  with  contemporaneous  field  measurements  of  sea 
ice  stress  and  ice  motion  obtained  from  sequential  SAR 
imagery  and  buoy  motions.  The  first  long-term  field 
program  to  obtain  these  data  will  be  the  Sea  Ice  Mechanics 
Initiative  (SIMI)  from  Fall  1993  to  Summer  1994, 
sponsored  by  the  Office  of  Naval  Research  (ONR). 

This  paper  gives  a  brief  description  of  the  new  pack  ice 
model  and  the  data  gathering  and  analyses  for  validation. 
Stress  data  from  the  ONR-sponsored  Coordinated  Eastern 
Arctic  Experiment  (CEAREX)  are  used  to  provide  a  first 
measure  of  the  strength  of  the  thick  ice  element  of  the  new 
model.  An  example  of  SAR  data  for  a  ridging  event  is 
presented.  In  addition,  the  data  collection  plan  for  the 
1993-1994  ice  season  is  outlined. 

IL  BACKGROUND 

The  goal  of  [1]  was  to  develop  a  constitutive  law  for  pack 
ice  which  represents  the  discontinuous  behavior  of  lead 
opening,  shearing,  and  closing.  The  continuum  model  must 
account  for  the  discontinuities  associated  with  leads  and 
have  the  same  rate  of  energy  dissipation  as  occurs  in  the 
discontinuous  sea  ice  deformations.  Therefore,  the  model 
will  have  an  anisotropic  behavior. 

The  behavior  of  the  pack  ice  element  shown  in  Fig.  1  is 
controlled  by  the  behavior  of  the  thick  ice  and  the  refrozen 
lead  ice.  The  plastic  failures  of  the  thick  ice  and  the  active 
leads  are  characterized  by  ice  strength  parameters  written  as 
stress  resultants  (the  integral  of  the  stresses  through  the  ice 
thickness)  which  are  independent  of  the  local  variations  in 
ice  thickness.  The  orientation  of  the  re  frozen  leads  and  the 
thickness  of  the  ice  therein  dramatically  affect  the  failure 
strength  of  the  pack  ice. 

The  thick  ice  is  assumed  to  be  a  perfectly-plastic  material 
with  a  Coulomb  yield  criterion  and  a  normal  flow  rule.  The 


Coulomb  yield  criterion  for  pack  ice  was  previously  used  in 
[2].  The  normal  flow  rule  has  been  used  in  many  sea  ice 
models  (e.g.,  [2].  [3],  [4],  and  [5]). 
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Fig.  1.  Element  of  pack  ice  with  a  refrozen  lead  and  a  new  lead  forming  in  the 
thick  ice. 

If  it  were  possible,  the  material  properties  of  the  pack  ice 
would  be  determined  by  laboratory  tests  in  which  the  loads 
and  deformations  imposed  on  the  test  sample  are  controlled. 
For  pack  ice,  however,  the  test  data  can  only  be  obtained  by 
interpreting  field  data.  Fig.  2a  provides  a  view  of  physical 
space  with  a  refrozen  lead  lying  along  the  x-axis  which  may 
close  (rafting  or  ridging),  shear,  open  (forming  open  water), 
or  appropriate  combinations  thereof.  Fig.  2b  shows  a. 
rectangular  plastic  failure  surface  for  this  refrozen  lead 
(expressed  in  terms  of  stress  resultants,  the  product  of  the 
stress  and  the  thickness  of  the  ice)  which  increases  in  size 
as  the  lead  ice  grows  in  thickness.  The  normal  stress 
resultant  parallel  to  the  lead  direction,  crx>  comes  out  of  the 
page  in  Fig.  2b,  forming  a  rectangular  prism.  The  yield 
surface  shown  in  Fig.  2b  is  thus  a  surface  in  3-D  stress 
space  with  the  three  stress  coordinates  being  <rx,  aty,  and 
Cjyj  the  failure  of  the  lead  is,  however,  independent  of  crx 
because  that  stress  is  taken  by  the  thick  ice.  The  normal 
flow  rule  from  plasticity  theory  then  defines  the  strain  rate 
direction  as  being  an  outward  normal  to  the  surface  of  the 
rectangular  prism,  as  illustrated. 

The  data  required  to  determine  a  point  on  the  failure 
surface  are  the  three  stresses  (cx,  Cy,  and  cr^y)  and  the 
three  corresponding  strain  rates.  To  find  the  strain  rates, 
the  ice  velocity  on  both  sides  of  the  lead  are  needed.  The 
velocity  data  could  come  from  sequential  satellite  SAR 
images  or  from  buoys  equipped  with  GPS  receivers. 
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(b)  Yield  surface  for  a  lead. 

Fig.  2.  The  stress  and  deformation  data  for  an  ice  failure  event  will  provide  a 
stress  data  point  on  the  failure  surface  and  the  strain  rate  vector  for  the  lead  (x- 
axis  is  along  the  lead) 

An  anisotropic  failure  envelope  has  been  introduced  [1], 
wherein  a  conical  failure  envelope  (in  3-D  stress  space)  is 
used  to  represent  the  thick  ice;  the  rectangular  prism  failure 
envelopes  for  the  lead  ice  are  at  45°  to  the  axis  of  the  cone. 
With  the  3-D  stress  space  axes  defined  by 


the  cone  lies  along  the  P  axis  and  the  cone  corresponds 
exactly  to  a  straight  line  in  stress  invariant  coordinates, 
since  the  pressure  invariant  is  P  and  shear  invariant  is  the 
radial  distance  out  from  the  pressure  axis  in  3-D  stress 
space,  i.e„ 

Shear  Invariant  =  VD2  +Sj  (2) 

Fig.  3a  shows  a  view  in  physical  space  of  a  lead  about  to 
be  created  in  the  x -direction.  Here  again  the  three  stresses 
and  the  corresponding  ice  velocities  are  needed  to  get  a 
point  on  the  yield  surface  and  a  strain  rate  vector.  The 
thick  ice  failure  surface,  a  cone  in  stress  three-space,  is 
shown  in  Fig.  3b,  along  with  a  2-D  cut.  The  normal  flow 
rule,  as  illustrated,  anticipates  a  combined  opening  and 
shearing  motion  when  the  thick  ice  fails,  thereby  creating  a 
new  lead. 


m.  DATA  GATHERING  AND  ANALYSIS 

To  calculate  the  failure  surfaces  for  various  ice 
thicknesses,  one  needs  to  have  measurements  of  directional 
stress,  concurrent  ice  strain,  and  the  thickness  of  the 
deformed  ice  during  periods  of  deformation.  A  pilot 
program  was  conducted  in  the  Beaufort  Sea  during  April 
1992  to  gather  such  data  as  part  of  the  ONR-sponsored 
LEADEX  Program.  Stress  was  measured  by  NWRA 
investigators  at  two  locations.  Ice  motion  was  available 
from  sequential  SAR  imagery  and  from  ARGOS  satellite 
positioning  system  buoys.  Unfortunately,  during  the 
measurement  period  no  geophysical  stress  events  were  seen 
in  the  data  and  no  significant  ice  failure  events  were 
observed.  The  gathering  of  concurrent  ice  stress  and  strain 
during  deformation  events  would  have  to  wait  until  the  field 
experiments  planned  for  the  1993-1994  ice  season. 

To  prepare  for  the  future  field  experiments,  techniques 
have  been  developed  to  analyze  stress  data  and  SAR 
imagery  individually.  The  procedures  and  results  of  these 
analyses  are  described  below. 

A.  Stress  Data  Analysis 

The  measurement  of  the  large  scale  failure  stress  of  thick 
pack  ice  is  difficult,  but  a  method  to  infer  it  is  presented 
here.  Pack  ice  in  the  Arctic  Ocean  is  generally  composed  of 
thick  ice  and  refrozen  lead  ice.  The  thick  ice  consists  of 
multi-year  (MY)  floes,  thick  first  year  (FY)  ice,  ridges,  and 
rubbled  areas.  This  ice  has  an  average  thickness  of  about  3 
meters,  varying  from  1  to  50  meters  or  so.  The  refrozen 
lead  ice  is  generally  thin  (less  than  2  meters  thick)  and 
saline. 


(a)  Physical  space. 
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(b)  Conical  yield  surface  for  thick  ice. 

Fig.  3.  A  lead  created  along  the  x-axis  results  from  the  stress  stale  being  on  the 
conical  failure  surface. 

The  stress  required  to  fail  the  refrozen  lead  ice  is  smaller 
than  the  stress  required  to  fail  the  thick  ice,  although  the 
thick  ice  has  cracks  and  other  flaws  so  that  it  is  not  as 
strong  as  a  homogeneous,  continuous  block  of  ice.  For 
certain  ice  stress  states,  however,  the  thick  ice  fails.  Thus, 
there  are  three  possible  ice  states,  listed  in  order  of 
increasing  stress: 


1.  No  ice  is  failing; 

2.  Refrozen  lead  ice  is  failing; 

3.  Thick  ice  is  failing. 

Since  pack  ice  cannot  exist  with  a  larger  stress  than  that 
required  to  fail  the  thick  ice,  the  plastic  failure  envelope  of 
the  thick  ice  should  provide  an  upper  bound  fa-  the  ice 
stress  data.  The  Coordinated  Eastern  Arctic  Experiment 
(CEAREX)  stress  data  ([6]  and  [7])  offers  an  opportunity  to 
identify  features  of  the  pack-ice  failure  envelope.  Fig.  4 
shows  the  hourly  averaged  stress  resultant  invariants  (from 
[6])  measured  in  the  middle  of  a  large  multi-year  floe  about 
1.6  m  thick  for  the  months  of  October  and  November,  1988, 
in  the  eastern  Arctic  Ocean.  Stress  resultant,  the  integral  of 
the  horizontal  stress  through  the  thickness,  is  taken  to  be 
the  measured  ice  stress  times  the  thickness.  With  the 
pressure  invariant  on  the  horizontal  axis  and  the  shear 
invariant  on  the  vertical  axis,  the  line  at  45°  is  the  boundary 
between  compression  in  all  directions  (below  the  line)  and 
tension  in  some  direction  (above).  Since  the  line  at  45°  also 
represents  a  uni-axial  compression  stress  state,  these  pack 
ice  data  indicate  that  pack  ice  is  near  a  state  of  uni-axial 
compression  most  of  the  time. 


Fig.  4.  Hourly  averaged  stress  resultant  invariants. 


For  the  anisotropic  pack  ice  constitutive  law,  the  x-axis 
of  Fig.  4  corresponds  to  the  axis  of  the  cone  and  the  y-axis 
corresponds  to  the  radial  distance  (in  stress  resultant  space) 
out  from  the  axis  of  the  cone  (refer  to  Fig.  3b).  The  cone  in 
Fig.  3b  is  thus  represented  by  the  conical  failure  surface  line 
in  Fig.  4  which  intersects  the  x-axis  at  about  50  kN/m  with 
a  slope  of  about  -0.65.  This  failure  envelope  is  insensitive 
to  the  data  averaging  period  over  the  range  from  2  minutes 
to  12  hours.  In  this  case,  the  points  below  this  line  (within 
the  cone  failure  envelope)  can  be  interpreted  as  a  time  of  no 
ice  Mure  or  a  time  when  a  refrozen  lead  has  offered  a 
failure  envelope  interior  to  the  cone.  This  failure  envelope 
corresponds  to  a  uni-axial  (unconfined)  compressive 


strength  of  about  180  kN/m.  In  contrast,  laboratory  sea  ice 
samples  have  a  uni-axial  compressive  strength  of  between  1 
and  6  Mpa,  which  corresponds  to  1600  to  9600  kN/m  in 
1.6-m  pack  ice.  Apparently  the  cracks  and  other 
inhomogeneities  in  the  thick  ice  greatly  reduce  its  strength. 

The  assertion  that  the  CEAREX  ice  stress  measurements 
represent  geophysical  stress  was  examined.  In  the 
installation  of  the  stress  gauges,  a  procedure  was  developed 
and  followed  to  achieve  a  measurement  of  geophysical  ice 
stress  which  is  consistent  with  a  procedure  outlined  in  [8]. 
They  suggest  placing  the  stress  sensors  in  the  middle  of  a 
large  floe,  in  a  relatively  flat  area  free  from  major  cracks,  at 
a  depth  from  one-third  to  one-half  the  ice  thickness  from 
the  surface,  such  that  the  whole  floe  acts  as  a  pressure 
transducer.  Stress  sensors  near  the  middle  of  a  large  floe 
avoid  the  higher  stresses  found  along  the  edges  of  the  floe 
[9].  A  relatively  flat  area  is  needed  to  avoid  stress  vari¬ 
ations  associated  with  thickness  variations.  Placing  the 
sensor  at  mid-depth  to  a  great  extent  avoids  tensile  and 
compressive  stresses  resulting  from  bending  stresses  caused 
by  temperature  changes  near  the  upper  surface  of  the  ice 
(the  so-called  thermal  stresses)  [10], 

There  are  other  potential  error  sources  in  the  data: 
inclusion  factor,  transducer  accuracy,  zero  drift,  zero 
accuracy,  tensile  stress  measurements,  and  direct  thermal 
effects.  The  inclusion  factor  is  the  ratio  of  pressure  in  the 
stress  sensor  to  the  far- fie  Id  ice  stress.  Based  on  the  stress 
gauge  calibration  measurements,  the  inclusion  factor  for 
these  gauges  exceeds  0.95,  which  includes  any  errors  in 
transducer  accuracy.  The  zero  drift  was  negligible,  but  the 
accuracy  of  the  zero  stress  measurement  is  not  as  well 
known.  For  the  data  shown  in  Fig.  4,  a  zero  stress  state  was 
assumed  to  occur  at  the  end  of  a  two-week  period  with  no 
winds,  no  ice  motion,  and  all  the  stress  gauges  reading  near 
their  minimum  compressive  value.  The  stress  data  look 
somewhat  different  if  the  zero  stress  state  is  assumed  to 
occur  at  the  beginning  of  the  same  time  period  when  the 
floe  was  surrounded  by  open  Wats’  since  one  of  the  stress 
gauges  reads  about  20  kPa  more  compression  at  the 
beginning  than  at  the  end  of  the  time  period.  In  addition,  it 
is  unclear  from  the  data  how  well  these  gauges  can  measure 
tensile  stress.  These  data  show  the  ice  experienced  sane 
tension  but  never  far  into  positive  pressure.  Direct  thermal 
effects  on  the  components  of  the  stress  sensor  above  the  ice 
surface  are  a  concern,  but  there  was  no  sunlight  during  the 
test  period  and  the  above-ice  components  wee  buried  in 
snow.  In  all,  we  estimate  that  the  sum  of  the  errors  is  about 
30  kN/m.  From  the  discussion  above,  it  seems  clear  this 
methodology  can  be  used  to  measure  the  large  scale  Mure 
stress  of  pack  ice.  In  this  way,  a  small  scale  measurement 
(ice  stress)  is  used  to  determine  a  large  scale  parameter  (the 
strength  of  thick  pack  ice). 

The  Mure  envelope  may  vary  somewhat  from  season  to 
season  and  from  region  to  region.  More  will  be  known 


when  sufficient  stress  and  deformation  data  are  available  to 
relate  the  failure  stress  to  the  ice  motion  at  failure. 

B.  Ice  Motion  Measurements  from  SAR  Imagery. 

The  Geophysical  Processor  System  at  the  Alaska  SAR 
Facility  (ASF)  automatically  generates  fields  of  sea  ice 
motion  using  pairs  of  HRS- 1  SAR  images  as  input  These 
low  resolution  image  pairs,  typically  separated  in  time  by 
three  days  (43  satellite  revolutions),  have  been 
radiometrically  corrected  and  geolocated  onto  a  polar 
stereographic  projection.  The  image  pixel  size  is  100 
meters  and  the  area  coverage  is  approximately  100  square 
kilometers.  The  ice  tracking  algorithm  automatically 
selects  the  image  pairs  form  the  ASF  archive  using  a  model 
of  ice  drift  based  on  the  geostrophic  wind  field.  A  regular 
array  of  grid  points  with  five  kilometer  spacing  is  also 
defined  on  the  same  polar  stereographic  projection.  The 
algorithm  proceeds  by  attempting  to  track  the  ice  at  each 
grid  point  to  its  new  location  in  the  later  image.  In  the 
central  pack  ice,  a  cross-correlation  technique  is  used  to 
identify  corresponding  points.  Near  the  coast  or  in  the 
marginal  ice  zone,  a  feature  matching  technique  is  used  as 
well  as  cross-correlation,  since  ice  floes  can  undergo  large 
rotations  in  those  areas.  Finally,  a  set  of  consistency  checks 
and  filters  is  applied  to  remove  bad  matches.  The  result  is  a 
set  of  corresponding  points,  or  equivalently,  displacement 
vectors,  on  a  regular  five  kilometer  grid  [11].  Each  ice 
motion  product  contains,  on  average,  over  200  matches  or 
vectors.  More  than  4000  ice  motion  products  have  been 
archived  at  the  ASF  since  the  Geophysical  Processor  System 
became  operational  in  February  1992.  These  products  cover 
the  period  from  September  1991  through  September  1992  in 
the  Beaufort  and  Chukchi  Seas  and  the  central  Arctic 
Ocean.  More  ice  motion  products  are  generated  every 
week. 

During  the  Commissioning  Phase  and  the  first  Ice  Phase 
of  ERS-1,  from  July  1991  through  March  1992,  the  satellite 
was  in  a  three-day  repeat  orbit,  providing  many  thousands 
of  image  pairs  for  ice  motion  analysis.  During  the  Multi¬ 
disciplinary  Phase,  from  April  1992  through  December 
1993,  the  satellite  is  in  a  35-day  repeat  orbit  with  a  drifting 
three-day  subcycle.  This  still  provides  good  swath  overlap 
at  far  northern  latitudes,  but  the  potential  number  of  ice 
motion  products  is  less  than  that  provided  during  the  exact 
three-day  repeat  orbit.  Finally,  ERS-1  will  be  in  the  second 
Ice  Phase  from  January  through  March  1994,  again 
providing  large  amounts  of  three -day  repeat  imagery. 

Reference  [12]  surveys  the  uses  of  ice  motion  data  from 
SAR  and  emphasizes  their  importance.  SAR-derived  ice 
motion  is  being  combined  with  data  from  drifting  buoys  to 
improve  our  understanding  of  the  circulation  of  the  Arctic 
Ocean.  Lead  dynamics  are  being  studied  in  connection  with 
ice  deformation.  Ice  motion  also  plays  an  important  role  in 


the  equation  governing  the  ice  thickness  distribution,  which 
describes  the  balance  between  ice  production,  advection, 
and  melting. 

The  data  from  an  ice  motion  product  could  be  displayed 
by  simply  drawing  all  the  displacement  vectors  on  a  map. 
But  sea  ice  tends  to  move  in  large  rigid  pieces  over  long 
distances,  so  it  is  difficult  to  see  the  small  differential 
motion  between  adjacent  vectors.  By  subtracting  the  mean 
displacement  vector  and  displaying  the  resulting  deviations 
from  the  mean,  the  differential  motion  stands  out  more 
clearly.  Fig.  5  shows  a  field  of  such  deviations,  magnified 
by  a  factor  of  five  relative  to  the  grid  spacing  of  five 
kilometers.  This  scene  is  located  just  off  the  southwest  tip 
of  Prince  Patrick  Island  in  the  Canadian  Archipelago,  and 
represents  the  ice  motion  between  March  1 1  and  March  14, 
1992.  The  mean  displacement  vector  is  shown  in  the  upper 
left  comer,  approximately  13.8  kilometers  in  length.  A 
linear  zone  of  convergence,  where  the  arrows  point  toward 
one  another,  runs  from  the  upper  left  to  the  lower  right. 

The  ice  deformation  in  this  and  2000  other  ice  motion 
products  have  been  analyzed  [13].  They  computed  area- 
averaged  values  of  the  strain  invariants  e /  (divergence)  and 
en  (shear)  as  well  as  I  e  I  =  (e{2  +  £//)1/2  and  Q  = 
arctan (euJefr  The  magnitude  of  deformation  is  measured 
by  lei,  and  the  type  of  deformation  is  characterized  by  q: 
pure  divergence  (q  =  0°);  uniaxial  extension  (q  =  45°);  pure 
shear  (q  =  90°);  uniaxial  compression  (q  =  135°);  and  pure 
convergence  (q  =  180°).  For  the  ice  motion  in  Figure  5, 
which  covers  an  area  of  5475  km2, 1  e  I  =  2.6%  and  q  =  129° 
for  the  three-day  period.  The  net  convergence  is  1.6%, 
representing  an  area  loss  of  89  km2  which  is  presumably 
accounted  for  by  the  closing  of  leads  and  ridge  building. 

On  a  smaller  scale,  strain  invariants  can  be  computed 
from  displacement  (or  velocity)  derivatives  that  in  turn  are 
computed  from  simple  finite  difference  expressions.  If 
and  z*2  are  the  x-components  of  displacement  (or  velocity)  at 
two  adjacent  grid  points,  then  (^  -  u^)fDx  is  approximately 
du/dx  at  the  midpoint.  It  is  not  possible  to  compute 
derivatives  at  scales  smaller  than  Dx  =  5  km  with  the 
Geophysical  Processor  System  data,  but  it  is  possible  to 
acquire  the  SAR  images  and  track  ice  features  whose 
separation  distance  is  very  small.  This  has  been  done  in 
[14]  using  a  semi-automated  tracking  program  based  on 
area  correlation.  As  a  practical  matter,  the  geolocation 
errors  in  the  images  and  the  uncertainty  in  tracking  the  ice 
features  limit  the  separation  distance  to  a  minimum  of 
several  pixels  (several  hundred  meters).  The  full  resolution 
SAR  images  offer  no  advantage  over  the  low  resolution 
images  since  the  geolocation  errors  are  on  the  order  of  100 
meters  in  both,  and  the  precise  identification  of  features  in 
the  full  resolution  images  is  complicated  by  the  presence  of 
speckle.  Thus,  the  smallest  scale  on  which  strain  invariants 
can  be  computed  is  on  the  order  of  half  a  kilometer,  using 


the  low  resolution  images. 


Fig.  5.  Example  of  differential  motion  vectors  derived  from  SAR  imagery, 
lire  differential  vectors  are  magnified  five  times  relative  to  the  5x5  km  grid. 

rv.  DATA  COLLECTION  DURING  1993-1994  ICE  SEASON 

During  September  1993,  three  ice  stress  stations  will  be 
deployed  by  NWRA  in  the  southern  Beaufort  Sea  as  part  of 
ONR's  Sea  Ice  Mechanics  Initiative  (SIMI).  These  stations 
will  provide  data  through  April  1994  and  perhaps  longer. 
During  the  same  period,  ice  motion  data  will  be  gathered  by 
fifteen  GPS  buoys  placed  on  floes  in  the  vicinity  of  the  SIMI 
camp  by  the  National  Oceanic  and  Atmospheric 
Administration's  Pacific  Marine  Environmental 
Laboratories.  Ice  strain  will  be  used  from  these  buoys  and 
from  ice  motion  vectors  from  the  ASF  which  will  be 
supplemented  by  more  detailed  motion  analyses  done  by  the 
Polar  Science  Center  for  selected  ice  motion/stress  events. 
The  combined  SAR  derived  motions  and  buoy  meas¬ 
urements  will  be  used  with  the  directional  stress  data  to 
construct  failure  surfaces  needed  by  the  anisotropic 
constitutive  model. 

V.  CONCLUSIONS 

1.  Ice  motion  and  field  stress  measurements  can  be 
combined  to  define  the  failure  surfaces  for  pack  ice  of 
various  thicknesses. 


2.  In  the  absence  of  motion  data,  the  failure  envelope  for 
thick  pack  ice  can  be  obtained  from  stress  measurements 
alone. 

3.  We  measured  the  unconfined  compressive  strength  of 
thick  ice  in  the  eastern  Arctic  Ocean  during  October  and 
November  1988  to  be  120  kN/m,  much  less  than  a  strength 
based  on  laboratory  tests  of  sea  ice  samples. 

4.  The  differential  ice  motion  data  generated  routinely  from 
sequential  ERS-1  SAR  imagery  can  be  used  to  help  define 
failure  surfaces  of  sea  ice  in  refrozen  leads. 
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ABSTRACT 

The  Sea  Ice  Mechanics  Initiative  (SIMI)  is  a  five-year 
U.  S.  Navy  Office  of  Naval  Research  (ONR)  program.  The 
main  SIMI  field  experiment  was  in  the  Beaufort  Sea  from 
September  1993  through  April  1994 ,  with  numerous  other  small 
field  experiments ,  laboratory  experiments ,  and  modeling 
efforts.  The  goals  of  this  program  are  to  understand  sea  ice 
constitutive  laws  and  fracture  mechanics  over  the  full  range  of 
geophysical  scales ,  to  determine  the  scaled  responses  to 
applied  external  forces,  and  to  develop  physically-based 
constitutive  and  fracture  models .  About  twenty  principal 
investigators  are  working  to  achieve  these  goals  along  with 
their  associates.  The  SIMI  experiments  include  ice  stress, 
strain,  strength,  tilt,  motion ,  temperature,  and  response  to 
controlled  load  experiments.  This  work  will  provide  new  ice 
information  for  loads  on  structures  and  for  Arctic  operations, 
both  private  and  government. 

INTRODUCTION 

The  goals  of  the  Sea  Ice  Mechanics  Initiative  (SIMI) 
Arctic  research  program  are  stated  in  the  SIMI  Summary  Plan 
(Curtin,  1993).  These  are: 

•  Understand  sea  ice  constitutive  laws  and  fracture  me¬ 
chanics  over  the  full  range  of  geophysical  scales  and 
determine  the  scaled  responses  to  applied  external 
forces. 


Nomenclature  and  references  at  end  of  paper. 


•  Develop  physically-based  constitutive  and  fracture 

models. 

These  goals  were  developed  at  the  Sea  Ice  Mechanics  Workshop 
(Curtin,  1991),  along  with  an  assessment  of  the  limits  of  current 
understanding,  a  list  of  priority  research  issues,  and  state-of-the- 
art  methods  for  the  many  aspects  of  sea  ice  mechanics.  Individ¬ 
ual  projects  within  the  program  were  defined  around  the  scale 
size  used  in  the  ice  mechanics  work. 

This  paper  describes  SIMI,  focusing  on  the  field  pro¬ 
gram  initiated  in  September  1993. 

BACKGROUND 

The  Office  of  Naval  Research  held  a  Sea  Ice  Mechan¬ 
ics  Workshop  at  Airlie,  Virginia  on  November  12-14,  1990, 
where  sea  ice  mechanics  was  examined  from  different  perspec¬ 
tives,  both  methodological  (theoretical  modeling,  laboratory 
experiments,  and  field  observations)  and  behavioral  (material, 
structural,  and  acoustical  aspects).  The  Workshop  Proceedings 
(Curtin,  1991)  outline  limits  of  current  understanding,  priori¬ 
tizes  research  issues,  and  documents  innovative  methods  of 
approach.  The  workshop  participants  included  representatives 
from  government  agencies  and  laboratories  as  well  as  from 
universities,  the  oil  industry,  and  other  private  sector  engineers 
and  scientists.  The  Workshop  Proceedings  were  used  to  shape 
the  five-year  (1992-96)  SIMI  program. 

ONR  sponsored  another  workshop  at  Dunsmuir  Lodge, 
Sidney,  British  Columbia,  Canada  on  October  12-14,  1992 
(Curtin,  1993).  During  this  workshop,  the  detailed  plan  for 
SIMI  was  worked  out.  The  main  field  experiment  included  a 
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drifting  camp  in  the  Beaufort  Sea  collecting  data  from  Fall  1993 
through  Spring  1994.  The  SIMI  plan  and  the  field  program  are 
described  below  in  a  quote  from  Curtin  (1991). 

"A  major  question  in  understanding  and  predicting  the 
mechanical  behavior  of  Arctic  sea  ice  involves  the  effect  of  size 
on  the  dynamics.  The  notion  that  the  scales  of  a  specimen  alone 
influences  the  operative  physics  has  fundamental  implications. 
The  prime  example  is  the  pressure  versus  contact  area  relation¬ 
ship,  suggesting  a  decrease  in  strength  with  an  increase  in  size. 
There  was  consensus  that  obtaining  a  comprehensive  quantita¬ 
tive  relationship  spanning  from  microscale  to  structural  scale 
behavior  is  too  ambitious  on  a  five-year  horizon;  however, 
establishing  the  linkage  between  a  few  scales  (e.g.,  centimeters 
to  meters,  meters  to  tens  of  meters)  is  achievable  and  should  be 
pursued. 

"Theoretical  modeling  is  essential  in  relating  behavior 
across  scales  and  developing  parameterizations.  Improved 
models  are  required  for  both  continuum  and  fracture  behavior. 
Continuum  and  fracture  models  should  be  linked  where  possi¬ 
ble.  Small-scale  laboratory  experiments  are  required  to  link 
micromechanical  properties  with  physical  behavior  and  to 
identify  mechanisms  that  may  be  important  across  scale. 
Intermediate  scale  laboratory  tests  under  carefully  controlled 
conditions  can  provide  a  crucial  link  between  small  and  large 
scales.  There  was  consensus  on  the  need  to  establish  the  critical 
scales  at  which  microscale  processes  become  measurably 
important  to  larger  scale  behavior. 

"Brittle  or  ductile  sea  ice  behavior  depends  to  a  large 
extent  on  loading  rates,  which  in  nature  typically  range  from 
10'7  to  10"4  s'1.  Continuum  behavior  at  different  rates,  particu¬ 
larly  at  larger  scales,  needs  further  investigation.  Rate  effects 
on  fracture  need  to  be  quantified. 

"While  it  is  well  known  that  temperature  can  deter¬ 
mine  elastic  or  viscous  plastic  behavior  under  isothermal 
conditions,  the  behavior  of  sea  ice  with  realistic  temperature 
gradients,  as  in  the  Arctic,  remains  to  be  characterized.  The 
effects  of  temperature  on  fracture  need  to  be  resolved.  Porosity, 
salinity,  and  microcrack  distribution  are  important  at  small  and 
medium  scales.  Gradients  in  these  properties,  in  addition  to 
those  in  thickness,  degree  of  ridge  consolidation,  and  macroflaw 
distribution  at  larger  scales,  need  to  be  ordered  and  quantified 
in  the  effects  on  the  governing  mechanics. 

"Coherent  and  economical  field  tests  designed  to  de¬ 
finitively  relate  forcing  and  response  under  the  range  and  limits 
of  conditions  in  nature  are  essential.  Controlled  loading  tests  in 
the  field  are  necessary  to  understand  large-scale  behavior.  In- 
situ  ice  properties  and  stresses  can  be  monitored  to  some  extent 
with  current  technology.  Methods  for  more  complete  spatial 
mapping  of  these  fields  are  needed.  Techniques  to  measure 
three-dimensional  flaw  generation  and  distribution  within  full 
scale  floes  need  development 


"The  roles  of  mixed  mode  and  non-simultaneous  fail¬ 
ure  in  reducing  large-scale  loads  needs  to  be  understood,  as  well 
as  the  important  mechanisms  in  the  initiation  of  long  fractures 
(leads)  in  pack  ice.  The  relative  importance  of  bending,  buck¬ 
ling,  crushing,  and  crack  nucleation/propagation  in  overall 
failure  needs  to  be  determined.  Factors  causing  ice  of  the  same 
thickness  to  ridge  in  some  cases  and  raft  (including  finger 
rafting)  in  others  need  to  be  determined.  Thermal  stresses  may 
be  as  important  as  mechanical  stresses  in  preconditioning  the 
ice  for  the  propagation  of  cracks.  The  role  of  existing  flaws  and 
the  residual  stress  state  in  the  generation  and  propagation  of 
long  fractures  requires  study. 

"The  transition  from  laboratory  to  floe  scales  is  best 
addressed  within  a  hypothesis-testing  framework.  One  hy¬ 
pothesis  is  that  the  behavior  of  samples  under  tensile  stress  and 
the  distribution  of  this  stress  is  most  important  in  floe-scale 
cracking  processes.  Therefore,  these  stress  conditions  should  be 
examined  experimentally  and  parameterized  in  larger  scale 
relationships  that  can  be  tested  in  floe  experiments. 

"Bridging  the  gap  in  understanding  of  ice  processes 
from  floe  to  regional  scales  may  come  from  relating  floe  meas¬ 
urements  to  a  viable  statistical  measure  of  the  structure  of 
regional  scale  distributions  of  ice.  The  measurements  of 
environmental  forcing  and  regional  strain,  for  instance,  can  be 
related  to  average  pack  ice  forces  measured  on  the  floe.  Careful 
study  of  the  relationships,  particularly  during  deformation, 
should  lead  to  improved  regional  constitutive  laws.  Also, 
fracture  models  developed  for  a  single  floe  may  be  applicable  to 
a  collection  of  floes.  Thus,  selected  geophysical  parameters 
obtained  from  aerial  or  satellite  imagery  of  a  region  combined 
with  environmental  forcing,  will  produce  time-dependent 
theoretical  and  numerical  models  that  can  predict  with  improv¬ 
ing  degrees  of  skill  the  spatial  distribution  of  new  leads  and 
ridges." 

MODELING  SUMMARY 

The  modeling  studies  included  in  the  SIMI  program 
are  summarized  below,  grouped  by  scale  size. 

100  KM  SCALE  (PACK  ICE1 

•  Coupled  air-ice-ocean 

•  Fracture  characteristics 

•  Deformation  characteristics 

10  KM  SCALE  fFLOE  CLUSTER) 

•  Granular  media  rheology 

•  Anisotropic  constitutive  law 

1  KM  SCALE  (FLOE.  RIDGE.  LEAD1 

•  Ridge  building  energetics 

•  Fracture  mechanics,  size  effect 

1-100  M  SCALE  (MACROCRACK) 

•  Thermal  fracture 

•  Physically-based  constitutive  laws 
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•  Fracture  mechanics,  size  effects 

1  CM  -  1  M  SCALE  fMICROCRACK-CRYSTALl 

•  Ductile-brittle  transition,  failure  surfaces 

•  Crack  nucleation 

MEASUREMENT  SUMMARY 

The  experimental  studies  in  the  SIMI  program  are 
summarized  below,  grouped  by  scale  size. 

100  KM  SCALE  fPACK  TCF.1 
FORCING 
Mechanical 

•  Regional  surface  winds  (NOGAPS,  NMC) 

•  Regional  surface  currents  (GPS  buoy  array) 
Thermodynamic 

•  Surface  temperature  (Argos  buoy  array) 
DEFORMATTON/FRACTURE 

•  Point  strain  field  (Argos  buoy  array) 

•  Image  strain  field  (satellite  AVHRR,  SSMI,  SAR) 
PROPERTIES 

•  Thickness  distribution  (low  frequency  acoustics) 

10  KM  SCALE  (FLOE  CLUSTER’! 

FORCING 

Mechanical 

•  Stress  on  floes  (GPS  buoy  array) 

Thermodynamic 

•  Surface  temperature  (15  element  GPS  buoy  array) 
DEFORMATION/FRACTURE 

•  Image  strain  field  (satellite  SAR/ASF) 

•  Image  strain  field  (satellite  AVHRR,  SSMI) 

•  Point  strain  field  (15  element  GPS  buoy  array) 
PROPERTIES 

•  Thickness  distribution  (AUV) 

•  Ridge  event  distribution  (SAR/ASF,  AVHRR) 

•  Ice  type  distribution  (SAR/ASF) 

1  KM  SCALE  (FLOE,  RIDGE.  LEAD1 
FORCING 
Mechanical 

•  Wind,  ocean  stress  (main  floe  station) 

•  Ice  stress  array  (24-30  stress  sensors) 

Thermodynamic 

•  Surface  temperature,  radiation  (main  floe  station) 
DEFORMATION/FRACTURE 

•  Strain,  tilt,  acceleration  (sensor  clusters) 

•  Ridging  event  detection,  localization  (hydrophone  ar¬ 
ray) 

PROPERTIES 

•  Morphology,  ridge-lead  distribution  (main  floe  survey) 

•  Thickness  distribution  (AUV,  submarine) 

•  Ridge  event  distribution  (SAR/ASF,  AVHRR,  buoys 

1-100  M  SCALE  (MACROCRACKl 
FORCING 


Mechanical 

•  Controlled  loads  (stress  gauge  array) 

Thermodynamic 

•  Ice  temperature  profiles  (thermistor  arrays) 
DEFORMATION/FRACTURE 

•  Event  localization  and  magnitude  (hydrophone,  geo¬ 
phone  arrays) 

•  Ridging  event  seismo-acoustic  inversion  (geophone, 
hydrophone  arrays) 

•  Ridging  deformation  (high  resolution  GPS,  stress  sen¬ 
sors,  tilt-strain  meters,  and  inflatable  air  jack) 

•  Fracture  energy,  failure  modes  (first-year  ice  controlled 
load) 

•  Fracture  energy  distribution  (first-year  ice  controlled 
load-acoustic,  seismic  response:  borehole  jack,  hydro¬ 
phone,  geophone  arrays) 

PROPERTIES 

•  Seasonal  evolution  of  elastic  and  shear  wave  properties 
and  fissures  on  flat  ice  (controlled  impulse  source) 

•  Seasonal  evolution  of  elastic  and  shear  wave  properties 
and  fissures  in  ridged  ice,  crack  event  localization, 
failure  mode  detection  (acoustic  tomography) 

•  Seasonal  evolution  of  ridge  elastic,  structural  proper¬ 
ties  (coring,  thermistor  arrays) 

•  Seismo-acoustic  scattering  from  ridges  (pencil  beam 
and  explosive  sources,  hydrophone  and  geophone  ar¬ 
rays) 

•  Block  friction  (in-situ  tests) 

1  CM  -  1  M  SCALE  ('CRYSTAL-MICROCRACK') 

FORCING 

Mechanical 

•  Experiments  under  controlled  loads 

Thermodynamic 

•  Surface  temperature  (NWS) 
DEFORMATION/FRACTURE 

•  Flexure,  fracture  related  to  microstructure  (in-situ 
samples) 

•  Flexure,  fracture  related  to  microstructure  (pit-pond 
tests) 

•  Fracture  related  to  microstructure  (testing  machine) 

•  Low  strain  rate  creep  (testing  machine) 

•  Internal  friction,  anelasticity  (cyclic  loading) 

•  Multiaxial  compressive  strength  (pit-pond,  testing  ma¬ 
chine) 

•  3-D  failure  surfaces  (testing  machine) 

PROPERTIES 

•  Microstructure  linked  with  growth  conditions  (in-situ 
samples) 

Figure  1  shows  the  SIMI  experiment  schedule,  includ¬ 
ing  tests  near  Resolute,  Northwest  Territories,  Canada,  near 
Barrow,  Alaska,  and  in  the  Beaufort  Sea.  A  workshop  to 
present  and  discuss  preliminary  SIMI  results  will  be  scheduled 
in  the  first  half  of  1995. 
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Figure  1.  SIMI  Experiment  Schedule 


BEAUFORT  SEA  FIELD  PROGRAM 

The  field  program  in  the  Beaufort  Sea  was  deployed 
from  an  icebreaker,  the  U.S.  Coast  Guard  Cutter  Polar  Star.  A 
temporary  camp  was  established  on  the  pack  ice  north  of 
Alaskaduring  the  latter  half  of  September  1993.  The 
approximately  25-30  tons  of  camp  cargo  moved  to  the  site  by 
Polar  Star  allowed  the  construction  and  subsequent  support  of  a 
20-person  temporary  ice  camp.  The  facilities  included  heated 
insulated  quarters,  electric  power,  communications,  navigation, 
surface  transportation,  and  food,  fuel  and  lumber  supplies  to  last 
beyond  the  scheduled  December  1  closure  of  camp.  While  the 
camp  was  being  set  up,  investigators  riding  the  ship  started 
setting  up  their  experiments  on  the  ice.  The  Polar  Star 
helicopters  supported  operations  to  select  the  floe  for  the  main 
camp,  deploy  the  position/stress  buoy  array,  perform  rafting 
studies,  and  obtain  a  photo  mosaic  of  the  buoy  area. 

Even  before  boarding  the  Polar  Star,  it  was  recognized 
that  the  camp  would  have  to  further  north  and  west  that  origi¬ 
nally  planned  due  to  an  unusually  light  ice  year  in  the  Beaufort 
Figure  2  illustrates  the  degree  to  which  1993  was  an  unusual  ice 
year  (adapted  from  Brower  et  al.,  1977).  The  jagged  solid  line 
is  the  ice  edge  taken  from  an  AVHRR  image  for  September  28, 
1993;  the  circle  is  the  initial  location  of  the  SIMI  ice  camp. 
The  ice  edge  at  140°W  was  100  km  further  north  than  the  limit 
shown  on  the  chart,  making  it  impossible  to  deploy  the  ice  camp 
at  the  desired  location  of  73°N,  138°W.  For  an  average  year, 
the  desired  camp  location  would  have  been  over  200  km  north 
of  the  ice  edge,  presumably  in  heavy  multi-year  ice.  The  new 
location  to  look  for  a  Main  Camp  floe  was  set  as  75°N,  141°W. 

The  science  crew  boarded  the  Polar  Star  by  helicopter, 
and  the  cargo  by  barge,  on  September  16,  1993,  from  Dead- 
horse,  Alaska.  Some  scientific  work  was  performed  on  the 


Polar  Star  in  the  wet  lab  and  science  room.  The  ice  edge  was 
crossed  near  74°00'N,  143°30’W.  An  ice  floe  was  drilled  at  74° 
50’N,  141°10rW;  it  was  found  to  be  rotten  ice  with  a  thin  layer 
on  melt  ponds.  SLAR  airplane  and  helicopter  overflights  (with 
landings  to  drill  the  ice)  were  used  to  search  for  an  adequate 
floe  in  the  vicinity  of  75°N,  141°W.  It  was  recognized  that, 
without  sufficient  ice  for  a  runway,  the  science  team  could  not 
be  left  at  the  ice  camp  as  planned.  After  another  day  of  helicop¬ 
ter  flights,  a  floe  was  selected  at  75°10'N,  142°10'W  with  a 
promising  future  runway  with  ice  30-90  cm  thick.  The  camp 
construction  was  begun  on  September  20,  in  balmy  -5°C 
weather.  Interested  members  of  the  science  team  were  provided 
site  selection  helicopter  overflights.  The  Polar  Star  helicopters 
were  also  used  to  deploy  twelve  GPS  and  four  stress  buoys  at 
distances  up  to  10  km  around  the  main  floe.  Alternative  run¬ 
ways  were  identified  and  marked  on  nearby  floes  and  radio 
communications  were  established  between  the  ice  camp,  ship, 
and  shore.  The  10-km  square  buoy  array  area  was  mapped  and 
photographed  from  the  helicopters.  A  few  science  experiments 
were  started  in  the  vicinity  of  the  camp.  Based  on  ice  condi¬ 
tions,  the  decision  was  made  to  close  the  camp  when  the  Polar 
Star  departed  and  reopen  it  on  October  25  when  colder  tempera¬ 
tures  should  prevail.  On  September  29,  a  crew  of  four  was  left 
to  complete  camp  construction  and  prepare  the  runway.  On  the 
southward  trip,  the  ice  edge  was  crossed  at  73°7'N,  143°W. 

During  the  interim  period,  the  crew  at  the  ice  camp  re¬ 
ported  temperatures  as  high  as  0°C  with  rain!  The  same  storm 
removed  the  20  cm  of  snow  on  the  floe.  By  the  time  the  camp 
was  reactivated  on  October  19,  there  was,  however,  20  cm  of 
new  snow  on  the  floe.  Twin  Otter  service  from  Deadhorse  to 
the  ice  camp  was  established  on  October  25.  At  that  time  there 
was  still  at  least  150  km  of  open  water  north  of  Alaska. 

The  camp  averaged  about  fourteen  persons  from  Octo¬ 
ber  25  through  December  4,  with  many  groups  in  camp  for  only 
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Figure  2.  1993  was  an  unusual  ice  year. 
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a  week  or  two.  This  allowed  the  relatively  small  camp  staff 
(five)  to  support  over  thirty  science  personnel.  Most  of  the 
planned  scientific  experiments  were  accomplished  or  set  up  for 
winter-over  data  collection.  The  camp  staff  personnel  left  the 
camp  on  December  9,  with  the  intent  to  return  in  the  spring. 

ICE  FLOE  CONDITION 

The  ice  floe  for  the  SIMI  main  camp  was  like  an  un¬ 
even  grillage  in  the  sense  that  it  had  multi-year  ice  perforated 
with  many  holes  formed  by  melt  ponds.  It  was  similar  to  six 
other  floes  that  were  drilled  except  for  the  runway,  in  which 
much  of  the  ice  was  90  cm  thick.  The  runway  had  thick  multi¬ 
year  ice  on  three  sides  of  it;  however,  there  were  places  where 
the  runway  ice  was  only  30  cm  thick.  The  90-cm  ice  was 
apparently  formed  by  a  triple  raft.  On  further  examination,  it 
was  determined  that  this  "floe"  was  formed  from  two  pieces  of 
multi-year  ice  which  came  together  and  formed  the  runway  by 
rafting.  The  west  end  of  the  floe  (about  1/3)  was  not  truly  part 
of  the  "floe"  but  was  separated  pieces  of  multi-year  ice  embed¬ 
ded  in  thin  (30  cm)  ice  which  also  connected  them  to  the  floe. 
The  floe  boundary  was  defined  by  a  "line"  of  pressured  ice 
(blocks  pushed  up)  but  not  ridged  ice. 

One  of  the  reasons  for  a  fall  SIMI  program  was  to  see 
how  a  multi-year  floe  changed  through  the  fall,  winter,  and 
spring.  To  this  end,  this  year  was  good  for  the  experiment. 

Some  observations  of  note.  There  were  many  'ridged' 
areas  that  were  not  winter  ridges  of  first-year  ice  with  typical 
sail  and  keel.  These  ridges  were  formed  from  single  (maybe  a 
few)  blocks  of  old  ice  pushed  up  to  form  a  ridge  with  a  sail  but 
no  keel.  This  floe  was  considered  to  be  'near  the  ice  edge'  even 
though  it  was  at  75°N  because  it  was  in  such  an  extreme  condi¬ 
tion  of  melt 

ICE  DRIFT 

The  ice  drift  of  the  main  camp  is  shown  in  Figure  3 
through  January  1994.  With  some  data  gaps,  the  '+’  marks 
represent  the  main  camp  GPS  buoy  position  every  six  hours. 
Each  part  of  the  figure  shows  the  whole  trajectory  through  the 
end  of  the  month  indicated.  The  cluster  of  points  encircling 
both  ends  of  each  trajectory  indicate  the  positions  of  fourteen 
GPS  buoys  deployed  by  PMEL  around  the  main  camp.  Al¬ 
though  there  are  some  data  gaps,  the  relative  motion  of  the 
surrounding  buoys  are  small,  less  than  3  km.  These  data  will  be 
combined  with  ice  stress  data  as  described  in  Coon  et  al.  (1993) 
to  analyze  the  energy  dissipation  in  the  pack  ice  surrounding 
the  main  camp  for  the  development  of  a  pack  ice  anisotropic 
constitutive  law  (Coon  et  al,  1992). 

DISCUSSION  AND  FUTURE  PLANS 

Since  the  plan  was  to  reoccupy  the  fall  experiment  site 
in  the  spring,  we  had  to  take  into  account  the  expected  drift  of 
the  ice  when  choosing  the  location  in  the  fall.  Since  the  fall 
effort  was  supported  and  recovered  by  airplane,  the  site  needed 


to  be  within  Twin  Otter  range  from  shore.  For  this  reason,  the 
fall  experiment  site  was  chosen  to  be  500  km  northeast  of 
Deadhorse.  Because  of  the  light  ice  year,  however,  the  camp 
was  established  further  north  and  west  than  planned  to  remain 
500  km  from  Deadhorse. 

The  floe  for  the  camp  had  a  location  for  an  ice  runway; 
it  was  not  possible,  however,  to  land  a  Twin  Otter  on  September 
29  when  the  Polar  Star  was  to  leave.  Therefore,  a  crew  of  four 
was  left  at  the  ice  camp  with  plans  to  resume  full  operations  in 
late  October  1993  when  the  runway  ice  was  colder  and  stronger. 
SIMI  has  demonstrated  that  it  is  possible  to  establish  an  ice 
camp  (which  can  be  later  supported  by  aircraft)  in  the  fall  from 
an  icebreaker,  even  in  a  light  ice  year. 

As  of  this  writing,  the  ice  camp  has  drifted  to  74°40'N, 
155°W,  roughly  north  of  Barrow,  Alaska.  Should  the  camp 
remain  within  500  km  of  a  useful  airfield,  the  spring  field 
operations  will  continue  as  planned.  If,  however,  the  camp 
drifts  too  far  out  or  the  main  floe  is  unusable,  the  scientific  and 
ice  camp  equipment  will  be  recovered  in  a  prudent  manner  and 
a  new  ice  camp  will  be  established  north  of  Deadhorse,  Alaska, 
for  the  spring  field  operations. 
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NOMENCLATURE 

ASF  Alaska  SAR  Facility 

AUV  Autonomous  Underwater  Vehicle 

AVHRR  Advanced  Very  High  Resolution  Radiometer 

GPS  Global  Positioning  System 

NMC  National  Meteorological  Center 

NOGAPS  Navy  surface  wind  analysis  code 

NWS  National  Weather  Service 

ONR  Office  of  Naval  Research 

SAR  Synthetic  Aperture  Radar 

SIMI  Sea  Ice  Mechanics  Initiative 

SLAR  Side  Looking  Airborne  Radar 

SSMI  Special  Sensor  Microwave  Imager 

PMEL  National  Oceanic  and  Atmospheric  Admini¬ 

stration's  Pacific  Marine  Environmental  Labo¬ 
ratory 

USCG  United  States  Coast  Guard 
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Figure  3.  SIMI  Ice  Camp  Drift  (with  six-hour  data  intervals) 
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Dr.  Max  Coon  (P.I.),  Mr.  Skip  Echert,  and  Dr.  Stu  Knoke  of  Northwest  Research 
Associates  (NWRA). 

Scope  of  Model  and/or  Data  Set 

The  NWRA  team  performed  a  variety  of  experiments  during  the  SIMI  field  program:  ice- 
stress  measurements,  controlled-load  and  noise-generation  tests,  ice-on-ice  friction  tests, 
ridge-strength  tests,  and  a  multiyear  floe  evolution  study. 

Tee-Stress  Measurements  -  NWRA  successfully  deployed  flatjack,  fluid-filled  stress 
sensors  with  an  automatic  pressurization-check  feature  in  an  autonomous  station 
configuration.  We  installed  these  sea-ice-stress  "buoys"  at  four  sites;  each  buoy  measured 
ice  stress  in  four  directions,  from  which  we  will  derive  the  ice-stress  state  in  the  horizontal 
plane.  A  ridging  event  buried  the  first  buoy,  installed  near  Mt.  Curtin  at  the  West  SIMI 
Ice  Camp  (Map  2,  E7),  after  one  month.  The  second  buoy,  installed  about  two  miles 
north  of  West  Camp  on  Floe  4  (Map  1,  E4),  operated  continuously  for  over  nine  months. 
The  crack  that  flooded  portions  of  the  West  Camp  disabled  the  third  buoy  after  three 
months  (Map  2,  F5),  but  we  repaired  it  and  operated  it  for  an  additional  three  weeks  at  the 
East  Camp  (Map  4,  C6).  Table  1  lists  the  operational  lifetime  of  the  four  buoy 
installations  and  the  number  of  thermistors  used  to  measure  ice  temperatures.  The  buoys 
wrote  ice  stress  to  an  internal  storage  module  at  five-minute  intervals  and  transmitted  ten- 
minute  data  through  the  Argos  satellite  relay  system.  The  buoys  also  stored  and 
transmitted  temperature  readings  from  the  higher  (colder)  two  thermistors  on  the  same 
schedule  but  stored  and  transmitted  the  lower  (warmer)  thermistor  readings  and  a  battery 
voltage  every  six  hours. 


Table  1.  Ice-stress  buoys 


Buoy 

number 

Test 

site 

Start 

date 

End 

date 

Number  of 
thermistors 

1 

Mt.  Curtin 

29  Oct.  93 

4  Dec.  93 

4 

2 

Floe  4 

31  Oct.  93 

3  Aug.  94 

4 

3 

Lake  Andy 

22  Nov.  93 

25  Feb.  94 

8 

4 

East  camp  runway 

1  Apr.  94 

23  Apr.  94 

8 

In  support  of  the  sea-ice-stress  buoys,  NWRA  performed  various  other  ice-stress  sensor 
experiments  at  22  locations  near  the  two  ice  camps: 

•  Individual  stress  sensor  measurements, 

•  "High-data-rate"  stress  tests  with  data  collection  rates  to  4  Hz, 

•  "Sensor/ice  contact"  tests  to  check  contact  between  the  sensor  and  the  ice. 
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•  "Inclusion”  tests  to  calibrate  the  response  of  the  stress  gauge  to  the  ice  stress  produced 
by  an  air  jack  a  few  decimeters  away,  and 

•  Comparison  tests  with  other  types  of  ice-stress  and  strain  sensors. 

We  installed  and  tested  individual  stress  gauges  at  the  22  locations  listed  in  Table  2. 
Dataloggers  generally  collected  this  stress  and  temperature  data  at  five-minute  intervals. 
The  two  winter-over  sites  also  collected  temperature  data  hourly.  The  first  site  (Map  2, 
F4)  tested  the  installation  of  a  stress  gauge  in  the  low,  flat,  multiyear  hummocks  in  the 
early  autumn.  The  second  site  (Map  2,  F5)  tested  the  freeze-in  of  gauges  installed  at 
various  depths  in  autumn  first-year  ice;  the  lowest  was  half  in  the  water  at  installation. 
The  sensors  at  the  third  site  (Map  2,  F5)  were  installed  identically  to  the  stress  sensors  on 
the  buoys  to  obtain  inclusion  test  data  to  support  the  interpretation  of  the  buoy  data.  The 
fourth  site  (Map  2,  F4)  again  tested  the  installation  of  stress  gauges  in  the  low  flat 
hummocks  of  multiyear  ice  in  the  autumn.  The  fifth  site  (Map  2,  E3)  was  used  for  the 
comparison  tests  in  the  West  camp  runway  ice.  The  sixth  site  (Map  2,  E5  and  F5)  studied 
the  horizontal  and  vertical  variation  of  ice  stress  in  a  large  area  of  flat  first-year  ice.  Lake 
Andy.  (We  hypothesize  that  Lake  Andy  had  been  a  large,  bottomless  melt  pond  during 
the  previous  summer.)  The  seventh  site  (Map  2,  F4)  studied  the  local  variations  of  ice 
stress  near  a  hummock.  The  final  site  (Map  4,  C6)  was  installed  for  inclusion  test  data  to 
support  the  nearby  stress  buoy. 


Table  2.  Individual  ice-stress  gauge  sites 


Site 

number 

Test 

site 

Start 

date 

End 

date 

Number  of 
stress  gauges 

Number  of 
thermistors 

1 

Study  area 

23  Sept.  93 

27  Sept.  93 

1 

0 

2 

Lake  Andy 

27  Sept.  93 

16  Nov.  94 

5 

0 

3 

Lake  Andy 

15  Nov.  93 

28  Nov.  94 

3 

2 

4 

Study  area 

17  Nov.  93 

21  Nov.  93 

2 

2 

5 

West  camp  runway 

18  Nov.  93 

26  Nov.  93 

2 

2 

6 

Lake  Andy 

29  Nov.  93 

27  Feb.  94 

4 

2 

7 

Study  area 

30  Nov.  93 

17  Mar.  94 

3 

2 

8 

East  camp  runway 

5  Apr.  94 

23  Apr.  94 

2 

0 

"High-data-rate"  is  a  relative  term;  we  use  it  to  mean  data  rates  of  one  to  four  readings  per 
second  as  compared  to  our  standard  interval  of  five  minutes.  We  performed  high-data- 
rate  tests  with  one  to  four  stress  sensors  at  a  time,  looking  for  stress  signals  we  otherwise 
would  have  missed  due  to  sampling  rate.  We  also  obtained  ice-stress  data  at  one-second 
intervals  for  six  15-hour  periods  and  at  1/4-second  intervals  for  one  three-hour  period. 

Two  methods  were  used  to  evaluate  stress  sensor  installations:  sensor/ice  contact  tests 
and  inclusion  tests.  Turning  a  screw-plunger  assembly  added  additional  fluid  to  each  ice- 
stress  sensor  to  assure  good  contact  with  the  surrounding  ice.  The  additional  fluid  caused 
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either  an  abrupt  increase  in  fluid  pressure  in  the  sensor  or  no  response  at  all.  The  stress 
buoys  performed  these  tests  automatically  on  a  schedule,  which  proved  useful  for  timing 
the  melt-out  in  the  spring.  Each  stress  sensor  received  one  or  more  sensor/ice  contact 
tests;  the  response  varied  with  season,  ice  temperature,  and  ice  type.  If  the  sensor  showed 
an  increase  of  at  least  50  kPa  during  the  test,  we  considered  the  sensor  to  be  in  good 
contact  with  the  ice.  Sensors  that  failed  the  test  also  failed  to  respond  to  air  jack  stress 
and  natural  ice  stress. 

Air  jacks  were  used  to  perform  inclusion  tests,  which  were  stress  sensor  response  tests  at 
geophysical  stress  levels  (0  to  300  kPa)  and  geophysical  loading  rates  (1  to  100  kPa/min.). 
We  used  a  square,  flat,  air  jack  about  0.6  meter  on  a  side  to  apply  pressure  to  ice 
containing  one  or  two  sensors.  We  performed  "square-wave"  tests,  for  which  we  adjusted 
the  air  jack  pressure  in  an  on-off  cycle,  and  "ramp"  tests,  for  which  we  adjusted  the 
pressure  in  a  slow-increase,  slow-decrease  pattern.  The  square-wave  test  was  a  quick 
gauge  calibration,  but  the  ramp  test  simulated  realistic  loading  rates  based  on  our  prior 
field  data.  We  measured  the  stress-sensor  inclusion  factor  on  20  sensor  installations. 
Both  15-cm  and  20-cm  diameter  sensors  were  tested  using  square-wave  and  ramp  tests. 
On  the  square- wave  tests,  the  air  jack  pressure  was  varied  from  20  kPa  to  200  kPa  above 
ambient;  on  the  ramp  tests,  the  air  jack  pressure  was  100  to  200  kPa  above  ambient.  The 
square-wave  tests  typically  had  a  cycle  period  of  a  minute,  while  each  up  and  down  ramp 
test  took  from  one  to  three  hours. 

In  cooperation  with  Jackie  Richter-Menge  from  the  Cold  Regions  Research  and 
Engineering  Laboratory  (CRREL)  and  Dr.  Peter  Wadhams  from  the  Scott  Polar  Institute 
(SPI),  we  performed  a  test  series  to  provide  a  head-to-head  comparison  of  the  CRREL 
three-wire  stress  sensor  with  the  NWRA  flatjack  sensor  and  of  the  SPI  strain  gauge  with 
the  NWRA  flatjack  sensor.  We  ran  a  set  of  square-wave  tests  and  a  ramp  test.  In  each, 
the  air  jack  pressure  was  varied  from  ambient  to  about  60  kPa  above  ambient,  and  the 
resulting  stress  and  strain  sensor  readings  were  recorded. 

Controlled-Load  and  Noise-Generation  Tests  —  NWRA,  working  with  Dr.  Robert 
Pritchard  of  IceCasting,  Inc.  (ICI),  performed  over  30  controlled-load  tow  tests  and  40 
other  tests  to  generate  various  types  of  the  noise  representing  potential  Arctic  acoustic 
sources.  These  tests  were  performed  at  the  Football  Field  (Map  4,  D6),  East  Pond  (Map 
4,  C6),  South  Camp  (about  18:  km  south  of  East  Camp),  and  the  IOS  site  (Map  4,  B6). 
Dr.  Henrik  Schmidt's  group  from  MIT/WHOI  and/or  Dr.  David  Farmer's  group  from  the 
Institute  for  Ocean  Sciences  recorded  the  acoustic  signals  generated  by  many  of  these 
tests  with  hydrophones  and/or  geophones.  The  tow  tests  consisted  of  towing  a  block  of 
ice  on  partially  submerged  ice  to  simulate  rafting,  on  snow-covered  or  bare  ice  to  simulate 
over-ride,  and  under  young  ice  to  simulate  under-ride.  On  each  of  the  tow  tests,  the  tow 
load  was  usually  recorded  at  16  Hz  and  showed  irregular  slip-stick  oscillations.  We  also 
measured  the  tow  speed,  the  block  dimensions  and  weight,  and  the  snow  thickness  so  that 
we  could  calculate  the  minimum  and  maximum  sliding  friction  coefficient 
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The  other  tests  included  dropping  blocks  onto  ice,  flopping  blocks  onto  snow-covered  or 
bare  ice,  tipping  blocks  into  water,  tipping  blocks  underwater  to  float  up,  releasing  blocks 
underwater  to  bounce  on  the  bottom  of  the  ice,  and  breaking  cantilever  beams.  Table  3 
lists  these  acoustic  experiments. 


Table  3.  Acoustic  Source  Tests 


Block  Motion 

Number  of 
Tests 

Number  of  Unique  Blocks 
Used  for  Test  Series 

Flop  onto  bare  ice 

5 

5 

Flop  onto  snow 

7 

5 

Drop 

10 

9 

Splash 

2 

2 

Underwater  flop 

2 

2 

Underwater  bounce 

7 

2 

Cantilever  beam 

3 

3 

Ice-on-Ice  Friction  Tests  -  During  April  1993,  NWRA  ran  a  series  of  tests  to  measure  the 
ice-on-ice  friction  coefficients  on  "natural"  cracks  generated  in  the  sea  ice  near  Resolute, 
NWT.  These  tests  were  conducted  in  conjunction  with  Dr.  John  Dempsey  of  Clarkson 
University,  who  conducted  sea-ice  fracture  and  ice  characterization  tests.  All  of  these 
tests  were  performed  in  the  1.8-meter  thick,  first-year  sea  ice  in  Allen  Bay.  The  test 
parameters  are  listed  in  Table  4. 

Table  4.  Test  Parameters  for  Ice-on-Ice  Friction  Tests 


Normal  air 
pressure  (kPa) 

Number  of 
tests 

Total  face 
displacement  (mm) 

40 

7 

13-65 

200 

8 

41-54 

400 

4 

38-93 

Varied  15-170 

2 

77-89 

Ridge- Strength  Tests  -  NWRA  ran  strength  tests  on  first-year  ridged  ice  using  the  air  jack 
to  simultaneously  load  the  ice  and  measure  the  displacement.  In  November,  we  obtained 
underwater  video  and  one  inflation  sequence  of  the  air  jack  in  the  keel  at  the  bog  site 
(Map  2,  H5).  In  April,  we  characterized  the  sail  and  keel  geometry  of  a  ridge  formed 
from  first-year  ice,  inflated  the  air  jack  in  four  slots,  and  recorded  the  acoustic  emissions 
during  the  air-jack  tests  at  our  ridge  site  (Map  4,  A4).  The  ridge  configuration  in  the 
vicinity  of  the  four  air-jack  slots  was  characterized  by  drilling  12  holes  and  through 
observations  with  an  underwater  camera  (but  no  video  tape).  Dr.  Robert  Pritchard 
recorded  the  noise  using  Dr.  David  Farmer’s  hydrophone  equipment. 
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Multiyear  Floe  Evolution  Study  -  We  studied  the  evolution  of  a  60-by-60-meter  area  of 
multiyear  ice  near  the  SIMI  West  Camp  by  measuring  ice  thickness,  snow  cover,  ice 
temperature,  ice  salinity,  and  ice  stress  (Map  2,  F4).  Our  dataloggers  recorded  ice 
temperatures  at  three  depths  at  four  sites  every  hour  from  September  1993  to  March 
1994.  We  measured  ice  salinities  from  core  samples  taken  in  September.  We  also  made 
aligned  ice-stress  measurements  at  three  sites  from  November  1993  to  March  1994.  Ice 
thickness  and  snow  cover  were  measured  in  September  and  March  at  nine  locations  in  the 
study  area.  We  obtained  underwater  video  of  the  bottom  of  the  ice  in  September  through 
four  holes  in  the  ice  using  an  underwater  camera  mounted  on  a  pole.  Also,  U.S.  Coast 
Guard  divers  from  the  USCGC  Polar  Star  took  underwater  videos  of  the  underside  of  the 
ice  in  September  with  a  hand-held  camera. 


Findings  to  Date 

Ice-stress  Measurements  —  NWRA  buoys  measured  comparable  ice-stress  events  lasting 
three  to  nine  days  in  December  (days  350-358),  February  (days  33-35),  and  April  (days 
91-95),  as  can  be  seen  in  Figure  1.  Figure  1  shows  six-hour  samples  of  the  stress-sensor 
oil-pressure  data  transmitted  via  Argos  satellite  for  the  full  life  of  the  buoy,  including 
spikes  resulting  from  the  auto-pressurization  cycles.  There  were  no  major  stress  events 
after  mid-May,  and  the  ice  apparently  melted  away  from  the  sensors  in  early  July,  totally 
relieving  all  residual  ice  stress  around  the  sensors. 

We  found  comparable  stress  in  the  first-year  ice  on  opposite  sides  of  a  hummock  but  much 
lower  stresses  in  the  hummock  itself.  Some  hummocks  were  strong  ice  (the  stress  gauges 
responded  as  though  they  were  in  lake  ice),  but  others  were  very  soft  or  porous  ice  since 
the  stress  gauges  did  not  respond  to  sensor/ice  contact  tests  or  inclusion  tests.  We  also 
found  significant  stress  variations  across  the  length  of  Lake  Andy,  a  large  area  of  first-year 
ice  surrounded  by  multiyear  ice.  We  hypothesize  that  a  slip-stick  ice  motion  event  caused 
the  most  interesting  signal  which  we  recorded  during  the  high-data-rate  tests  (Figure  2). 
From  the  inclusion  tests,  we  found  that,  at  these  load  values  and  load  rates,  the  sensor 
response  is  significantly  lower  than  prior  laboratory  measurements  indicate,  implying  that 
ice-stress  levels  may  be  two  to  five  times  more  than  the  measured  sensor  fluid  pressure. 

Controlled-load  and  Noise-generation  Tests  -  Blocks  of  ice  sliding  over  or  under  flat  ice 
generally  exhibited  slip-stick  behavior  with  large  fluctuations  in  the  sliding  friction,  as 
illustrated  for  the  tow  load  history  in  Figure  3.  Figure  3  shows  the  first  12  seconds  of  Test 
8-13,  a  dry  over-ride  test;  the  signal  almost  immediately  settles  into  a  repetitive  slip-stick 
pattern.  The  load  fluctuations  on  the  ramp-ups  are  caused  by  cable  oscillations  since  the 
ice  block  was  stationary  during  those  times.  The  ice-deformation  mechanisms  that 
generated  more  acoustic  noise  were  blocks  dropping  or  flopping  onto  bare  ice,  blocks 
sliding  over  snow-covered  ice,  and  blocks  falling  into  the  water.  Quiet  ice-deformation 
mechanisms  included  blocks  sliding  or  bumping  underwater,  rafting  (lubricated  by  water), 
and  blocks  dropping  or  flopping  onto  snow-covered  ice. 
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Day  of  Year  (1993-1994) 


Figure  1.  Ice-stress  Buoy  2  (Floe  4)  data  via  Argos,  Nov.  1, 93  -  Aug.  3,  94 
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Figure  2.  Ice-stress  signal  recorded  near  East  Camp 
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Figure  3.  Tow  load  time  series  from  an  over-ride  test  (Test  8-13) 


Ice-on-Ice  Friction  Tests  —  In  the  configuration  tested  at  Resolute,  the  shear  force  to 
cause  shear  motion  along  a  fresh  crack  was  surprisingly  steady  and  smooth,  except  for  the 
initial  start-up  motion.  This  may  have  resulted  from  having  an  air  jack  provide  the  normal 
force  and  a  hydraulic  jack  provide  the  shear  force.  These  dynamics  are  very  different  from 
an  ice  block  being  pulled  with  a  winch.  The  measured  coefficients  of  friction  of  natural 
cracks  in  first-year  sea  ice  are  consistent  with  the  friction  angle  of  pack  ice  calculated  from 
ice-stress  measurements  made  during  the  CEAREX  program. 

Ridge-strength  Tests  —  The  strength  of  the  fused  ice  in  the  ridge  was  highly  variable, 
ranging  from  essentially  solid  ice  near  the  surface  to  loose  pieces  floating  underwater  at  or 
near  their  melting  point  temperature.  Figure  4  shows  an  intermediate  case  of  air  jack 
pressure  versus  time.  We  interpret  the  down  steps  during  the  long  rise  as  ice-to-ice 
contacts  breaking  and  slipping.  We  stopped  the  displacement  of  the  air  jack  when  it  had 
reached  maximum  stroke.  We  caused  the  large  down  step  at  the  end  by  releasing  the 
pressure  in  the  air  jack.  These  tests  suggested  that  a  better  test  of  ridge  keel  strength 
would  require  a  larger  volume  displacement,  on  the  order  of  8  cubic  meters  or  more. 
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This  panel  met  Wednesday,  April  26,  1995,  in  the  morning.  The  discussion  of  the  100  km 
scale  (pack  ice)  centered  around  the  anisotropic  behavior  of  the  pack  ice  as  indicated  and  driven 
by  the  formation,  deformation,  and  closure  of  leads. 

GOALS 

The  SIMI  goals  are  reported  in  Appendix  B,  Vol.  1,  of  these  Proceedings. 

•  Understand  sea  ice  constitutive  laws  and  fracture  mechanics  over  the  full  range  of 
geophysical  scales  and  determine  the  scaled  responses  to  applied  external  forces 

•  Develop  physically  based  constitutive  and  fracture  models 

One  of  the  aims  of  the  Sea  Ice  Mechanics  Initiative  (SEMI)  was  to  bridge  scales  in  sea  ice 
mechanics.  This  can  and  has  been  done  by  considering  pack  ice  on  a  scale  from  3-100  km  as 
made  up  of  two  types  of  ice  as  shown  in  Figure  1.  The  two  types  of  ice  are  not  necessarily 
represented  by  presently  used  nomenclature.  The  scene  in  Figure  1  is  intended  to  represent  a  sea 
ice  element  system  with  an  active  lead  in  it.  A  lead  system  is  made  of  parallel  leads  that  are 
shown  as  one  lead  in  the  element  of  Figure  1. 

So,  what  is  the  Type-1  material  in  Figure  1?  The  ice  in  material  Type-1  includes  multiyear 
ice,  first-year  ice,  ridged-ice,  rafted-ice,  and  inactive  refrozen  leads.  The  material  has  a  wide 
distribution  of  sea  ice  thickness  and  is  made  up  of  closely-fitting,  angular  particles,  grains,  or  cells, 
often  diamond  shaped  in  the  winter.  Type-1  material  is  deformable,  controlled  by  a  2-D  Coulomb 
behavior  with  no  significant  thickness  redistribution  due  to  deformation.  The  particles,  grains,  or 
cells  that  make  up  Type-1  material  are  of  order  100  meters,  and  they  have  been  formed  by  out-of¬ 
plane  loadings  such  as  thermal  loading  or  isostatic  imbalance.  These  loadings  have  led  to  cracks 
which  have  produced  the  grains  or  cells,  and  the  cracks  have  only  partially  healed.  This  cracking 
and  rehealing  are  the  small-scale  mechanisms  which  cause  the  Type-1  material  to  have  its 
properties.  When  this  material  deforms,  it  produces  new  areas  for  Type-2  material  (leads). 
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Fig.  1.  Element  of  pack  ice  with  a  refrozen  lead. 


Type-2  material  is  what  we  typically  think  of  as  lead  ice  or  rafted-lead  ice  or  ridged  ice.  It 
is  the  ice  that  was  formed  in  the  openings  made  in  Type-1  material  and  then  deformed.  Therefore, 
the  constituents  of  Type-2  ice  are  open  water,  thin  new  ice,  thick  new  ice,  rafted  ice,  ridged  ice, 
and  floating  broken  pieces  of  the  lead  ice.  This  material  has  a  thickness  distribution  which  is 
associated  with  a  direction,  0  in  Figure  1,  and  the  redistribution  is  caused  by  deformation.  This 
material  has  two  discrete  kinds  of  behavior,  both  of  which  are  the  small-scale  mechanisms.  One  is 
the  failure  or  initialization  of  dynamics  of  the  Type-2  ice.  This  initialization  is  controlled  by  small- 
scale  mechanisms  such  as  the  strength  of  an  intact  new  ice  sheet  in  tension,  compression,  or  shear. 
This  type  of  behavior  is  affected  by  ice  thickness,  temperatures,  salinity,  grain  size,  and  strain 
rates.  Another  type  of  failure  of  this  material  is  controlled  by  buckling  instabilities  which  might  be 
elastic  under  short-term  loadings  or  visco-elastic  under  longer  duration  loadings.  A  third 
intermediate  type  of  initialization  of  failure  for  this  material  is  the  buckling  of  individual  (ridged) 
blocks  that  might  be  floating  in  the  lead.  The  second  kind  of  behavior  here  is  that  of  the 
mechanical  redistribution  of  the  ice,  and  this  is  then  controlled  by  either  the  rafting  process  which 
depends  on  ice  thickness  and  friction  or  the  ridging  process  which  is  controlled  by  thickness, 
friction,  and  change  in  energy  state. 

Type-1  ice  is  isotropic  and  has  a  failure  surface  which  can  be  expressed  in  terms  of  stress 
resultant  invariants  which  change  very  slowly  and  may  be  approximated  as  constant  in  time. 

The  Type-2  ice  is  oriented,  and  its  failure  surface  cannot  be  expressed  in  terms  of 
invariants  and  its  failure  surface  varies  over  shorter  time  scales  (hours)  driven  by  the  ridging  and 
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rafting  pressures.  Pack  ice  is  a  combination  of  Type  1  and  Type  2  and  has  an  anisotropic  behavior 
whose  properties  can  change  over  time  periods  of  hours. 

Coon  et  al.  [1992]  introduced  a  method  for  incorporating  the  behavior  of  individual 
leading,  rafting,  and  ridging  events  into  the  large-scale  constitutive  behavior  of  pack  ice.  This  was 
not  expressed  in  terms  of  Type  1  and  Type  2  materials,  but  that  is  what  it  represents.  The 
concept  is  that  the  large-scale  stress  state  must  lie  within  all  of  the  yield  surfaces  that  describe 
failure  of  individually  resolved  ice  features.  Ice  conditions  are  described  in  terms  of  an  onented 
thickness  distribution  including  the  thickness  of  ice  in  leads.  In  resolving  behavior  on  scales  from 
one  to  hundreds  of  kilometers,  processes  of  leading,  rafting,  and  ridging  that  describe  and  control 
the  behavior  of  sea  ice  are  explicitly  considered.  Previous  models  [e.g.,  Coon  et  al.,  1974]  are 
applicable  only  over  scales  large  enough  that  the  multiplicity  of  ice  features  give  rise  to  an 

isotropic  response. 

SIMI  included  nested  arrays  of  drifting  buoys  to  measure  deformation  on  different  scales. 
The  large,  pack  scale  Arctic  buoy  array,  with  resolution  of  order  tens  of  kilometers,  is  not  useful 
for  locating  individual  events.  The  floe  cluster  array  of  GPS  and  the  stress  buoys  was  deployed  at 
roughly  3  km  intervals  over  a  uniform  grid  15  to  20  km  on  a  side,  allowing  the  deformation  of 
individual  events  to  be  tracked.  Vol.  1  of  these  Proceedings  discusses  the  data  taken. 

The  SIMI  data  collection  plan  for  this  scale  is  given  in  Appendix  B  of  Vol.  1  which  said: 
"An  anisotropic  model  enables  measurements  of  stress,  deformation,  and  ice  conditions  associated 
with  individual  leading,  rafting,  and  ridging  events  to  be  used  as  components  of  a  large-scale 
model.  During  SIMI,  these  variables  will  be  measured  during  specific  events.  At  times,  isolated 
events  with  length  scales  of  order  100  m  will  likely  occur  within  the  floe  cluster  array  while  the 
remainder  of  the  array  will  be  nearly  rigid.  Such  events  will  be  described  by  observations  on  the 
scale  of  3  to  15  km  by  the  deformation  of  the  array  and  by  direct  ground  observations.  Airborne 
imagery  of  the  region  near  the  end  of  October  will  help  identify  earlier  events  and  provide  a 
baseline  for  comparison  in  spring.  Ice  motion  vectors  will  be  derived  from  SAR  and  AVHRR 
imagery."  A  test  of  this  model  and  data  is  given  by  Coon  et  al.  [1995]  in  this  Volume.  For  the 
first  time,  simultaneous  measurements  of  sea  ice  stress  on  both  sides  of  a  lead,  ice  motion  from 
sequential  SAR  imagery,  and  ice  motion  from  drifting  GPS/Argos  buoys  are  compared.  The  lead, 
which  appeared  early  in  Day  35  (Feb.  4),  1994,  was  observed,  relative  motions  were  calculated 
across  and  along  the  lead  direction,  and  stresses  were  transformed  to  this  same  coordinate  system. 
The  lead  orientation  and  approximate  location  were  determined  from  SAR  images  of  the  pack  ice 
surrounding  the  SIMI  camp,  as  well  as  net  ice  motion  at  3-day  intervals.  The  ice  motion  around 
the  lead  was  provided  in  greater  detail  by  an  array  of  drifting  GPS/Argos  buoys  producing  hourly 
positions.  Sea  ice  stress  was  measured  on  each  side  of  the  lead  using  flapjack  sensors.  The  ice 
motion  normal  to  the  lead  and  the  ice  stress  normal  to  the  lead  are  compared  as  well  as  the  ice 
shear  motion  and  the  ice  shear  stress. 

Shortly  before  the  lead  opening  motion  is  apparent  in  the  ice  motion  data,  the  ice  stress 
state  abruptly  changed  from  biaxial  compression  to  uniaxial  compression  with  traction  on  the  lead 
near  zero.  We  thus  measured  the  stress  conditions  at  lead  creation  as  sheared  compression,  there 
was  no  tension  at  the  time.  The  Type-1  ice,  made  up  of  closely  fitting  angular  particles,  begins  to 
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fail  and  create  a  lead  under  shear  and  compressive  stress  just  as  sand  can  expand  under  shear  and 
compressive  stress  by  going  from  a  closely  packed  configuration  to  a  loosely  packed 
configuration.  As  the  lead  continues  to  open  the  stress  normal  to  the  lead  is  zero,  the  stress 
parallel  to  the  lead  is  compression,  and  there  are  no  tensile  stresses.  Excellent  correlation  was 
observed  between  ice  motions  measured  with  SAR  products  and  drifting  buoys.  The  SAR 
products  provide  a  more  precise  definition  of  the  orientation  of  the  lead.  It  was  also  concluded 
that  a  properly  installed  20-cm  stress  gauge  in  a  single  floe  can  provide  a  reasonable  measurement 
of  the  large-scale  (50-km)  pack  ice  stress  resultant. 

The  needs  of  two  user  groups  were  expressed  during  the  meeting.  Ruth  Preller  expressed 
needs  from  the  National  Ice  Center  for  methods  of  predicting  leads  (interpreted  to  mean  direction 
and  area)  and  a  need  for  smaller  scale,  i.e.,  a  few  kilometer  resolution.  Ruth  went  on  to  describe 
the  plans  for  the  Navy  in  its  Polar  Ice  Prediction  System  (PIPS)  to  have  a  model  which  would 
cover  all  areas  in  the  northern  hemisphere  covered  by  sea  ice  with  a  resolution  down  to  20 
kilometers.  There  was  discussion  and  concern  about  one  model  being  acceptable  for  the  multi¬ 
year  pack  ice  in  the  central  Beaufort  Sea  and  the  marginal  sea  ice  zones.  One  model  can  be  made 
to  work  if  the  ice  properties  in  the  model  vary  sufficiently  such  that  the  multi-year  and  marginal 
regions  can  be  described. 

The  needs  of  the  oil  industry  were  articulated  by  Albert  Wang.  Albert  explained  that  for 
the  oil  industry  real-time  information  was  important.  Forecasting  of  the  pressure  in  the  ice  was  a 
key  as  well  as  the  ice  conditions  including  leads  and  thickness  distribution.  There  is  a  desire  to 
use  large-scale  ice  models  to  get  large-scale  forces  that  would  be  valuable  for  loads  on  structures. 
There  is  a  need  for  ground  truthing  of  the  models  such  that  confidence  can  be  developed.  For  the 
oil  industry  and  offshore  operations,  a  three-to-four-day  forecast  of  extreme  events  would  be 
desirable  but  a  one-day  forecast  could  lead  to  a  form  of  orderly  shutdown. 

Kate  Hedstrom  discussed  ice  modeling  work  at  Rutgers  University  where  they  are 
coupling  an  ocean  model  to  an  ice-atmosphere  climate  model  (ARCSyM,  from  U.  of  AK, 
Fairbanks).  They  also  have  done  ice-ocean  modeling  for  use  in  oil  spill  risk  assessment  studies  for 
MMS.  A  third  project  involves  the  development  of  a  spectral  finite  element  ocean  model  to  be 
used  on  a  global,  variable-resolution  grid.  Several  of  these  projects  require  an  ice  model  with 
better  small-scale  ice  physics.  For  additional  information,  see  Hedstrom  (1995)  in  this  Volume. 

Igor  Appel  of  Fairweather  Forecasting  discussed  a  model  they  use.  The  approach  used  is 
described  in  Appel  (1990)  and  results  of  the  model  are  presented  in  Appel  and  Samuels  (1995). 
In  this  model  the  anisotropic  behavior  (interpreted  as  leads)  is  described  by  a  viscous  model.  The 
viscosity  depends  on  the  value  of  principal  strain  rate.  In  this  model  the  leads  would  open  across 
the  direction  of  the  principal  extension  strain  rate  direction.  This  is  a  difference  between  this 
model  and  that  of  Coon  et  al.  (1992)  when  there  would  be  an  angle  between  the  direction  of  the 
lead  and  principal  strain  rate  direction.  In  the  work  of  Cunningham  et  al  (1994),  there  is  evidence 
of  such  an  angle. 
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There  is  also  a  model  by  Moritz  (1995)  that  accounts  for  leads  in  a  statistical  manner.  In 
this  model  and  the  viscous  model  discussed  by  Appel  there  is  no  direct  way  to  bring  in  the  small 
scale  ice  mechanics  as  there  is  with  the  Coon  et  al.  model. 

The  need  to  be  able  to  account  for  leads  and  their  behavior  is  a  theme  that  comes  up  in 
many  of  the  other  panel  discussions  on  ice  behavior  and  applications  in  this  Volume.  Accounting 
for  this  behavior  and  using  small  scale  mechanics  to  do  so  will  meet  the  goals  of  SIMI  as  stated 

above. 


The  discussion  of  the  panel  is  summarized  by  the  two  lists  below.  The  first  item  in  each 
list  is  that  of  most  importance. 

LIST  OF  THINGS  WE  HAVE  LEARNED: 

1)  From  GPS  buoys,  stress  buoys,  and  SAR,  the  importance  of  lead  formation  and 
direction,  on  the  large  scale  behavior  of  pack  ice,  has  been  verified.  New  leads  have 
been  shown  to  be  created  from  states  of  compression  and  shear  stress  as  measured 
on  the  scale  of  a  floe  but  representing  stress  on  tens  of  km.  Also  when  leads  are 
open  the  measured  stress  state  is  compression  parallel  to  the  lead  and  zero 
perpendicular  to  the  lead.  These  findings  are  new  and  very  important  to  the  large 
(100  km)  behavior  of  the  ice. 

2)  We  have  developed  a  formulation  for  an  anisotropic  behavior  (constitutive  law)  for 
pack  ice  which  can  account  for  creation  and  direction  of  leads 

3)  Examples  of  improved  ice-ocean  models  are  northern  hemisphere  PIPS,  ARCSyM 
(U.  of  AK,  Fairbanks  and  Rutgers  University),  etc.  Desirable  features  are: 

•  Higher  resolution  coverage, 

•  Ability  to  assimilate  satellite  data,  and 

•  Coupling  to  the  atmosphere. 

WHAT  NEEDS  TO  BE  DONE  NEXT: 

1)  Need  formulation  and  verification  of  model  against  deformation  (strain  rates),  leads 
(direction  and  deformation),  and  ice  concentration  using  buoy  and  satellite  data. 

2)  Need  to  put  new  rheology  into  improved  models. 

3)  Need  to  use  SIMI  data  to  quantify  the  terms  in  the  new  constitutive  law. 

4)  Need  to  incorporate  results  of  smaller-scale  models  in  large-scale  model. 

5)  More  aware  of  user  needs. 
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The  10  km  scale  is  critical  to  understanding  and  modeling  spatial  and  temporal  variability 
of  sea-ice  dynamics  (Figure  1).  Applications  are  ice  forecasting,  ice  loads  on  structures,  and 
pollution  transport.  At  this  scale  there  are  two  important  transitions  as  the  scale  increases,  from 
understanding  ice  as  discrete  floes  to  a  continuum  medium  and  from  a  stochastic  view  of  ice 
processes  to  deterministic  parameterizations.  In  the  early  1970s,  during  AIDJEX,  researchers 
developed  the  concept  of  sea  ice  on  the  regional  scale  as  a  strain-hardening  plastic  material.  This 
remains  the  prevailing  approach.  In  a  plastic  response,  as  the  external  forcing  increases,  the  stress 
in  the  ice  pack  will  increase  up  to  a  certain  failure  limit,  at  which  point  the  ice  will  deform.  As  the 
pack  ice  converges,  its  strength/compactness/thickness  increases  and  thus  its  failure  limit  increases 
and  the  pack  "hardens".  This  process  results  in  the  development  of  regional  stresses  in  the  ice 
cover,  which  probably  relates  to  sea-ice  features  such  as  the  long,  linear  leads  seen  on  AVHRR 
images.  The  question  exists  as  to  whether  these  sea-ice  dynamics  assumptions  are  valid  based  on 
measurements  at  the  10  km  scale,  where  we  consider  a  cluster  of  floes. 

The  essential  processes  that  control  ice-dynamics  behavior  are  leading,  rafting,  and 
ridging.  The  present  generation  of  isotropic  sea-ice  models  were  originally  developed  to  describe 
ice  dynamics  on  scales  averaged  over  100  km  and  larger.  Researchers  hypothesized  that  this  scale 
was  large  enough  to  include  the  effects  of  many  leads  and/or  ridges  at  random  orientation  within 
the  sampling  scale.  There  are  times,  however,  when  an  individual  large  feature  can  dominate 
behavior  even  on  these  scales;  for  example,  if  one  huge  lead  opens  or  narrow,  linear 
discontinuities  propagate  in  from  boundaries.  The  isotropic  model  may  adequately  represent 
behavior  averaged  over  10  km,  such  as  during  conditions  just  after  freeze-up  when  the  floes  are 
small  (<1  km)  with  randomly  oriented  thin  ice  between  them.  As  the  scale  is  reduced  from  100 
km,  however,  the  likelihood  of  encountering  strongly  oriented  features  increases,  particularly  for 
winter  pack  ice  conditions.  Thus,  to  resolve  regional  sea-ice  dynamics  behavior  it  appears 
necessary  to  account  for  anisotropy.  Deriving  more  accurate  solutions  by  reducing  grid  size  in 
existing  models  is  beneficial,  but  these  solutions  must  still  be  interpreted  as  averages  over  the 
larger  scale.  A  major  conclusion  of  SIMI  research  is  that  it  is  important  to  understand  and  model 
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anisotropic  sea-ice  processes  on  the  10  km  scale  to  correctly  interpret  processes  on  the  100  km 
scale.  The  final  summary  discussion  of  the  meeting  noted  a  need  to  merge  the  10  and  100  km 
scales  in  terms  of  understanding  ice  mechanics. 

The  findings  of  LEADEX  and  SIMI  Experiments  indicate  a  mismatch  in  scales  between 
the  forcing  on  the  ice  field,  primarily  wind  and  coastal  orientation,  and  pack  ice  itself. 

Atmospheric  storms  have  large  spatial  scales  (order  of  800  km),  and  small  temporal  scales  (order 
of  days).  Regional  sea-ice  motion  has  a  temporal  scale  of  a  week  or  longer  and  a  spatial  scale  of 
kilometers.  Therefore,  the  cluster  scale  response  of  the  sea  ice  is  not  directly  forced  by  a  similar 
cluster  scale  variability  in  the  wind.  The  cluster  scale  response  of  the  ice  is  determined  by  the 
aggregate  mechanical  properties  of  the  ice  field  in  response  to  nearly  uniform  wind  fields,  acting 
over  3-5  days,  before  or  after  passages  of  a  storm  front.  Figure  2  shows  a  S  AR  motion  vector 
plot  for  a  shear/opening  event  near  day  35  in  early  February  1994.  The  vectors  are  5  km  apart. 
There  are  large  regions  (plates?)  which  behave  with  apparent  solid  body  motion  separated  by  long 
(100  km),  narrow  (<5  km)  shear  lines.  The  sea-ice  problem  on  the  10  km  scale  is  the  aggregate 
problem:  describe  a  material  which  under  uniform  winds,  subject  to  boundary  forcing,  will 
respond  by  the  formation  of  plates,  which  are  large  compared  to  the  floe  size  and  separated  by 
discontinuities.  Also,  what  pack-ice  mechanical  processes  result  in  this  material  description  when 
aggregated?  A  goal  of  SIMI  and  post  SIMI  research  is  to  describe  and  model  this  internal  sea-ice 
mechanics. 

During  the  SIMI  field  program,  researchers  examined  sea  ice  dynamics  concepts  at  a 
variety  of  spatial  and  temporal  scales.  The  10  km  studies  concentrated  on  identifying  and 
studying  cluster-scale  events  during  an  extended  period  with  both  Synthetic  Aperture  Radar 
(SAR)  data  and  Global  Positioning  System  (GPS)  buoy  data  (Days  23  to  75).  Three  technical 
advances  facilitated  this  work  at  the  10  km  scale.  The  first  was  the  availability  of  ice  motion 
vectors  derived  from  SAR  data  processed  at  the  Alaska  SAR  Facility.  This  product  uses 
sequential  SAR  images  to  track  ice  motion  at  a  5  km  resolution  over  an  area  of  80  X  80  km  per 
product.  Taken  together,  these  products  can  define  regional  (80  x  400  km)  deformation  fields  at 
a  three  day  sampling  interval.  The  second  technical  advance  was  that  the  accuracy  of  the  GPS 
allowed  researchers  to  resolve  relative  motions  between  adjacent  floes  in  the  vicinity  of  the  SIMI 
ice  camps.  These  GPS  measurements  provided  a  comparison  between  the  differential  GPS  ice 
motion  with  the  SAR-derived  ice  motion  fields.  The  third  advance  was  the  deployment  of  arrays 
of  stress  sensors  that  had  been  developed  and  tested  on  earlier  experiments.  Two  major  events 
were  studied  during  the  SIMI  period.  One  significant  event  occurred  in  early  February  1994 
(days  38  to  41)  with  uniform  northerly  winds  in  which  a  convergence/shear  event  propagated 
northward  away  from  the  coast,  reaching  the  location  of  the  SIMI  ice  camp  400  km  north  of 
Alaska  in  3  days.  A  major  opening/shear  event  occurred  between  days  63  to  74. 

Preliminary  results  from  a  20  km  stress  and  deformation  array,  centered  about  the  west 
camp  floe,  suggest  a  correspondence  between  stress  and  deformation  at  the  10  km  scale,  but  not 
at  the  5  km  scale.  Further,  a  comparison  of  in  situ  stresses  on  the  central  ice  floe  show  a  strong 
linear  correlation  in  time,  indicating  that  the  floe  responds  to  regional  (10  km),  rather  than  local  (1 
km)  deformation  events.  Together  these  results  imply  that  there  is  a  lower  bound,  between  10 
and  5  km,  for  considering  the  ice  floes  in  a  region  as  a  continuum.  There  is  a  higher  frequency  of 
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major  stress  and  regional  deformation  events  as  compared  to  large-scale  (50  to  100  km) 
deformation  events,  defined  by  the  formation  of  shear  line  patterns  observed  in  SAR  images.  This 
may  indicate  several  levels  of  continuum  behavior,  with  the  basic  unit  being  a  cluster  of  floes  on 
the  scale  of  10  km.  Once  the  floes  within  these  clusters  have  moved  into  a  close-packing,  the 
clusters  themselves  can  lock  together  in  response  to  a  forcing  system.  This  structure  can  then 
form  regional  plates  separated  by  discontinuities. 

There  are  several  promising  approaches  to  describe  the  ice  interaction  and  the  spatial 
evolution  of  stress  at  the  10  km  scale.  One  approach  is  to  model  the  pack  ice  as  a  classical 
granular  material  with  an  effective  stiffness  and  strength.  Results  of  such  models  suggest  that 
variation  in  size,  shape,  and  orientation  of  the  granular  elements,  or  ice  floes,  are  important  to 
producing  fracture-like  features.  Application  of  this  technique  could  lead  to  the  development  of 
an  effective  (i.e.,  large-scale  parameterization)  constitutive  response  (i.e.,  stress/strain  rate 
relationship).  A  very  promising  approach  is  to  directly  observe  the  aggregation  properties  of  pack 
ice  through  the  analyses  of  SAR  imagery.  Two  state  parameters,  texture  and  structure,  can  be 
derived  from  measures  of  nearest  neighbor  pixel  intensity  differences  on  the  SAR  data.  Sharp 
changes  in  the  slopes  of  observational  points  on  a  texture/structure  axis  plot  have  an  analog  to 
phase  transitions  and  may  correspond  to  boundaries  between  different  ice  types.  There  is  a 
relationship  of  "information  content"  between  these  analyses  and  similar  1-D  analysis  of  ice 
thickness  estimates  from  sonar.  There  are  subsets  of  SAR  pixels  which  possess  a  similarity 
property  which  are  multi-fractal.  These  sets  display  a  critical  connectivity  tentatively  associated 
with  ice  fracture  nucleation,  i.e.,  similar  to  avalanche  phenomena.  A  value  exists  of  the  nearest 
neighbor  probability  measure  where  pixels  are  mobilized  into  a  space-spanning  path  called  the 
percolation  cluster.  For  large-scale  compression  of  a  mixture  of  stiff,  brittle  inclusions  (i.e., 
multi-year  ice  floes?),  and  the  percolation  cluster  (i.e.,  the  rubble  field?),  the  least  unconsolidated 
part  absorbs  the  strain.  In  materials  other  than  ice,  observations  show  that  stress  is  transmitted 
along  such  poorly  bonded  paths.  If  the  analog  of  sea  ice  to  other  stochastically  described 
materials  is  consistent,  then  the  percolation  cluster  hypothesis  associated  with  FY  rubble/ndges, 
may  relate  to  the  load-bearing  capacity  of  the  ice  field. 

A  third  approach  to  model  pack  ice  at  the  10  km  scale  is  the  oriented  pack  ice  constitutive 
law.  This  law  includes  the  creation  of  new  leads  in  thick  ice  and  the  effects  of  leads  and  their 
orientations  on  the  strength  of  the  ice  cover.  The  behavior  requires  an  oriented  ice  thickness 
distribution  to  account  for  the  growth  and  deformation  of  ice  in  refrozen  leads.  The  shear  and 
compressive  strength  of  each  lead  are  functions  of  the  lead  s  ice  thickness  distribution.  This 
approach  predicts  the  orientation  of  lead  systems  and  accounts  for  the  highly  anisotropic  strength 
of  the  pack  ice  resulting  from  lead  systems. 

In  summary,  wind  and  coastal  boundaries  have  large  spatial  scales  (order  of  500  km)  while 
ice  response  is  at  all  scales,  with  a  tendency  to  form  long,  linear  discontinuities  at  the  10-100  km 
scale.  There  are  distinct  patterns  apparent  in  currently  active  leads  and  ridges  which  depend  on 
boundary  conditions  and  the  type  of  forcing.  Historical  ridges  and  leads  may  appear  more 
random.  There  appears  to  be  a  scale  discontinuity  in  understanding  between  the  regional  and 
smaller  scales.  At  the  0.1-100  m  modelers  can  provide  scaling  laws  for  strength.  However,  at  the 
1  km  scale  much  of  the  work  on  ridging  is  focused  on  linear  compression,  where  evidence  from 
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SAR  analyses  suggests  that  shear  is  a  major  yield  mechanism  and  that  failures  follow  lightening 
bolt,  fault-like  structures.  Studies  of  strengths  and  failures  of  1  km  structures  need  to  be 
coordinated  with  the  horizontal  2-D  percolation  ideas  discussed  above.  It  is  important  to 
understand  how  stress  propagates  through  the  pack  based  on  coastal  orientation,  wind  forcing  and 
ice  property  distribution  on  a  10  km  scale. 

The  group  stated  that  the  most  important  SIMI  accomplishment  was  to  acquire  the  first 
major  data  set  of  stress  and  strain  rate  of  a  pack  ice  cover  on  the  10  km  scale,  and  to  correlate  this 
with  results  from  both  smaller  and  larger  scale  observations.  There  was  good  correlation  between 
stress  and  strain  rate  at  the  10  km  scale,  but  not  at  small  scales,  even  though  the  stress  and  strain 
measurements  existed  down  to  a  1  km  scale.  We  believe  that  ice  can  be  treated  as  a  continuum 
at  10  km  as  long  as  the  medium  allows  plate/discontinuity  type  structures. 

The  group  recommended  that  the  main  goal  for  future  research  is  to  understand  what 
mechanical  processes  control  the  spacing  of  discontinuities  (slip  lines,  leads,  etc.).  On  AVHRR 
imagery  there  seems  to  be  a  preference  for  a  15  km  and  a  35-40  km  separation  scale.  There  is  a 
need  for  more  case  studies  of  the  stress-strain  relationship  for  ice  deformation  events  and 
modeling  of  buckling,  rafting,  and  ridging  in  a  3-D  context  to  characterize  the  anisotropic  strength 
parameters  of  pack  ice  in  terms  of  an  oriented  thickness  distribution. 

Figure  1.  Distribution  of  large  floes  in  the  vicinity  of  the  SIMI  camp.  Space  between  the  larger 
floes  was  filled  with  rubble  and  small  floes. 

Figure  2.  Relative  motion  of  an  ice  field  in  the  vicinity  of  the  Ice  Camp  on  day  35-38, 1994  from 
5  km  spacing  SAR  motion  data.  The  large  arrow  is  north.  (Figure  courtesy  of  H.  Stem). 
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Map  1.  Vicinity  of  SIMI  West  Camp,  September  1993. 
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The  large-scale  composition  of  the  ice  pack  is  the  aggregate  of  the  basic  1-km-scale 
building  blocks,  discrete  parcels  of  multiyear  and  first-year  ice,  and  leads.  Similarly,  large-scale 
deformation  is  the  aggregate  of  the  basic  1-km-scale  processes  of  pressure  ridging,  rafting,  m- 
plane  sliding,  and  opening,  which  redistribute  energy  from  larger  scales  to  smaller  scales.  By 
determining  stress  levels  in  the  surrounding  pack,  1-km-scale  processes  help  define  design  loads 
on  offshore  structures.  The  1-km-scale  work  in  the  SIMI  program  should  be  assessed  on  the 

basis  of  the  following  criteria: 

•  Contribution  to  understanding  of  the  processes  of  ridging,  rafting,  floe 
splitting,  and  in-plane  sliding  between  ice  parcels; 

•  Contribution  to  determination  of  realistic  estimates  of  forces  on  offshore 
structures; 

•  Contribution  to  large-scale  modeling  efforts. 


1-KM-SCALE  PROCESSES 


Pressure  Ridging 

A  realistic  computer  model  of  the  pressure  ridging  process  [Hopkins]  was  developed  in 
SIMI.  This  two-dimensional  ridging  model  is  based  on  a  concept  of  ridge  growth  in  which  an 
intact  sheet  of  thin  ice  is  driven  against  a  thick  ice  floe.  The  key  assumption  is  that  the  ice  sheet 
breaks  in  flexure.  The  important  features  of  the  ridging  model  are:  a  dynamic  linear  viscous- 
elastic  model  of  a  floating  ice  sheet,  flexural  failure  (including  buckling)  of  the  ice  sheets,  blocks 
broken  from  the  parent  sheet  at  points  where  tensile  stress  exceeds  strength,  secondary  flexural 
breakage  of  rubble  blocks,  inelastic  contacts  between  rubble  blocks,  frictional  sliding  contacts 
between  blocks,  separate  friction  coefficients  for  submerged  and  above  water  contacts,  buoyancy 
of  the  ice  sheets  and  rubble,  and  water  drag.  Simulations  were  performed  with  the  model  to 
determine  the  energetics  and  dynamics  of  the  ridging  process.  Major  unknowns  in  modeling  the 
ridging  process  are  the  various  friction  coefficients. 
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In-situ  friction  measurements  between  blocks  of  sea  ice  were  made  during  SIMI  by  Coon 
et  al.  and  Pritchard.  Preliminary  analysis  [Pritchard]  of  data  from  dry  override  tests  points  to  a 
coefficient  of  sliding  friction  of  approximately  0.8+-0.2  and  a  static  friction  coefficient  of 
approximately  1.2+-0.2.  Corresponding  values  in  wet  rafting  tests  were  approximately  0.2  and 
0.4.  In  pressure  ridge  simulations  [Hopkins],  the  various  friction  coefficients  were  found  to 
control  the  relative  volumes  of  the  ridge  sail  and  keel.  The  coefficients  which  gave  the  most 
realistic  balance  between  sail  and  keel  volume  were  in  the  range  0.8  to  1.0  -  very  similar  to  the 
values  found  by  Coon  and  Pritchard.  The  ridge  simulations  were  less  sensitive  to  the  underwater 
(wet)  friction  coefficient,  but  again  values  of  0.3  to  0.6  used  in  the  simulations  were  similar  to  the 
values  reported  by  Coon  and  Pritchard.  Because  sliding  friction  is  the  dominant  energy  sink  in  the 
processes  of  ridging,  rafting,  and  the  compression  of  floating  rubble,  this  data  set  should  be 
studied  carefully. 

Based  on  the  results  of  long-term  ridging  simulations  [Hopkins],  the  pressure  ridging 
process  can  be  divided  into  three  phases.  The  first  phase  begins  with  an  intact  sheet  of  relatively 
thin  ice  impacting  a  floe  and  ends  when  the  maximum  sail  volume  is  reached.  In  the  second 
phase,  the  ridge  keel  maintains  a  triangular  shape,  deepening  and  growing  in  the  leadward 
direction,  eventually  creating  a  rubble  field  of  more  or  less  uniform  thickness.  The  second  phase 
continues  as  long  as  there  is  a  sufficient  driving  force  and  a  supply  of  thin  ice.  The  third  phase  is 
the  compression  of  the  rubble  field  between  floes.  In  this  scenario,  the  formation  of  rubble  fields 
naturally  follows  the  formation  of  ridges. 

The  effect  of  consolidation  on  the  strength  of  first-year  pressure  ridges  and  rubble  fields 
was  the  subject  of  numerous  discussions  during  the  workshop.  In  coastal  areas,  consolidated 
rubble  fields  and  pressure  ridges  are  an  important  consideration  in  the  design  of  offshore 
structures.  Consolidated  pressure  ridges  are  a  major  impediment  to  ship  passage  and  may 
increase  the  severity  and  depth  of  scour  events  in  near-shore  areas.  Future  pressure  ridge 
modeling  efforts  should  focus  on  the  effects  of  consolidation  on  the  strength  of  ridges  and  rubble 
fields.  The  study  of  the  consolidation  of  a  first-year  ridge  by  Shapiro  et  al.  will  be  useful  in  such 
an  effort. 

Opening 

Macro-cracks  in  first-year  and  multiyear  ice  precede  lead  formation.  Deformation  and 
fracture  processes  at  the  scale  of  1  km  are  characterized  by  formation  of  ordered  systems  of 
cracks  as  well  as  coupled  systems  of  cracks  and  ridges.  Tensile  crack  formation  perpendicular  to 
a  ridge  was  detected  acoustically  during  the  SIMI  field  program  [MIT].  Systems  of  cracks 
perpendicular  to  ridge  lines  have  been  theoretically  predicted  by  Goldstein  and  Osipenko  (1995). 

Acoustics  provide  the  only  viable  remote  sensing  capability  for  detecting,  characterizing, 
and  identifying  failure  mechanisms  in  natural  sea  ice  on  the  1-m  to  10-km  scale.  A  major  SIMI 
accomplishment  is  the  development  of  a  layered  approach  to  seismo-acoustic  mapping  of  sea  ice 
fracture  processes.  The  concept  uses  a  32-element  horizontal  hydrophone  surveillance  array  with 
a  500  aperture.  Focus  beamforming  is  performed  using  parallel  processing,  producing  a  real-time 
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map  of  ice  seismicity  over  a  floe  cluster,  permitting  early  detection  of  active  floe  boundaries. 
When  an  active  region  is  identified,  a  mobile  array  of  five  3-component  geophones  and  a  srngle 
hydrophone  is  rapidly  deployed  to  provide  near-field  measurements  of  the  acoustic  emission  from 
the  ice  fractures.  An  important  component  of  this  concept  is  the  development  o  a  nove  lgi 
telemetry  system,  providing  a  wireless  extension  to  a  local-area  network  in  the  camp.  This 
permits  wireless  data  transfer  from  geophone  clusters. 

During  the  Fall  and  Spring  SIMI  experiment  the  layered  concept  successfully  detected, 
located,  and  monitored  several  ice-mechanical  events  [MU].  A  particularly  interesting  event  was 
a  large,  1-km-scale  fracture,  mentioned  above,  which  propagated  transversely  from  a  new  ridge 
located  N  of  East  Camp.  Another  event  was  the  formation  of  a  major  lead  2  km  NE  o  ast 
Camp  in  the  latter  part  of  the  experiment  The  dynamics  and  refreezing  of  this  lead  was 
monitored  by  two  geophone  clusters  for  several  days,  providing  near-field  measurements  of 
acoustic  emissions  from  rafting  and  early  ridge  building  processes.  At  the  same  time  strong 
signals  corresponding  to  what  are  believed  to  be  "edge  waves"  were  recorded.  A  preliminary 
analysis  of  this  data  suggests  that  the  signals  were  associated  with  a  20  s  suck/slip  behavior  along 
a  lead  which  was  part  of  a  large,  30  -  40-km  long,  lead  structure.  Similar  waves  were  recorded 
by  tilt-meters  and  accelerometers  at  a  ridge  on  the  perimeter  of  the  West  Camp  floe  [Wadhams]. 

Another  accomplishment  on  the  floe  scale  was  the  execution  and  monitoring  of  plate  wave 
experiments  performed  at  Resolute  as  part  of  SIMI  using  automated  hammer  blows  and  geophone 
arrays  The  analysis  of  the  data  collected  [MU]  suggests  strongly  frequency  dependent, 
transverse  anisotropy  of  the  sea  ice.  The  precise  nature  of  the  frequency  dependence  is 

unresolved. 


The  analysis  of  SIMI  seismographic  data  will  address  the  following  issues: 

•  Spatial  and  temporal  statistics  of  ice  events  in  the  vicinity  of  the  West  Camp. 
The  spatial  distribution  of  fractures  in  relation  to  old  and  new  ridges  is  of 
particular  interest; 

•  Analysis  of  data  from  geophone  arrays  for  fracture  type,  size,  and  speed.  The 
analysis  is  expected  to  reveal  the  dominant  fracture  mechanisms  associated 
with  rafting  and  ridging.  The  results  should  be  compared  with  recent 
theoretical  predictions  of  compression  induced  ice  fractures  [Goldstein]; 

•  High-resolution  analysis  of  plate  wave  data  has  revealed  frequency  dependent 
anisotropy  of  sea  ice.  This  behavior  may  be  associated  with  unknown 
phenomena  having  a  length- scale  of  20  -  40  m; 

•  Edge  waves  were  observed  in  SIMI  along  leads  [MIT/WHOI]  and  along  ridges 
on  the  floe  perimeter  [Wadhams].  The  floe  recordings  of  Wadhams  indicate 
narrow-line  spectrum,  whereas  the  lead  recordings  by  MIT/WHOI  are  of  a 
more  wide  band  nature  which  suggest  a  20  s  period,  stick/slip  dynamic 
behavior  of  motion  parallel  to  leads.  This  observation  may  shed  light  on  the 
dynamics  of  in-plane  sliding  between  ice  sheets. 
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ABSTRACT 

For  the  first  time,  simultaneous  measurements  of  sea  ice  stress  on  both  sides  of  a 
lead,  ice  motion  from  sequential  SAR  imagery,  and  ice  motion  from  drifting  GPS/Argos 
buoys  are  compared.  The  lead,  which  appeared  early  in  Day  35  (Feb.  4),  1994,  was 
observed,  relative  motions  were  calculated  across  and  along  the  lead  direction,  and 
stresses  were  transformed  to  this  same  coordinate  system.  The  lead  orientation  and 
approximate  location  were  determined  from  SAR  images  of  the  pack  ice  surrounding  the 
SIMI  camp,  as  well  as  net  ice  motion  at  3-day  intervals.  The  ice  motion  around  the  lead 
was  provided  in  greater  temporal  detail  by  an  array  of  drifting  GPS/Argos  buoys 
producing  hourly  positions.  Sea  ice  stress  was  measured  on  each  side  of  the  lead  using 
flatjack  sensors.  The  ice  motion  normal  to  the  lead  and  the  ice  stress  normal  to  the  lead 
are  compared  as  well  as  the  shear  ice  motion  and  the  ice  shear  stress. 

Shortly  before  the  lead  opening  motion  is  apparent  in  the  ice  motion  data,  the  ice 
stress  state  abruptly  changed  from  biaxial  compression  to  uniaxial  compression  with 
traction  on  the  lead  near  zero.  We  thus  measured  the  stress  conditions  at  lead  creation  as 
sheared  compression;  there  was  no  tension  at  the  time.  We  observe  excellent  correlation 
between  ice  motions  measured  with  SAR  products  and  drifting  buoys.  The  SAR  products 
provide  greater  spatial  detail  and  a  more  precise  definition  of  the  orientation  of  the  lead 
while  the  drifting  buoys  provide  greater  temporal  detail.  We  also  conclude  that  a  properly 
installed  20-cm  stress  gauge  in  a  single  floe  can  provide  a  reasonable  measurement  of  the 
large-scale  (50-km)  pack  ice  stress  resultant. 

We  thank  R.  Colony,  R.  E.  Moritz,  and  L.  McNutt  for  their  helpful  technical 
discussions  concerning  these  data. 
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These  results  closely  follow  the  experiment  plan  described  in  Section  IV  of  "Contemporaneous  Field 
Measurements  of  Pack  Ice  Stress  and  Ice  Strain  Measurements  from  SAR  Imagery,"  by  Coon,  Knoke,  Echert,  and  Stern 
(in  Proceedings  of  OCEANS'93,  Victoria,  B.C.,  October,  1993).  The  schematic  below,  first  published  in  1991,  accurately 
describes  the  overview  of  this  1994  experiment.  Note  that  three  satellite  systems  (ERS-1,  GPS,  and  Argos)  were  utilized 
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Ice  Stress  Measurement  Sites  and  In-Situ  Stress  Measurements 

Buoy-  and  SAR-Derived  Ice  Vectors 
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Drawing  by  J.  Overland 


Synthetic  Aperture  Radar  (SAR)  image  with  SIMI  Camp 
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This  figure  shows  the  positions  of  the  PMEL  GPS/Argos  buoys  at  the  SIMI  camp  and  on  the  floes  around  it  at 
0000Z  on  Day  35  (February  4),  1994.  We  also  include  the  orientation  (though  not  necessarily  the  accurate  position)  of  a 
lead  which  separated  Floes  4,  6, 10,  and  12  from  the  camp. 
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Positions  on  Day  35,  1994 


These  figures  depict  3-day  displacements  of  the  GPS/Argos  buoy  array  relative  to  the  displacement  of  the  SIMI 
camp  (plus  sign).  The  left  figure  is  Day  35  to  Day  38;  the  right  figure  is  Day  38  to  Day  41;  the  figure  axes  and  the  buoy 
positions  are  latitude  and  longitude.  The  velocity  vectors,  scaled  by  the  0.2  km  arrow  in  the  lower  left  corner,  indicate 
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9403500-9403800  Relative  Displacement  9403800-9404100  Relative  Displacement 


Vector  Ice  Motion  Comparison 
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JS  positions  and  ]ump  in  shear 
SAR  images 


Day  of  year  1994 
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Day  of  Year  1 994 


The  ice  stress  normal  to  the  lead  and  the  jump  in  ice  velocity  normal  to  the  lead  are  plotted  versus  time.  Positive 
motion  is  opening,  and  positive  stress  is  tension.  Ice  stress  normal  to  the  lead  plummets  shortly  before  the  lead  opening 
is  apparent  in  the  ice  motion  data.  The  stress  normal  to  the  lead  remained  small  (less  than  20  kN/m)  for  the  next  1 1 
days. 
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ABSTRACT 

During  the  April  1994  Sea  Ice  Mechanics  Main  Field  Experiment,  we  conducted  a  set  of  controlled-load  tests 
to  mimic  the  small-scale  processes  that  are  active  during  the  formation  of  rafts  and  ridges.  Each  test  approximated 
the  behavior  of  in  situ  sea  ice,  except  that  each  process  was  isolated  and  the  loading  was  artificially  applied.  In  this 
paper  we  report  results  of  the  tests  in  which  blocks  of  ice  were  dropped  onto  the  ice  sheet. 

The  acoustic  signals  emitted  were  separated  into  the  forced  acoustic  wave  generated  by  the  first  impact  and  the 
flexural  oscillation  of  the  floating  ice  sheet  For  the  forced  impact,  the  total  radiated  acoustic  energy  is  proportional 
to  the  kinetic  energy  at  impact  When  the  ice  sheet  is  bare  (no  snow  cover),  the  ratio  is  about  11x10  .  With  a  1 .3- 
1.5  cm  snow  cover,  the  ratio  is  about  2x10**.  Variations  in  block  size  and  ice  sheet  thickness  are  less  important 
than  the  impact  energy.  The  spectra  were  also  determined  for  all  tests.  The  snow  cover  reduced  sound  amplitude 
above  200  Hz  by  about  10  dB. 

INTRODUCTION 

As  a  ridge  builds,  there  are  numerous  small-scale  processes  that  occur  the  ice  sheet  is  bent,  buckled,  crushed, 
and  cracked,  blocks  pile  up  and  down,  these  same  blocks  tumble  down  and  up  along  the  sail  and  keel,  blocks 
splash  into  the  water,  blocks  are  pushed  over  and  under  the  ice,  and  the  ice  rafts.  This  last  process  is  separated 
from  over  and  under-ride  because  the  ice  blocks  (and  ice  sheet)  are  pushed  up  and  downhill,  and  the  amount  of 
wetness  differs  between  over-ride  and  under-ride.  It  is  difficult  to  quantify  the  physical  processes  that  contribute 
to  the  noise  generated  during  formation  of  a  ridge  because  these  processes  occur  simultaneously  and  in  unknown 
amounts. 

Of  these  processes,  cracking  has  received  extensive  study  (e.g.  Stein,  1 993).  Our  focus  is  on  the  other  processes. 
Furthermore,  when  we  listen  to  under  ice  noise,  it  is  not  cracking  that  dominates,  but  instead  it  is  blocks  falling 
against  the  ice  sheet,  and  blocks  sliding  over  and/or  under  the  ice  sheet,  and  rafting  of  ice  sheets. 


As  one  component  of  the  April  1994  Sea  Ice  Mechanics  Main  Field  Experiment  (e.g.,Coon  et  al  1995  Yie  and 
Farmer,  M  FnM,  1995),  wc  conduced  ,  *  of  c^mono.  E«h 

f»‘(T  o  ’"  'X"P'  *"  «*  P™““  «■»  isoldlcd  cud  the  loading  wee  amiict.ll. 

applied.  The  controlled- load  experiments  included  over-ride,  under-nde,  rafting,  flopping,  dropping  underwater 

bumping  and  splashing.  Each  of  these  processes  occurs  during  the  formation  of  a  ridge  or  raft.  Each  type  of  test 
is  now  described  more  completely:  over-ride  -  towed  a  cold  ice  block  over  the  surface  of  the  intactiw  sheer 
mder-nde  -  towed  a  warm  ice  block  under  the  intact  ice  sheer,  raft  -  towed  a  warm  ice  block  (freshly  excavated 
nf?ZeH  UP  Md  °VCr  a  cant,levercd  beam  ^on  of  the  ice  sheer,  drop  -  dropped  a  block  of  ice 

£2  fa  i  at  3  ,  ,egr  a"Sle)  0nt0  the  intact  icc  sheer’  -  tipped  a  standing  block  offce  over  onto  the 

TficI  m  not  SPAW  PPe  aStand,ng  bl0Ck  0fice  °ver  int0  the  water-  bump  -  allowed  a  submerged  block 
of  ice  to  float  up  and  bump  against  the  underside  of  the  ice  sheer,  cantilever  -  lifted  the  free  end  of  a  floating 

cantilever  bram  until  the  beam  failed  in  flexure;  and  fracture  -  split  a  floating  ice  sheet  with  air  jacks  8 

J°™  °£’hcse  pr0CCSSe;  afTeCt'he  mechanical  response  of  the  ice,  whereas  all  processes  generate  acoustic 

iuL  iZ'”  3VeS  r  m  'he  imf>aCt  iS  3  VaiU3blc  Si8nal  f0r  Prob'n8  the  interaction  process  even 
,?°  ffi  ?  £  CTergy  accounts  for  an  extremely  small  amount  of  the  input  mechanical  energy.  The  signal  is 

ZlTn  Tn‘S  “ ,S  readi,y  mCaSUrable’ and  ite  si8natures  reflect  intetactio^pr^ss  The 

°W  P°SSlbIe  reconstruction  ofsuch  Processes.  Hydrophones,  geophones,  ice  strain 
sensors,  and  ice  tilt  sensors  were  installed  in  arrays  around  the  controlled-load  tests 

d^rVlt  WOfk,iS  t0  Underfn£  thC  Iar8e'SCale  behavior  ofs*a  ice-  ar|d  the  noise  emitted  as  it  moves  and 
deforms.  The  large-scale  process  that  dominates  ice  behavior  on  this  scale  are  rafting  and  ridging  Therefore  we 

mechanical  ^  Ra?lng  13  also  a  P^ursor  to  ridging.  With  better  understanding  of  the  small-scale 

hone  to  ,tha' C®ntntiute  t°  ndgmg  and  of  the  acoustic  signals  that  are  emitted  by  these  processes,  we 

hope  to  develop  better  models  of  ambient  noise  under  the  ice  and  inversely  to  use  sound  measurements  as  a  gage 

to  rr  C  “  a  llbrary  0f  pr0babili(y  distnbut,on  signatures  in  phase  space  corresponding 

benartmn  T'T*  "T1™'?5'  A  8eneral  Arct'c  noise  time  series  of  an  ice  deformation  event  codd  then 
be  partitioned  into  the ;  library  of  events.  The  converse  may  also  be  possible,  that  is,  to  generate  artificial  or 
predictive  time  senes  of  ambient  noise  that  include  the  appropriate  proportions  of  library  signatures 
Here  we  present  the  acoustic  signatures  of  these  drop  and  flop  tests.  We  estimate  the  acoustic  emissions  from 

ene  gy  Mnv  bl  ^  °f ‘0tal  radiated  aC0UStic  enerSy t0  the  mechanical  input 

energy.  The  falling  blocks  studied  here  result  from  rearrangement  of  blocks  in  the  ridge  sail  during  the  ridge 

?eD0se8TPh^'  red,Stnb#U‘‘on^as  first  modeled  by  Parmerter  and  Coon  (1972)  to  match  the  angle  of 

mode^  rt  ,  ?!ener8y  °f  he  ndge  W3S  then  incorP°rated  explicitly  by  Rothrock(\91 5)  into  the  Aidjex 
model.  This  and  other  energy  measures  were  then  used  by  Pritchard  (1984,  1990)  as  a  proxy  variable  for 

S™8 amb,en!n0'S* 87erated  ^  ndSmg-  Finai'y.  Hopkins  ( 1994)  estimated  the  energy  baLS  dunng  ridge 
building  using  a  discrete  element  model  of  the  forces  acting  on  individual  ice  blocks  in  the  ridge. 

ACOUSTIC  DATA  ACQUISITION  AND  ANALYSIS 

deolovedZborh  ZZZTZ™"1  smal|-scale  ice  Processes,  an  array  of  geophones  and  hydrophones  were 
t£|  a  p^r  E  a"d.  test  3,tes  from  the  Main  camp.  Figure  1  is  a  sketch  of  aerial  view  of  the  arrav 
omSit  fH3'  1  a7y  TS,StCd  0f  five  tn'axial  seophones  (frozen  into  the  ice  surface)  and  the  five 
ZmZr  Tn  hydr0pti°nes  (deployed  at  5-m  depth  in  the  water).  The  maximum  spacing  of  the  array  was 
We  arc  reportin8  on  twe*ve  tests  at  this  site.  A  similar  array  was  deptoyed  on  the  West 
.  owever,  the  West  site  array  had  a  much  greater  maximum  spacing  of  up  to  200  meters-  the  hydrophones 
were  deployed  at  20-meter  water  depth.  We  are  reporting  on  four  tests  a.  this  she.  Detailed 

aZZTa^  MmXieandFarmer^^  1995),  Coon  et  al.  (1995),  and  Pritchard  ( 1 995)  For  each  test 
acoustic  and  seismic  emissions  were  digitally  recorded  on  VHS  tape. 

unde^Sh" thc  mechanioal  waves  recorded  on  a  sensor  some  distance  away  from  the  impact  point,  one  must 
understanh  how  the  seismic  and  acoustic  waves  propagate  in  the  Arctic  environment.  There  have  teen  extensive 
wave  prWtton  mechanisms  in  the  Arctic  ocean  (see  Stein,  1993 ,  Are  and  Farmer,  1994b). 
There  are  two  propagation  channels  for  mechanical  waves  generated  by  the  ice  process:  the  floating  ice  sheet  and 
acoustically  up-refractive  water  column.  The  ice  sheet  guides  vanous  elastic  waves  in  the  ice  Tte  w“er Sum 
provides  signal  channels  for  acoustic  waves  coupled  into  the  water  from  the  impacted  local  ice-  it  also  allows  the 
e^age  ofe.asbc  waves  from  the  ,ce  sheet  into  the  ocean.  It  has  teen  shown  J the  charactUZ ^ o'the  gu  ded 
elasttc  waves  are  mainly  determined  by  the  mechanical  properties  of  the  ice  and  its  geometry,  especially  the 
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thickness.  The  elastic  waves  relevant  to  the  present  study  are  flexural  waves,  P-waves  and  SH-waves.  The 
acoustic  wave  that  is  coupled  directly  from  the  ice  impact  has  a  signature  that  is  affected  both  by  the  local  ice 
properties  and,  most  importantly,  by  the  interaction  process  between  the  falling  block  and  the  ice  sheet.  It  is  the 
acoustic  wave  that  provides  more  valuable  information  on  the  interaction  process. 

Both  the  SH-wave  and  flexural  waves  are  trapped  in  the  ice  and  upper  water  column  due  to  horizontal 
polarization  and  subsonic  phase  velocities,  respectively.  However,  P-waves  leak  into  the  water  column  in  an 
efficient  manner  governed  by  Snell’s  law.  Acoustic  waves  can  propagate  for  a  long  distance  (of  order  10  km  in 
a  500  meter  deep  water  column)  before  being  reflected  from  the  bottom  of  the  ice  sheet.  A  layer  of  sea  ice 
enriches  the  physics  ot  underwater  sound  field  in  the  Arctic  ocean  by  projecting  various  elastic  waves  into  the 
ocean.  The  co-existence  of  elastic  (leaked  from  the  ice)  and  acoustic  (coupled  directly  from  local  failure  processes 
in  the  ice)  waves  complicates  the  underwater  sound  field  especially  in  the  upper  water  column.  From  a  signal 
detection  point  of  view  (in  the  water),  one  observes  different  signal  structures  arriving  at  different  times, 
corresponding  to  the  above  described  wave-process.  With  these  theoretical  descriptions  and  understandings  of 
seismic  and  acoustic  processes,  we  can  examine  the  observed  signals  . 

The  total  radiated  acoustic  power  assuming  a  dipole  distribution  is 

RA([)  =  f—ds  (1) 

J  pc 

where  p(r,B)  =  pocos0/r,  pc  is  acoustic  impedance,  0  is  the  radiation  angle,  p0  is  source  pressure  observed  at 
range  r,  and  s  is  the  spherical  radiation  area  at  range  r.  The  total  radiated  acoustic  energy  is 

TRAE  =  |  RA(f)dt  m  (2) 


which  can  be  rewritten  using  ParsevaFs  equality  as 


TRAE  =  / 


cos20 

r2pc 


ds ' 


0 


(3) 


where  pJJ)  is  the  Fourier  transform  of  p0(t). 

RESULTS 

On  8  April  1994  twelve  drop  tests  were  conducted  at  the  East  Camp  and  on  22  April  1994  four  tests  were 
conducted  at  West  Camp.  In  each  test  a  block  of  ice  was  either  set  on  its  edge  and  tipped  so  that  it  flopped  onto 
its  side  or  it  was  dropped  onto  the  ice  sheet.  On  8  April  the  ice  sheet  was  0.25  m  thick,  it  was  covered  by  1 .3  cm 
of  snow  (although  this  was  scraped  off  for  some  of  the  tests),  temperature  of  the  upper  surface  (under  the  snow) 
was  -16.4  C,  and  the  air  temperature  was  -20  C.  On  22  April  the  ice  sheet  was  0.51  m  thick,  covered  by  1.5  cm 
of  snow,  and  air  temperature  was  -12.3  C. 

A  falling  ice  block  dissipates  its  initial  potential  energy  by  interacting  with  the  surrounding  ice  fields  before 
reaching  the  minimum  energy  level.  One  of  the  common  interacting  modes  is  the  impact  of  these  free-fall  ice 
blocks  onto  the  surface  of  the  underlying  ice  sheet.  Xie  and  Farmer  (1994b)  found  that  a  weight  (an  i  l-kg  lead 
ball)  dropped  onto  a  2-meter  thick  first-year  ice  sheet  created  a  dipolar  acoustic  source  beneath  the  ice.  Very  small 
fractions  of  the  total  potential  energy  were  radiated  out  from  the  impact  as  acoustic  waves  in  the  water  column  and 
flexural  waves  in  the  ice  sheet.  This  implies  that  most  potential  energy  of  a  falling  weight  is  dissipated  by 
deforming  the  local  ice  at  the  impact. 

Figure  2  shows  a  time  series  of  the  acoustic  signals  for  Test  8- 1 :  an  ice  block  flopping  onto  a  snow-covered  ice 
sheet  These  signals  were  recorded  by  the  array  of  five  hydrophones  deployed  at  the  East  site  (see  Figure  1).  Since 
the  impact  occurred  very  close  to  Ht  (about  6  meters  in  range),  the  signal  saturated  this  channel  precluding 
estimation  of  acoustic  source  levels  from  this  phone.  Figure  3  is  the  acoustic  signal  due  to  a  similar  flopping  test 
(Test  8-2)  except  that  snow  was  cleared  off  from  the  ice  surface  before  the  test.  These  signals  are  discussed  in 
the  following  three  paragraphs. 

The  nearest  hydrophone  (Ht  )  was  saturated  for  most  tests.  Channel  H2  was  used  for  most  analyses  because  it 
was  farthest  from  the  impact  site  and  its  signal  had  the  best  separation  of  acoustic  and  flexural  waves.  However, 
to  see  the  flexural  wave,  we  must  look  at  H3  or  another  hydrophone  because  the  large,  heavy  data  acquisition 
boxes  were  located  along  the  path  to  H2,  which  reduced  its  flexural  wave  signal.  We  are  interested  primarily  in 
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“ih^LT t^ZrmTi  by  ,h'.ta  ta““  *  W=*  farther  ta  fa  ofa,  sig„l, 
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oresented  in  the  T?w  r  S’  Th  V3lue  ,S  used  t0  heiP  isoiatc  thc  fl«ural  wave.  The  TRAE  levels  are 
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were  excited  by  this  impact.  reStmS  f°  "°te  ,hat  no  S18n,ficant  fl«ural  waves 

Pr°Pcrties  include: .  block  size  and  shape 

°m  Tht  r  ™ wtre  '■£***” 

81).  The  total  energy  input  dunng  each  test  is  defined  by  the  potential  energy  change  PE  =  PP  PP  ’ 
potermal  energy  at  *e  begtnning  of  the  test  was  PE„  potential  Lr^at  the  endofeache^tCl  p7=  SpTh 
kmehc  energy  at  the  moment  of  first  .mpact  was  KE.  For  the  flop  tests,  initial  potential  encr^  was 


«?,  -  fa 


2+H2)m 


(4) 
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V 


(5) 
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KE  =  W(Z  -  fi(L+Hy4).  (6) 
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Table  I.  Parameters  for  Drop  and  Flop  Tests. 


Test 

Test 

L 

B 

H 

M 

Snow 

initial  CG 

Final  CG 

No. 

Type 

(m) 

(m) 

(m) 

(kg) 

Z(m) 

(m) 

8-1 

Flop 

0.50 

0.50 

0.25 

60 

Yes 

0.28 

0.13 

8-2 

Flop 

0.50 

0.50 

0.25 

60 

No 

0.28 

0.13 

8-3 

Flop 

1.00 

0.50 

0.25 

120 

No 

0.52 

0.13 

8-4 

EH 

0.53 

0.51 

0.10 

24 

No 

0.27 

0.05 

8-5 

0.53 

0.51 

0.24 

58 

Yes 

0.43 

0.12 

8-6 

Drop  Flat 

0.53 

0.51 

0.24 

58 

No 

0.73 

0.12 

8-7 

Drop  45 

0.48 

0.51 

0.24 

56 

No 

0.61 

0.12 

8-8 

Drop  Flat 

0.53 

0.51 

0.25 

62 

No 

0.43 

0.13 

8-9 

Drop  Flat 

0.43 

0.36 

:  0.06 

9 

No 

1.25 

0.03 

8-10 

Drop  45 

0.51 

0.46 

0.10 

24 

No 

1.22 

0.05 

8-11 

0.53 

0.51 

0.08 

24 

No 

0.61 

0.04 

8-12 

Drop  45 

0.99 

0.53 

0.10 

49 

Yes 

0.91 

0.05 

22-4 

599H 

0.51 

0.48 

0.07 

15 

Yes 

0.26 

0.03 

22-5 

sn 

0.49 

0.50 

0.15 

38 

Yes 

0.26 

0.08 

22-6 

Flop 

0.51 

0.48 

0.07 

15 

No 

0.26 

0.03 

22-7 

0.49  |  0.50  |  0.15  | 

38 

No  1  0.26  | 

0.08 

Table  H  Total  Radiated  Acoustic  Energy  ( TRAE ),  Potential  Energy  (PE),  and  Impact  Energy  (KE)  for  Drop  and 
Flop  Tests. 


Test 

Test 

PE 

KE 

TRAE 

TRAE/PE 

TRAE/KE 

Comments 

No. 

Type 

w 

(j> 

(j> 

xlO" 

xl0‘ 

8-1 

Flop 

89 

89 

0.000444 

5.00 

5.00 

AE+Flexure 

8-2 

Flop 

89 

89 

0.001034 

11.64 

11.64 

AE+Flexure 

8-3 

E39H 

459 

459 

0.004053 

8.83 

8.83 

AE+Flexure 

8-4 

559H 

52 

52 

0.001171 

22.50 

22.50 

8-5 

jDrop  Flat 

172 

172 

0.000419 

2.43 

2.43 

Double  Impact  ! 

8-6 

344 

344 

0.002572 

7.47 

7.47 

8-7 

Drop  45 

268 

194 

0.001907 

7.13 

9.83 

3  Imoacts 

8-8 

186 

186 

0.002341 

12.60 

12.60 

AE+Flexure 

8-9 

114 

114 

0.000961 

8.44 

8.44 

AE+Flexure 

8-10 

275 

237 

0.000786 

2.85 

3.32 

AE  +  Multi-bump 

8-11 

S3SE91 

135 

93 

0.000277 

2.06 

2.98 

AE  +  Multi-bump 

8-12 

419 

256 

0.000264 

0.63 

1.03 

Multiple  Impact 

22-4 

Flop 

34 

34 

0.000010 

0.29 

0.29 

AE 

22-5 

5KHH 

67 

67 

0.000055 

0.82 

0.82 

AE 

22-6 

559  II 

34 

34 

0.000536 

15.15 

15.15 

22-7 

Flop  | 

67 

67 

0.000985 

14.70 

14.70 

AE  +  small  bump  f 
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Parameters  listed  in  the  Table  I  were  varied  in  an  attempt  to  understand  the  physical  behavior  of  falling  blocks 
of  ice.  It  was  our  desire  to  understand  which  variables  affected  the  sound  emitted  and  to  quantify  these  effects. 
All  blocks  dropped  onto  bare  ice  broke  (tests  8-2  through  8-4,  and  8-6  through  8-11).  This  was  true  whether 
flopped  or  dropped  flat  or  at  an  angle.  One  result  of  such  breaking  was  to  generate  multiple  impacts  from  all  pieces 
of  the  block,  rather  than  clean  single  or  double  impacts  when  dropped  flat  or  at  an  angle.  Although  the  block  in 
test  8-5  was  initially  level,  it  was  dropped  at  an  (unknown)  angle,  which  caused  a  double  impact.  Test  8-7  was 
dropped  at  an  angle  and  produced  two  distinct  impacts;  first  when  the  edge  hit  the  ice  sheet,  and  second  when 
it  flopped  over  flat. 

We  assume  that  a  relationship  exists  between  the  total  energy  radiated  by  the  acoustic  waves  and  the  energy  input 
by  dropping  the  block,  and  we  test  this  hypothesis.  The  energy  available  from  dropping  the  block  is  easily 
estimated  by  calculating  the  change  in  potential  energy  before  dropping  the  block  and  the  instant  it  impacts  the 
ice  sheet.  This  calculation  is  somewhat  different  for  the  various  geometries. 

The  TRAE  is  compared  directly  with  the  impact  kinetic  energy  (KE)  in  Figure  6.  Open  symbols  indicate  ice 
blocks  falling  onto  snow-covered  ice;  black  symbols  onto  bare  ice.  Different  symbol  shapes  have  been  used  to 
identify  the  different  kinds  of  tests  (flat  drop  -  square,  drop  at  45  degrees  -  diamond,  and  flop  -  triangle).  For  all 
drops  on  a  bare  ice  sheet,  the  average  ratio  of  TRAE  to  potential  energy  was  1 1  xl  O-6.  Two  data  points  (tests  8- 1 0 
and  8-1 1  had  lower  ratios.  Both  of  these  ice  blocks  shattered,  which  likely  cause  the  lower  TRAE  values.  When 
snow  is  present,  the  average  ratio  of  TRAE  to  potential  energy  is  reduced  to  2x1  O’4.  These  results  appear  to  be 
consistent  over  the  full  range  of  impact  energy. 

The  data  were  also  examined  for  correlations  with  block  area,  mass,  test  type,  and  angle  of  the  drop.  None  were 
as  strong  as  variations  with  snow  cover.  We  recognize  of  course  that  this  is  a  small  test  set,  and  definitive 
correlations  are  difficult.  For  the  lead  ball  used  by  Xie  and  Farmer  (1994a),  the  average  ratio  of  TRAE  to 
potential  energy  was  about  20x1 0"6,  so  that  ice  blocks  are  less  than  half  as  efficient  at  generating  acoustic  energy 
as  the  idealized  spherical  lead  ball.  We  noted  earlier  that  these  block  sizes  are  representative  of  ice  blocks  found 
in  ridges.  We  note  here  that  the  range  of  distances  from  which  they  were  dropped  (0.15-1.22  m)  is  also 
representative  of  distances  that  blocks  fall  during  ridging.  We  have  also  compared  the  total  change  in  potential 
energy,  but  the  ratio  between  TRAE  and  PE  is  not  are  consistent  as  is  the  ratio  of  TRAE  to  KE. 

CONCLUSION 

Sixteen  tests  were  conducted  in  which  ice  blocks  were  either  dropped  or  flopped  onto  the  ice  sheet.  The  acoustic 
signals  emitted  were  separated  into  the  forced  acoustic  wave  generated  by  the  first  impact  and  the  seismic  wave 
generated  as  a  flexural  oscillation  of  the  floating  ice  sheet.  We  focused  analysis  on  the  acoustic  wave  that  has 
approximately  a  dipole  pattern.  It  propagates  to  longer  ranges  than  does  the  flexural  wave. 

When  the  ice  blocks  were  dropped  at  an  angle,  several  distinct  impacts  observed.  When  the  blocks  were  dropped 
onto  bare  ice  they  broke,  which  also  generated  multiple  impacts. 

The  total  radiated  acoustic  energy  was  found  to  be  proportional  to  the  kinetic  energy  at  impact.  When  the  ice 
sheet  was  bare  (no  snow  cover),  the  ratio  was  about  1  lxl  O'4.  With  a  1.3-1 .5  cm  snow  cover,  the  ratio  was  about 
2x1  O'4.  Variations  in  block  size  and  shape  appeared  to  be  less  important  than  the  impact  energy.  There  was  no 
apparent  dependence  on  ice  sheet  thickness  (over  the  range  0.25-0.51  m)  for  the  first  impact  acoustic  signal 
studies  here.  We  would  expect  such  a  dependence  to  exist  for  the  flexural  wave.  The  spectra  were  also 
determined  for  all  tests.  The  snow  cover  reduced  sound  amplitude  above  200  Hz  by  about  10  dB. 

These  results  suggest  that  a  model  could  be  developed  to  describe  the  acoustic  emission  from  the  blocks  falling 
during  formation  of  a  ridge.  The  model  would  require  that  we  estimate  the  kinetic  energy  available  for  blocks 
falling  onto  the  ice  sheet  and  the  probability  distribution  of  events,  but  these  variables  can  be  estimated  from  in 
modem  ice  dynamics  models. 
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Figure  1.  Hydrophone  Array  at  SIMI  East  Camp. 
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Figure  4.  Sound  Emitted  by  Ice  Block  Dropping  Eat  onto  Bare  Ice  Sheet  (Test  8-9). 
Multiple  Impacts  Result  from  Shattering  of  the  Ice  Block.  Top:  Pressure  History  at  Five 
Hydrophones.  Bottom:  Sound  Spectrum  Level  (SSL)  at  Hydrophone  H2. 
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Figure  5.  Sound  Emitted  by  Ice  Block  Dropping  at  45  Degree  Angle  onto  Bare  Ice  Sheet 
(Test  8-10).  Multiple  Impacts  Result  from  Initial  Angle  of  Drop  and  Shattering  of  the  Ice 
Block.  Top:  Pressure  History  at  Five  Hydrophones.  Bottom:  Sound  Spectrum  Level 
(SSL)  at  Hydrophone  H2. 
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ABSTRACT 

The  strengths  of  the  consolidated  and  partially  consolidated  parts  of  a  first-year  sea  ice  ridge  keel 
may  be  the  design  parameters  for  structures  to  be  placed  in  regions  covered  seasonally  by  first-year  ice. 
These  strengths  are  also  important  for  large-scale  ice  dynamics  models  which  can  account  for  leads  and 
discontinuous  ice  motion.  The  strength  of  a  first-year,  partially  consolidated  keel  has  not  been  meas¬ 
ured  with  certainty  because  of  the  heterogeneous  nature  of  the  broken  ice  and  the  large  size  of  the  keel. 
The  authors  have  developed  and  have  performed  an  initial  field  test  of  a  system  to  measure  the 
strength  of  a  first-year  partially  consolidated  keel.  The  system  consists  of  a  wheeled  chain  saw,  a 
60x60x2 -cm  air  jack  (a  stiffened  rubber  bladder),  an  air  compressor,  and  an  air  pressure  recorder. 
After  removing  a  portion  of  the  sail  of  the  test  ridge,  the  wheeled  chain  saw  is  used  to  cut  a  vertical 
slot  into  the  keel  for  the  air  jack.  Inflating  the  jack  provides  the  pressure-displacement  curve  which  in 
turn  provides  the  limit  load  of  plastic  deformations  of  the  keel,  an  index  of  keel  strength.  Four 
measurements  were  made  on  a  small,  first-year  ridge  in  the  Beaufort  Sea  during  April  of  1994.  The 
maximum  pressures  were  variable:  20,  40,  210  and  greater  than  700  kPa.  Calculations  of  the  load  for 
a  continuing  motion  solution  based  on  laboratory  test  data  on  unconsolidated  rubble  resulted  in  a 
significantly  smaller  value  than  the  measured  limit  load.  We  believe  significant  consolidation  of  the 
keel  occurred  in  the  two  weeks  between  the  ridging  event  and  our  tests.  The  fully  consolidated  layer 
between  the  sail  and  keel  generally  provides  the  majority  of  the  ridge  strength.  Future  work  should 
include  tests  of  a  ridge  in  fall,  winter,  and  spring  to  track  strength  with  consolidation. 


INTRODUCTION 

First-year  ridges  can  control  the  total  load  on  structures  in  parts  of  the  Arctic  (Wang  and  Poplin, 
1995).  First-year  ridges  are  also  involved  in  the  dynamics  of  pack  ice  through  the  ice  thickness  distri¬ 
bution  (Thorndike  et  al.,  1975)  and,  in  models  such  as  the  oriented  pack  ice  constitutive  law,  the  large- 
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scale  strength  (Coon  et  al.,  1992).  The  strength  of  ridges  is  as  important  as  the  ridge  shape,  frequency 
of  occurrence,  and  keel  depth  because  the  strength  is  required  to  determine  the  load  on  a  structure  or 
ship  as  well  as  the  amount  of  force  that  pack  ice  can  transmit.  Field  measurements  of  ridge  strength 
provide  scale  for  laboratory  measurements  done  on  the  friction  and  cohesion  of  broken  ice  in  a  shear 
box  (Ettema  and  Schaefer,  1986;  Ettema  and  Urroz,  1989;  Keinonen  and  Nyman,  1978;  Urroz  and 
Ettema,  1987;  Weiss  et  al.,  1981;  and  Wong  et  al.,  1987).  For  example,  Prodanovic  (1979)  ran 
laboratory  model  tests  to  measure  the  shear  strength  and  the  simulated  loads  on  cylindrical  structures. 
Allyn  and  Wasilewski  (1979)  analyzed  the  strength  of  nibble  fields  in  connection  with  artificial 
islands.  The  laboratory  tests,  however,  are  nearly  always  performed  on  unconsolidated  ice  blocks.  As 
a  ridge  gradually  changes  from  being  unconsolidated  through  partially  consolidated  (as  ours  was)  to  a 
fully  consolidated  multi-year  ridge,  its  strength  increases. 

This  paper  presents  a  description  of  the  ridge  tested  (block  size  distribution,  porosity,  temperature, 
and  salinity)  followed  by  the  test  procedure,  the  test  results,  hindcast  analysis,  and  discussion. 


EXPERIMENT 

The  force- displacement  measurements  of  a  first-year  pressure  ridge  keel  were  made  about  500  km 
north  of  Prudhoe  Bay,  Alaska  on  April  14-16,  1994,  with  multi-year  ice  on  both  sides  of  the  ridge  and 
rc frozen  lead  studied.  The  equipment  for  this  experiment  was  tested  in  the  fall  of  1993,  and  data  were 
collected  in  the  spring  of  1994  with  modified  equipment  (see  Coon  et  al.,  1994).  The  strength  of  first- 
year  ridged  ice  was  measured  using  an  air  jack  (a  stiffened  rubber  bladder)  to  simultaneously  load  the 
ice  and  measure  the  displacement.  In  addition,  we  characterized  the  sail  and  keel  geometry  of  a  ridge 
formed  from  first-year  ice,  and  inflated  the  air  jack  in  four  slots.  The  ridge  configuration  in  the  vicin¬ 
ity  of  the  four  air-jack  slots  was  characterized  by  drilling  12  holes  and  with  observations  with  an 
underwater  camera. 

The  ridge  was  selected  by  searching  the  area  around  the  SEMI  East  Camp  from  a  helicopter,  look¬ 
ing  for  a  crisp,  new,  first-year  ridge.  The  ridge  was  about  2  km  from  the  SIMI  camp,  which  was 
manned  during  the  formation  of  the  lead  from  which  the  ridge  was  built  and  during  the  ridge  building. 


Rldae  data 

Holes  were  drilled  to  establish  the  extent  of  the  keel  and  refrozen  lead  ice.  The  lead  ice  had 
ridged  against  a  multi-year  floe,  causing  a  large  block  of  multi-year  ice  to  fracture  and  tilt  down.  The 
refrozen  lead  ice  had  continued  to  grow  thicker  after  the  ridging  event,  which  probably  took  place  on 
March  30.  After  looking  at  the  underside  of  the  test  area  with  an  underwater  video  camera  using  poles 
inserted  through  other  holes  in  the  ice  for  visual  reference,  we  cleared  off  some  sail  blocks  and 
selected  the  locations  for  cutting  slots  for  the  air  jack. 


Block  size  distribution.  Nearly  all  the  blocks  in  a  100  block  survey  of  the  ridge  sail  had  a 
thickness  of  0.14  m.  Their  longest  dimension  averaged  0.49  m  with  a  standard  deviation  of  0.29  m. 
Based  on  their  average  width,  the  block  area  averaged  0.20  m2  with  a  standard  deviation  of  0.22  m2 
These  results  are  typical  for  ridges  built  from  first-year  ice,  falling  between  the  two  empirical  fits 
provided  in  Tucker  and  Govoni  (1981)  on  their  Figure  10. 


Porosity.  To  establish  the  location  of  voids  in  the  keel  near  the  air  jack  slots,  three  sets  of  holes 
were  drilled  and  logged.  Figure  1  shows  the  plan  view  of  the  slots  and  the  holes;  there  are  two  rows  of 
holes  outside  the  test  area  and  a  set  of  holes  drilled  alongside  the  slots  after  the  force-deformation 
tests.  Figure  2  illustrates  the  locations  of  the  voids  in  the  two  rows  of  holes.  The  sail  is  to  the  left  in 
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these  figures  while  the  0.46  m  refrozen  lead  ice  is  to  the  right.  Each  hole  drilled  started  with  a  hard, 
solid  layer  which  averaged  0.56  m  thick;  at  the  bottom  of  that  layer,  we  hit  softer  ice  or  a  void.  After 
the  force-deformation  tests,  the  holes  shown  in  Figure  3  were  drilled  and  logged.  Hole  numbers  12 
and  13  were  drilled  alongside  the  slots  for  Tests  1  and  2;  holes  14  and  15  alongside  Test  3;  and  holes 
16  and  17  alongside  Test  4.  The  average  porosity  for  all  the  holes  was  25  to  35%,  depending  upon 
whether  mushy  soft  ice  is  considered  as  solid  ice  or  void.  These  results  are  similar  to  those  of 
Frederking  and  Wright  (1982). 


Sail  of  first-year  ridge 
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FIGURE  1.  PLAN  VIEW  OF  TEST  SITE 

The  measured  ice  thicknesses  in  the  refrozen  lead  and  the  ridge  blocks  were  compared  to  values 
computed  from  ice  camp  temperatures  using  the  formula  for  mean  ice  growth  developed  by  Zubov 
(1944), 

I2+50I-8R  =  0  (1) 

in  which  I  is  the  ice  thickness  in  cm  and  R  is  the  cumulative  freezing  degree  days  in  °C-days  using  0°C 
as  a  reference.  As  indicated  in  Figure  4,  the  lead  was  most  likely  formed  during  an  ice  motion  event 
observed  at  the  East  Camp  on  March  26  and  the  ridge  was  probably  formed  during  the  storm  and  ice 
motion  event  observed  at  East  Camp  on  March  30.  The  resulting  ice  thicknesses  for  the  ridge  blocks 
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are  consistent  with  Zubov's  formula.  The  formula  under-predicted  the  lead  ice  thickness  on  April  14 
by  17%.  r 
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FIGURE  2.  TWO  KEEL  CROSS-SECTIONS  SHOWING  VOIDS  AS  WHITE 


2 


Hole  Number 


FIGURE  3.  RELATIVE  LOCATIONS  OF  THE  AIR  JACK  TEST  LOCATIONS  AND  THE 
SIX  HOLES  DRILLED  ALONGSIDE  THE  AIR  JACK  SLOTS  (VOIDS  ARE  SHOWN  AS 

WHITE  IN  THE  HOLES) 
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Tempera  ture/saiinity,  Temperatures  of  blocks  from  the  keel  were  measured  with  a  thermistor; 
the  results  are  provided  in  Table  1.  The  temperature  measurements  were  obtained  by  drilling  small 
holes  into  the  ice  blocks  pulled  from  the  first  air  jack  slot,  which  extended  down  into  the  keel. 


TABLE  1.  TEMPERATURE  AND  SALINITY  DATA 


Location  of  reading 

Temperature  (°C) 

Salinity  (%o) 

Sample  from  15  cm  below  surface 

-9.4 

5.4 

Sample  from  40  cm  below  surface 

^,1 

4.4 

_ Sample  from  75  cm  below  surface 

-1.4 

5.0 

Salinity  samples  were  taken  from  the  same  blocks;  these  results  are  also  given  in  Table  1.  A  field 
study  by  Frederking  and  Wright  (1982)  indicate  that  4-5  %c  salinity  is  typical  for  submerged  rubble. 


Load  tests 

The  load  tests  in  the  keel  were  performed  with  an  air  compressor  inflating  an  air  jack  in  a  slot  in 
the  keel  while  recording  the  air  pressure  with  a  pressure  transducer  and  a  datalogger  (see  Figure  5). 
The  air  jack  allowed  a  compact,  light-weight  system  to  deform  the  ice  up  to  30  cm  (15  cm  in  each 
direction  from  the  centerline)  and  measure  the  deformation.  The  test  is  analogous  to  the  ice  strength 
measurement  test  procedure  described  by  Iyer  and  Masterson  (1987)  who  froze  thin-walled  metal  flat- 
jacks  into  multi-year  ice  to  determine  its  strength. 

<0 


(c)  Cross  Section  Parallel  to  the  Axis  of  the  Ridge 


FIGURE  5.  RIDGE  STRENGTH  EXPERIMENT 
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The  jack  displacement  was  derived  from  the  compressor  performance  curve  and  the  recorded  pres¬ 
sure  history  by  a  series  of  laboratory  calibration  tests  on  the  air  compressor  and  the  air  jack.  The  slot 
for  the  air  jack  was  cut  with  a  modified  wheel-mounted  chain  saw  with  a  1.5-m  bar,  as  shown  in 
Figure  6,  made  by  the  Ditcher  Saw  Company,  Inc.  Two  cuts  were  made  for  each  slot;  the  blocks  of  ice 
between  were  lifted  from  the  slot  for  temperature  and  salinity  measurements.  From  probing  the  slots, 
we  know  the  slot  cutting  did  not  dislodge  or  free  up  the  keel  blocks. 


FIGURE  6.  WHEEL-MOUNTED  CHAIN  SAW  FOR  CUTTING  AIR  JACK  SLOTS 


The  Power  Team  air  jacks  used  were  60x60x2  cm  thick  with  a  maximum  rated  pressure  of  800  kPa 
and  a  displacement  of  33  cm.  We  had  no  trouble  dropping  the  air  jack  into  these  precut  slots  in  the 
manner  shown  in  Figure  7.  Once  in  the  slot,  the  air  jack  was  inflated  by  a  small  Black  and  Decker  air 
compressor  (for  automobile  tires)  which  operated  off  a  12V  lead-acid  gel-cell  battery.  The  air  pressure 
was  measured  with  a  factory-calibrated  PSI-tronix  pressure  transducer  and  recorded  with  a  Campbell 
Scientific  CR-10  data  logger  and  SM-192  storage  module,  which  operated  off  a  second  12V  gel-cell 
battery. 

We  performed  four  inflation  tests.  As  illustrated  in  Figure  3,  the  first  test  was  in  the  consolidated 
layer,  the  second  test  was  in  the  partially  consolidated  ice  directly  below  it  in  the  same  slot,  and  the 
third  and  fourth  tests  were  in  the  partially  consolidated  ice  in  separate  slots.  The  top  of  the  jack  was 
60  cm  below  the  ice  surface  for  the  last  three  tests. 
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FIGURE  7.  AIR  JACK  INSERTION 


RESULTS 

Press ure -displacement  curves  for  the  Northwest  Research  Associates  tests  from  the  spring  of  1994 
are  shown  in  Figure  8.  Both  graphs  in  the  figure  show  the  same  four  tests,  with  different  pressure 
scales.  As  seen  in  the  figure,  the  strength  of  the  partially  consolidated  ice  in  the  ridge  was  highly 
variable  on  the  scale  of  the  air  jack.  Test  1  (in  the  top  layer  of  consolidated  ice)  saw  the  highest 
pressures.  It  was  stopped  because  of  pressure/volume  limitations  of  the  air  jack.  Given  the  precision 
of  the  calculated  displacement,  we  interpret  Test  1  as  having  negligible  displacement.  The  other  three 
tests  showed  varying  degrees  of  resistance  to  the  jack  and  were  stopped  because  of  displacement 
conditions.  All  three  show  numerous  displacement  events  at  erratic  time  intervals.  Test  3  had  a  large, 
abrupt  displacement  event  with  a  corresponding  rapid  pressure  drop  (due  to  the  expansion  of  the  air  in 
the  jack)  just  before  the  compressor  was  stopped. 

The  fluctuations  in  air  jack  pressure  (ram  force)  seen  here  are  similar  to  those  found  in  laboratory 
tests  (e.g.,  Urroz  and  Ettema,  1987)  despite  the  fact  that  the  ratio  of  the  ram  size  to  ice  block  size  is 
much  smaller  in  these  air  jack  tests.  These  fluctuations  result  from  a  combination  of  slip-stick  block 
morion,  breaking  of  cohesion  bonds,  and  the  ram  stiffness.  It  is,  therefore,  necessary  to  average  the 
results  of  numerous  tests  in  the  same  ridge  to  increase  the  precision  of  the  test  results. 
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FIGURE  8.  PRESSURE-DISPLACEMENT  CURVES  FOR  RIDGE  TESTS 


RIDGE  STRENGTH  HINDCAST  ANALYSIS 

Laboratory  test  data  from  Weiss  et  al.  (1981)  were  used  with  an  analysis  method  from  Coon  et  al. 
(1984)  to  hindcast  the  rigid-plastic  limit  load  for  Tests  2-4  in  the  partially  consolidated  portion  of  the 
keel.  As  shown  in  Table  2,  the  field  test  parameters  matched  quite  well  with  the  laboratory  test 
parameters.  The  ridge  ice  blocks  were  slightly  wanner  with  lower  salinity,  which  is  consistent  with 
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being  immersed  in  sea  water  for  two  weeks.  Weiss  reported  their  test  results  in  terms  of  the  Coulomb 
failure  parameters:  the  effective  angle  of  internal  friction  and  the  apparent  cohesion. 


TABLE  2.  ANALYSIS  PARAMETERS 


Parameter 

Laboratory  data  from 
Weiss  et  al.  (1981) 

SEV1I  ridge  on  April  14, 1995 

Ice  sheet  thickness  (cm) 

15 

14 

Void  ratio,  *n  (porosity) 

0.24  to  0.36 

0.35 

Ice  block  salinity  (%o ) 

10  to  11 

5 

Ice  block  temperature  (°C) 

-2.5  to  -5 

-1.4 

Displacement  rate  (m m/s) 

3 

0.4  | 

Load  area  (m2) 

0.55 

0.31  J 

Consolidation  state 

Unconsolidated 

Partially  consolidated  J 

Effective  angle  of  internal  friction,  <j) 

26° 

I 

Apparent  cohesion,  C  (kPa) 

2.3 

1 

Figure  9a  shows  the  plastic  failure  zone  assumed  for  this  analysis,  and  Figure  9b  shows  the  cor¬ 
responding  Coulomb  failure  surface  in  Mohr  space.  The  vertical  dimensions  in  Figure  9a  are  based  on 
the  average  values  from  the  six  drilled  holes  shown  in  Figure  3.  The  air  jack  or  ram  applies  a  uniform 
pressure,  Pa,  over  a  square  area  with  the  dimension,  a,  on  a  side  just  below  the  consolidated  layer. 
The  plastic  failure  zone  extends  down  and  away  from  the  ram. 


As  an  approximate  solution  to  the  problem  shown  in  Figure  9,  a  closely  related,  plastic  solution  for 
Coulomb  failure  caused  by  a  rectangular  indentor  will  be  used,  after  Coon  et  al.  (1984).  For  this 
problem,  the  loaded  area  extends  to  the  bottom  of  the  keel  as  shown  in  Figure  10,  in  which  the  ice 
sheet  is  inverted  to  expose  the  bottom  of  the  keel.  A  large  volume  of  ice  has  been  failed,  and  it  is 
assumed  that  the  load  is  on  a  half  space.  In  Figure  10,  there  are  three  failing  regions.  ABCDEF 
moves  in  the  direction  of  the  load  and  up.  The  regions  BHCF  and  ADGE  move  to  the  side  and  up.  If 
it  is  assumed  that  the  failure  load  of  Figure  10  is  the  same  as  that  of  Figure  9a,  then  the  indentor 
pressure,  P;.  is  smaller  than  the  equivalent  air  jack  pressure  because  of  the  load  area  difference 

Failure  load  =  P;ah  =  Paa2  (2) 

The  indentor  failure  pressure  is  given  in  Equation  17  of  Coon  et  al.  (1984)  as 

Pi  =CTc(l  +  C1(<{))|)  (3) 


where  ctc  is  the  unconfined  compressive  strength  of  the  Coulomb  material  and  C  [(<)>)  is  a  correction 
coefficient  derived  in  Coon  et  al.  (1984)  and  presented  graphically  in  Figure  1 1.  From  the  geometry  of 
the  triangle  in  Figure  9b  with  <(>  at  one  apex  and  2k  at  the  other,  the  unconfined  compressive  strength  is 
related  to  the  apparent  cohesion  and  effective  angle  of  internal  friction  given  in  Table  2: 

_  2C  COS  ()> 
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mm 


z 


a)  Plastic  failure  zone  in  the  keel 


x 


b)  Mohr  circles  for  stresses  showing  Coulomb  failure  surface 


FIGURE  9.  RIGID-PLASTIC  COULOMB  FAILURE  ANALYSIS 
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Indentation  F 


FIGURE  10.  ASSUMED  FAILURE  MECHANISM  FOR  MOHR-COULOMB  MATERIAL 
(ADAPTED  FROM  FIGURE  11  IN  COON,  EVANS,  AND  GIBSON,  1984) 


0  10  20  30  40  50  60 

Friction  Angle  of  Ice,  <)> 


FIGURE  11.  CORRECTION  COEFFICIENT  FROM  COON  ET  AL.  (1984) 
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FIGURE  11.  CORRECTION  COEFFICIENT  FROM  COON  ET  AL.  (1984) 


Combining  Equations  2, 3,  and  4  yields 

h|Ccos£( 

Table  3  includes  a  selection  of  solutions  to  Equation  5  using  the  dimensions  from  Figure  9a  (167 
cm  for  h,  and  56  cm  for  a).  Table  3  is  examined  row  by  row  in  the  following.  In  the  first  row,  the 
value  of  26°  for  <j)  from  Weiss  et  al.  (1981)  and  the  curve  in  Figure  11  are  used  to  find  a  value  of  1.8 
for  Cj(<j))  and  to  predict  a  value  of  140  kPa  for  the  air  jack  pressure.  The  average  jack  pressure 
achieved  in  the  three  tests  in  the  partially  consolidated  keel  was  90  kPa.  In  the  second  and  third  rows 
in  Table  3,  the  values  for  friction  angle  and  cohesion  were  adjusted  one  at  a  time  to  match  the  air  jack 
pressure  from  the  tests.  These  all  represent  reasonable  agreement  considering  the  scatter  in  the  three 
tests  and  the  approximations  made  in  the  prediction. 

In  addition,  Figure  12  compares  our  data  to  those  of  Weiss  et  al.  (1981)  for  tests  on  similar  size 
blocks.  Figure  12  is  adapted  from  Figure  3  in  Weiss  et  al.  (1981);  the  solid  circles  are  the  data  from 
their  slow  tests.  The  solid  line  marked  "slow"  is  their  best  fit  through  the  slow  data,  and  corresponds 
to  the  top  row  in  Table  3.  We  have  added  four  more  lines,  corresponding  to  the  bottom  four  rows  in 
Table  3.  The  last  two  rows  in  Table  3  are  a  lower  bound  and  an  upper  bound  for  the  Weiss  "slow" 
data.  The  predicted  air  jack  pressures  from  these  bounds  are  30  and  290  kPa;  these  two  pressures  are 
shown  as  horizontal  lines  in  Figure  8.  Tests  2  and  3  clearly  fall  within  this  range,  and  Test  4  may  have 
as  well  if  the  test  had  been  run  to  larger  displacements. 


TABLE  3.  ANALYSIS  RESULTS 


Identification 

Line 

number 

Effective  angle  of 
internal  friction,  d> 

Apparent  cohesion, 
C  (kPa) 

Correction 
coefficient,  C,(to) 

Air  jack  pressure  from 
Eqn.  5  (kPa) 

Weiss  et  al. 

1 

26 

2.3 

1.8 

140 

2 

26 

1.5 

1.8 

90 

3 

15 

2.3 

1.4 

90  | 

Lower  bound 

4 

10 

1 

1.1 

30  | 

Upper  bound 

5 

27 

4.5 

1.9 

290  1 

CONCLUSIONS  AND  DISCUSSION 

The  strength  of  the  ridge  keel  in  this  experiment  was  tested  at  only  one  time;  it  is  possible,  how¬ 
ever,  to  develop  a  clear  picture  of  the  keel  strength  from  its  formation  to  the  time  of  the  strength  tests, 
as  listed  in  Table  4.  On  March  30,  the  keel  was  formed  with  an  overall  depth  of  220  cm  comprised  erf 
a  14-cm  ice  sheet  and  a  206-cm  unconsolidated  keel.  On  April  14,  the  sheet  was  53  cm  thick  near  the 
test  slots,  leaving  167  cm  of  partially  consolidated  keel.  For  a  strength  comparison,  it  is  assumed  that 
the  unconfined  compressive  strength  of  the  ice  sheet  is  700  kPa,  which  is  the  horizontal  unconfined 
compressive  strength  of  unoriented  columnar  sea  ice  at  a  low  strain  rate  taken  from  American 
Petroleum  Institute  (API)  Bulletin  2N  (1988).  Although  we  stopped  our  test  in  this  ice  at  700  kPa,  the 
limit  load  pressure  of  an  air  jack  confined  in  this  ice  is  related  to  the  out-of-plane  failure  strength  of 
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CONFINING  PRESSURE  (  cr  ,KPa) 


FIGURE  12.  UNCONSOLIDATED  ICE  RUBBLE  SHEAR  STRENGTH  (ADAPTED  FROM 

WEISS  ET  AL,  1981) 


TABLE  4.  RIDGE  STRENGTH  EVOLUTION 


Parameter 

March  30 

April  14 

Ice  sheet  thickness  (m) 

0.14 

0.53 

Keel  thickness  (m) 

2.06 

1.67 

Ice  sheet  unconfined  compressive  strength 
(kPa) 

700 

700 

Keel  unconfined  compressive  strength  (kPa) 

7.4 

7.4 

Ice  sheet  strength  resultant  QcN/m) 

98  (87%) 

371  (97%) 

Keel  strength  resultant  (kN/m) 

15  (13%) 

12  (3%) 

383  (100%)  ! 

Total  strength  resultant  (kN/m) 

113  (100%) 

The  strength  resultant  is  the  product  of  the  strength  and  the  depth.  The  values  for  compressive 
strength  resultant  in  Table  4  are  calculated  using  the  unconfined  strengths  and  depths  in  Table  4.  It 
shows  that  the  total  unconfined  strength  of  a  unit  width  of  the  keel  has  increased  by  more  than  a  factor 
of  three  from  March  30  to  April  14  and  that  the  percentage  contribution  from  the  keel  has  decreased 
from  13%  to  only  3%.  Therefore,  it  appears  that  the  unconsolidated  part  of  the  keel  contributed  a 
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minor  portion  of  ridge  strength.  However,  it  is  important  to  remember  there  are  very  little  field  data 
on  which  to  base  this  conclusion. 

The  variation  of  the  stress-displacement  curves  for  the  partially  consolidated  keel  in  Tests  2-4  in 
Figure  8  imply  that  the  jack  size  of  0.31  m2  is  not  large  enough  to  measure  the  average  strength  in  a 
single  test.  This  is  not  surprising  given  that  the  blocks  had  a  thickness  of  0.14  m  with  an  average  area 
of  0.20  m2  with  a  standard  deviation  of  0.22  m2.  It  is  desirable  to  have  the  jack  size  be  much  larger 
than  the  block  size;  this,  however,  requires  very  large  equipment  to  make  the  slots  and  to  install  and 
inflate  the  large  jack.  It  is  reasonable  to  increase  the  size  of  a  rectangular  jack  to  about  1  m2  and  use  it 
with  a  system  that  can  be  transported  by  helicopter.  For  field  tests  such  as  these,  the  transportability  is 
very  desirable.  The  load  rate  can  be  increased  with  a  larger  compressor.  The  larger  jack  also  allows 
for  a  larger  displacement,  which  helps  determine  when  steady  state  motion  is  achieved  in  the  tests.  If 
many  tests  are  run,  then  the  average  should  be  same  as  the  result  from  a  larger  jack. 

Our  original  plan  was  to  test  a  ridge  in  the  fall,  winter,  and  spring  to  track  strength  with  consoli¬ 
dation.  That  was  not  possible  in  this  experiment  but  still  seems  to  be  a  worthy  objective. 
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ABSTRACT 

This  paper  describes  the  kinematics  and  ice-thickness  redistribution  of  pack  ice  while 
accounting  for  lead  and  ridge  orientation.  The  pack  ice  deformations  are  described  in  terms  of 
velocity  discontinuities  rather  than  invariants,  allowing  direct  interpretation  of  open-water 
production  and  ridging.  This  description  of  pack  ice  allows  an  oriented  ice-thickness  distribution 
to  account  for  the  growth  and  ridging  of  ice  in  refrozen  leads.  The  shear,  compressive,  and 
tensile  strengths  of  each  refrozen  lead  are  functions  of  the  lead  ice  thickness.  Results  of  simple 
model  calculations  for  ice-thickness  distribution  appear  reasonable.  Results  for  kinematics 
indicate  that  most  of  the  apparent  scatter  in  isotropic  calculations  of  open-water  production  can 
be  explained  by  accounting  for  lead  orientation.  The  theory  presented  in  here  is  a  small  part  of 
that  which  is  required  to  implement  a  large-scale,  anisotropic  model;  a  companion  paper  addresses 
the  architecture  of  such  a  model.  We  conclude  that,  to  have  high  resolution  with  cell  size  of  5  to 
100  km,  an  ice  dynamics  model  must  account  for  lead  orientation. 

INTRODUCTION 

In  polar  regions,  the  sea-ice  cover  acts  as  a  barrier  between  the  atmosphere  and  the  ocean, 
curtailing  the  transfer  of  thermal  and  mechanical  energy  as  well  as  curtailing  mass  transport.  The 
large-scale  motion  of  the  pack  ice  creates  a  small  amount  of  open  water,  in  the  form  of  leads.  In 
the  Arctic  winter  and  spring,  leads  are  linear  features  often  several  hundred  kilometers  long.  As 
these  leads  close  and  shear,  rafted  ice,  pressure  ridges,  shear  ridges,  and  rubble  are  created  along 
these  same  linear  features.  S AR  ice  motion  products  show  active  leads  as  linear  slip-lines  in  the 
pack  ice  (Stem  et  al.,  1995).  For  convenience,  we  will  refer  to  all  these  linear  features  as  leads. 

An  ice  dynamics  model  that  assumes  isotropy  in  the  ice  cover  can  adequately  characterize 
an  ice  cover  in  the  absence  of  leads,  e.g.,  summer/autumn  ice  with  multi-year  floes  surrounded  on 
all  sides  by  open  water  or  thin  ice.  However,  once  a  lead  forms,  the  ice  becomes  anisotropic  and 
is  best  described  by  a  model  that  accounts  for  this  property.  Careful  modeling  of  leads  in  sea  ice 
is  important  not  only  because  leads  provide  the  opportunity  for  the  heat  exchange  necessary  for 
ice  production  but  also  because  the  refrozen  lead  ice  controls  the  ice  dynamics  through  its 
strength  (Coon  et  al.,  1992).  The  dependence  of  lead  ice  strength  on  lead  ice  thickness  must  be 
established  from  experimental  data  supported  by  theories  such  as  Coon  et  al.  (1989). 

In  order  to  account  for  lead  and  ridge  orientation  for  pack  ice  on  the  scale  of  5  to  100  km, 
we  suggest  two  types  of  large-scale  pack  ice: 
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•  Isotropic  ice  is  consolidated  multi-year  and  thick  first-year  ice.  Much  of  it  is  the  jointed 
ice  described  above;  isotropic  ice  is  assumed  rigid  in  this  paper. 

•  Oriented  ice  includes  open  water,  new  ice,  nilas,  young  ice,  first-year  ice,  rafted  ice,  and 
ridged  ice  that  occur  in  long,  narrow  features.  Based  on  Advanced  Very  High  Resolution 
Radar  (AVHRR)  images,  these  features  are  assumed  linear  for  (grid,  element,  or  cell) 
scales  of  5  to  100  km.  All  pack  ice  deformations  are  assumed  to  occur  in  the  oriented  ice. 
In  addition,  the  linearity  assumption  implies  that  a  shearing  deformation  of  a  lead  does  not 
cause  rafting  or  ridging  of  the  oriented  ice  and  does  not  cause  any  failure  of  the  isotropic 
ice. 

In  the  next  sections,  we  describe  the  pack-ice  kinematics  accounting  for  lead  and  ridge 
orientation,  the  redistribution  of  oriented  ice-thickness  categories,  and  the  results  of  model 
calculations.  The  theory  presented  in  this  paper  is  only  a  small  portion  of  that  which  is  ultimately 
required  for  implementation  in  a  large-scale  ice  dynamics  model.  A  companion  paper  (Coon  et 
aL,  1997a)  includes  discussions  of  the  following: 

•  The  problems  to  be  addressed  before  an  oriented  ice  thickness  distribution  can  be 
implemented  in  an  ice  forecasting  or  climate  model. 

•  The  mechanical  response,  including  the  formation,  advection,  and  dissolution  of  the  leads. 

•  Methods  to  verify  model  results,  such  as  comparing  them  with  satellite-SAR  ice-motion 
products  for  lead  orientation,  ice  motion,  and  ice  thickness  distribution. 


KINEMATICS 


To  account  for  the  anisotropy  resulting  from  leads  and  ridges,  we  developed  a  new 
kinematic  description  of  the  pack  ice  motion  which  is  consistent  with  the  general  kinematic 
relations  for  pack  ice  presented  in  Pritchard  (1975).  The  orientation  of  the  velocity-discontinuity 
lines  where  leads  and  ridges  are  formed  is  referred  to  as  the  lead  direction.  Each  active  lead 
direction  contributes  to  the  plastic  deformation  of  the  pack  ice.  The  pack  ice  motion  is  therefore 
the  sum  of  deformations  along  individual  lead  directions  as  well  as  drift  and  rotation.  In  this 
paper,  we  assume  that  there  are  negligible  deformations  of  the  isotropic  ice  between  the  leads 
(e.g.,  elastic  deformations). 

As  illustrated  in  Figure  1,  we  align  a  coordinate  system  (i;,r|)  to  each  lead  direction  such 
that  the  £-axis  for  the  j-th  lead  direction  lies  at  the  angle  0j  to  the  x-axis  of  the  reference  (x,y) 
coordinate  system.  The  velocities  (u,v)  correspond  to  the  (x,y)  directions,  the  velocities  (U,V) 
correspond  to  the  (£,t|)  directions,  and  (Ax,Ay)  and  (A^,Aq)  define  the  size  of  the  cell  in  the 
reference  and  lead  coordinate  systems,  respectively.  In  the  lead  coordinate  system,  the  plastic 
stretching  parallel  to  the  lead  is  assumed  zero,  and  the  lead  width  is  assumed  locally  constant 
along  its  length.  We  can  express  velocity  derivatives  normal  to  the  lead  direction  (with  respect  to 
T|)  in  terms  of  smooth  gradients  or  in  terms  of  any  number  of  steps. 
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where  L  (=  A^)  is  the  average  length  of  many  parallel  leads  in  the  reference  element  or  the  length 
of  a  single  lead  in  the  reference  element  and  A  (=  AxAy  =  A^Ar))  is  the  area  of  the  reference 
element.  This  model  does  not  distinguish  whether  there  is  a  single  lead  crossing  the  element  or  a 
continuum  of  parallel  leads  that  cover  the  element. 

For  each  lead  direction,  we  form  the  plastic  velocity  gradient  L},  a  second-order  tensor 
with  four  components  in  the  coordinate  system  parallel  to  the  j-th  lead  direction: 
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The  components  of  the  corresponding  plastic  stretching  tensor  for  the  j-th  lead  direction 
are  as  follows: 


=i(Xj+x.T;= 


0 


V  ~  ' '  J\ 

f  T  ^  f^\l\ 


* 


au 


an 
av 

j 


(3) 


where  the  superscript  T  indicates  the  transpose. 

To  find  the  contribution  of  each  lead  motion  to  the  plastic  stretching  in  the  reference 
coordinate  system,  we  introduce  the  rotation  tensor  Qj  for  the  j-th  lead  direction.  The  rotation 
tensor  is  written  in  terms  of  0j,  the  angle  from  the  reference  x-axis  to  the  %-axis  parallel  to  the  j-th 
lead  direction: 
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To  form  D,  the  plastic  stretching  in  the  reference  coordinate  system,  we  rotate  each  of  the 
N  lead  tensors  to  the  reference  coordinate  system  to  form  Dj  and  add  the  contributions  of 
individual  lead  directions,  so 


D  =  ID,  =£0^01 

j=i  j=i 


(5) 


The  sum  of  the  diagonal  terms  is  Di,  the  plastic  divergence,  which  depends  only  on  the 
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In  words,  the  total  divergence  is  equal  to  the  sum  of  the  normal  velocity  derivatives  associated 
with  each  lead  direction.  For  AxAy  equal  to  A^Ar|,  the  rate  of  area  change  in  the  cell  is  therefore 
equal  to  AxAyDi. 

In  this  anisotropic  model,  the  motions  of  the  leads  are  not  uniquely  determined  in  terms  of 
the  plastic  stretching  tensor;  additional  information  is  needed.  As  is  discussed  in  a  companion 
paper  (Coon  et  al.,  1997a),  the  solution  to  an  anisotropic  plasticity  ice  dynamics  problem  will 
accomplish  the  following: 

1.  Place  the  stress  state  on  the  anisotropic  failure  surface  when  plastic  failure  occurs. 

2.  Identify  the  specific  failure  mechanisms  involved,  e.g.,  closing  lead  1  while  shearing  lead  2. 

3.  Allocate  the  stretching  D  to  the  contribution  of  the  j-th  lead  direction  to  the  plastic 
stretching  Dj. 

The  inverse  rotation  of  Equation  (5)  gives  the  plastic  stretching  components  in  lead 
coordinates: 

D]  =Q[DjQj  (7) 


Using  Equation  3,  the  components  of  the  matrix  yield  the  velocity  derivatives  normal  to 
the  lead  direction  in  terms  of  the  components  of  the  tensor  Dj: 


=  2  Dj  xycos  20J  -(D^-Dj>yy)sin  20j 


(8) 


=  Dj,xx+Dj,yy 


In  this  model,  the  velocity  derivative  normal  to  the  lead  directions  Vj  is  the  only  part  of  the 
kinematics  that  affects  ice-thickness  redistribution  in  the  leads.  This  result  depends  on  the 
assumption  that  the  leads  are  locally  straight  (Le.,  within  a  cell)  so  that  shear  motions  across  the 
lead  do  not  redistribute  the  ice  thickness,  as  opposed  to  the  usual  interpretation  expressed  in 
Thorndike  (1975). 


ICE-THICKNESS  REDISTRIBUTION 

The  failure  of  isotropic  ice  results  in  open  water  in  leads  or  at  least  slip  lines  with  little 
open  water.  In  our  experience,  rafting  2-m  ice  sheets  is  rare  and  ridges  built  from  a  2-m  ice  sheet 
are  unheard  of.  Although  an  occasional  2-m  ice  block  is  kicked  up,  we  have  never  seen  ridges 
with  10-m  sails  and  40-m  keels  made  of  2-m  blocks.  Hence,  we  do  anticipate  the  need  for  ridging 
isotropic  ice  in  this  model. 

Subfreezing  temperatures  will  cause  oriented  ice  to  form  in  leads,  and  ice  deformations 
will  cause  a  distribution  of  ice  thickness  within  these  leads.  The  distribution  of  ice  thickness  at 
each  lead  orientation  can  be  described  much  as  others  have  described  the  distribution  of  ice 
thickness  (e.g.,  Thorndike  et  al.,  1975).  Figure  2  illustrates  how  we  might  characterize  these  ice- 
thickness  distributions  assuming  isotropic  ice  has  one  ice-thickness  category  and  oriented  ice  has 
four  ice-thickness  categories  in  two  lead  systems.  The  thickest  oriented  ice  category  may  be 
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thicker  than  the  isotropic  ice,  since  the  majority  of  the  isotropic  ice  is  multi-year  ice.  Certain 
properties  define  the  i-th  ice-thickness  category  in  the  j-th  lead: 

•  Lead  or  lead  orientation  0j 

•  Ice  thickness  separating  adjacent  ice-thickness  categories  hy 

•  Area  fractions  of  each  ice-thickness  category  AGy 

•  Mean  ice  thickness  in  each  ice-thickness  category  hy 

We  assume  the  refrozen  lead  ice  to  be  homogeneous  (within  each  ice-thickness  category). 
In  addition,  we  assume  the  leads  are  straight  so  that  a  shearing  deformation  of  the  lead  deforms 
neither  isotropic  ice  nor  oriented  ice.  We  associate  four  ice-thickness  categories  with  each  lead 
direction,  applying  the  term  ’oriented  thickness  distribution’  to  distinguish  it  from  the  usual 
(isotropic)  ice-thickness  distributions  (e.g.,  Thorndike  et  al,  1975). 

When  the  air  temperature  is  below  freezing,  new  ice  grows  in  the  leads;  its  growth  rate  is 
often  modeled  simply  as  a  function  of  air  temperature  (Anderson,  1961;  Zubov,  1944;  Lebedev, 
1940),  e.g.,  from  Lebedev, 

h  =  0.0 133(  FDD) 0-58  (9) 


where  h  is  the  ice  thickness  in  meters  and  FDD  is  the  number  of  freezing-degree  (C)  days  since 
the  open  water  was  formed.  Alternatively,  Thorndike  et  aL  (1975)  give  the  ice  growth  rate  as  a 
function  of  time  of  year  for  the  central  Arctic,  and  Maykut  and  Untersteiner  (1971)  compute  it 
from  the  surface  energy  balance  and  ice  thermal  conductivity. 

Based  on  Pritchard  and  Coon  (1981),  we  have  developed  the  mathematical  equations  and 
calculation  procedures  for  the  ice-thickness  distribution  with  four  ice-thickness  categories  in  each 
lead.  The  pertinent  parameters  for  each  ice-thickness  category  are  in  Table  1.  We  define  ho  as 
the  (fixed)  ice  thickness  separating  the  open-water  and  thin-ice  categories,  hjj  as  separating  thin 
ice  and  medium  ice  in  the  j-th  lead,  and  h2j  as  separating  medium  ice  and  thick  ice  in  the  j-th  lead. 
We  let  fo  be  the  growth  rate  of  ice  of  thickness  ho  calculated  from  the  time-derivative  of  Equation 
(9).  Except  for  the  open-water  category,  the  definition  of  the  ice-thickness  category  itself 
accounts  for  thermal  growth;  Le.,  the  category  limits  hqj  and  h2j  increase  at  the  growth  rate  of  ice 
of  the  same  thickness.  By  assuming  a  linear  ice-thickness  distribution  within  the  open-water  and 
thin-ice  categories,  Pritchard  and  Coon  derived  the  rate  of  area  change  due  to  thermal  growth,  d> 
0j.  For  the  j-th  lead. 


'foAG0j 

h0 

f0AGy 


if  f0  >0  (freezing  conditions 
if  f0  <0  (melting  conditions 


(10) 


where  AGoj  is  the  area  fraction  of  open  water  in  the  j-th  lead  and  AGy  is  the  area  fraction  of  thin 
ice  in  the  j-th  lead. 
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When  the  j-th  lead  is  opening,  i.e.,  when  the  lead’s  velocity  gradient 


is  positive. 


the  following  equations  describe  the  total  time  rate  of  change  of  area  fraction  in  the  four  ice 
categories: 
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where  AG2j  is  the  area  fraction  of  medium  ice  in  the  j-th  lead  and  AG3j  is  the  area  fraction  of  thick 
ice  in  the  j-th  lead.  Equation  (11)  provides  that  the  source  for  open  water  is  opening  the  lead  and 
that  area  is  transferred  between  the  first  two  categories  by  freezing  or  melting.  These  Lagrangian 
time  rate  of  change  equations  can  be  easily  converted  to  the  usual  Eulerian  form  by  using  the 
Eulerian  derivative  (see,  e.g.,  Liepmann  and  Roshko,  1957): 


Pf  _  df 
Dt  dt 


+  v-Vf 


The  ice  category  area  fractions  and  the  components  of  the  lead  direction  (cos0,  sin9)  are  advected 
along  with  the  ice. 

fdv} 

When  the  j-th  lead  is  closing,  i.e.,  when  — —  is  negative,  the  following  equations 

Wi 

describe  the  rate  of  change  of  area  fraction  in  the  four  ice  categories: 

DAG,,  h„  ^ 

— T — =  TT  r - f^uolAGoj]  -<E>oj 

dt  ^Jj^u-ho  j 

3=  u0(AG0j) 

dag,  rav^  h,-h0  ^ 

dt  ^  +?A-h_U^AG^  t-uo(AGo,)l 

)  2'-  “  '  .  (13) 

nAr  [l-u0(AGoj)] 

DAG2j  _[  3V  h2j“hij 

^Tl  I.  h,,  /  \  r  m  r  /  yt 


+  d>oj 


V  h  +r-^-u0(AG2j)  &-u0(AG,)]  [l-u0(AG0j)] 

^  n3j  n2j 

&-u0(ag„)]  [i-u„(ag0j)] 

dt  (^drij.l  h3j  -h2j 


7 


where  h-bar  indicates  the  average  ice  thickness  within  the  ice-thickness  category.  Equation  (13) 
provides  the  following: 

1.  Freezing  and  melting  transfer  the  area  between  the  first  two  categories. 

2.  Rafting  transfers  the  area  from  the  open-water  category  to  the  thin  ice  category. 

3.  Rafting  transfers  the  area  from  the  thin-ice  category  to  the  medium-ice  category  when  the 
open  water  is  gone. 

4.  Ridging  transfers  the  area  from  the  medium-ice  category  to  the  thick-ice  category  when 
the  open  water  and  thin  ice  are  gone. 

For  convenience,  we  have  introduced  the  unit  step  function  ua(x)  such  that 


ua(x) 


fO  when  x  <  a 
|l  when  x  >  a 


(14) 


The  average  ice  thickness  in  each  category  increases  at  the  growth  rate  of  ice  of  the  same 
thickness,  except  when  opening  occurs.  In  that  instance,  the  average  thickness  of  the  ice  in  the 
open- water  category  is  reduced  by  the  contribution  of  true  open  water  to  the  average. 


RESULTS  OF  MODEL  CALCULATIONS 


The  behavior  of  these  ice-thickness  redistribution  equations  is  illustrated  by  a  simple 
example  for  one  lead;  the  results  are  shown  in  Figure  3.  In  this  example,  a  cell  of  pack  ice  initially 
has  only  ice  in  the  thickest  category  (white).  At  an  air  temperature  of  -40°  C,  the  lead  area 
fraction  increases  steadily  at  the  rate  of  10'5  s'1  for  two  days,  then  closes  steadily  at  the  same  rate 
for  two  days.  Initially,  the  open-water  area  fraction  (black)  increases  rapidly  at  the  expense  of  the 
thickest  category.  Gradually,  freezing  moves  ice  from  the  open-water  category  to  the  thin-ice 
category  (dark  gray).  During  closing,  the  open-water  ice-thickness  category  is  rafted  into  the 
thin-ice  category.  When  the  open-water  category  is  fully  consumed,  the  thin-ice  category  is 
ridged  into  the  medium-ice  category  (light  gray).  When  the  thin-ice  category  is  consumed,  the 
medium-ice  category  is  ridged  into  the  thick-ice  category  (white).  If  the  cell  size  is  fixed 
(Eulerian),  then  there  is  mass  flux  into  and  out  of  the  cell. 

We  also  compared  open-water  production  using  the  anisotropic  model  to  open-water 
production  using  the  usual  isotropic  models.  Isotropic  models  approximate  the  open-water 
production  as  a  single  function  of  the  ratio  of  the  stretching  invariants  (Thorndike  et  aL,  1975; 
Stem  et  aL,  1995).  For  example,  Thorndike  et  aL  approximated  open-water  production 
normalized  by  the  stretching  magnitude  as  a  function  of  the  angle  from  pure  divergence  to  the 
invariant  stretching  vector: 

°E=£  =  ao=l(l  +  cos0)  (15) 


where 


Di=Dx* 

0  =  tan1 


+  D 


D, 

D, 


Dn=V(D»-DJ+4D«> 

|d|=Vd?+d; 


(16) 


8 


Figure  6  in  Stem  et  aL  shows  considerable  scatter  in  the  opening  found  from  342  strips  of 
SAR  data  for  which  the  time  interval  was  three  days  when  compared  with  three  isotropic 
functions,  including  Equation  (15).  They  explained  the  scatter  in  the  data  as  “measurement  error” 
because  restricting  the  data  to  the  more  active  scenes,  defined  as  those  having  ID!  >  0.07,  reduced 
the  scatter.  We  suggest  that  accounting  for  lead  orientation  may  explain  much  of  the  apparent 
scatter. 

Based  on  the  actions  of  one  or  two  lead  systems,  we  generated  families  of  curves  for 
opening  as  a  function  of  0  with  a  fixed  angle  between  two  lead  systems  of  34°  (angle  from 
Cunningham  et  aL,  1994).  These  results  are  applicable  over  all  scales  for  which  the  leads  are 
linear,  e.g.,  from  5  to  100  km;  these  results  are  therefore  comparable  to  a  SAR  pixel  or  an  entire 
SAR  strip.  From  Equations  (3),  (4),  (5),  the  total  stretching  is  the  sum  of  N  lead  stretching 
contributions  rotated  to  the  reference  coordinate  system: 
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Applying  Equation  (17)  to  the  case  of  two  leads  at  an  angle  to  each  other  can  be  used  to 
demonstrate  that  any  stretching  can  be  accomplished  from  opening,  closing,  and  shearing  the  two 
leads.  This  result  implies  that  a  third  lead  cannot  be  formed  until  one  of  the  first  two  leads  is 
totally  closed. 


The  open-water  production  is  calculated  directly  from  the  positive  values  of 


Lead  1  opening  with  Lead  2  closing: 


Both  leads  opening: 
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Figures  4,  5,  and  6  show  that  open-water  production  depends  strongly  on  the  presence  of 
lead  systems.  Figure  4  shows  curves  on  which  the  ratio  of  lead-2  opening  velocity  to  lead-1 
opening  velocity  is  varied  over  a  wide  range.  The  curve  corresponding  to  the  opening  function  in 
Equation  (15)  is  shown  for  reference  as  Thorndike.  The  case  in  which  only  lead  1  is  active  and 
opening  is  shown.  Only  a  limited  range  of  values  for  0  are  possible  because  the  stretching  parallel 
to  the  lead  is  zero,  forcing  Dn/Di  to  be  greater  or  equal  to  one  and  therefore  45°  <  ©  <  135°.  The 
case  in  which  lead  1  is  opening  and  shearing  to  the  right  along  with  lead-2  closing  with  no  shear 
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passes  close  to  Equation  (14).  The  case  in  which  lead  1  is  opening  with  no  shear  and  lead  2  is 
closing  with  no  shear  produces  small  values  of  opening.  The  case  in  which  lead  1  is  opening  and 
shearing  to  the  left  along  with  lead-2  closing  with  no  shear  produces  large  values  of  opening. 
Such  large  values  of  opening  can  occur  only  until  lead  2  is  closed,  which,  depending  upon  the 
rates  involved,  may  take  hours  or  days.  This  limitation  suggests  that  the  reduced  scatter  seen  by 
Stem  et  aL  when  they  considered  only  the  more  active  scenes,  results  from  closing  the  second 
lead  during  the  3-day  period. 

Figure  5  shows  curves  on  which  the  ratio  of  lead-1  shearing  to  opening  is  varied  over  a 
wide  range.  Lead-2  shearing  is  zero.  The  opening  function  from  Thorndike  et  al.  is  again  shown 
for  reference.  The  cases  have  the  ratio  of  lead-2  opening  to  lead-1  opening  fixed  at  the  indicated 
value.  This  family  of  curves  shows  opening  even  at  very  high  values  of  0.  Figure  6  shows  curves 
on  which  the  ratio  of  lead-2  shearing  to  opening  is  varied  over  a  wide  range.  Lead- 1  shearing  is 
zero.  Equation  (14)  from  Thorndike  et  aL  is  again  shown  for  reference.  The  cases  have  the  ratio 
of  lead-2  opening/closing  to  lead-1  opening  fixed  at  the  indicated  value;  lead  2  is  opening  when 
the  ratio  is  positive  and  closing  when  the  ratio  is  negative.  This  family  of  curves  shows  extremely 
large  opening  at  low  values  of  0.  The  three  figures  together  show  that  reasonable  values  for  the 
various  velocity  ratios  (less  than  ±10:1)  produce  a  broad  range  of  opening  normalized  by  the 
stretching  magnitude  for  most  values  of  0. 

The  results  of  these  calculations  indicate  that  the  apparent  scatter  in  open-water 
production  can  be  explained  by  accounting  for  lead  orientation.  There  is  no  longer  the  one-to-one 
mapping  of  stretching  to  open-water  production  as  there  is  with  an  isotropic  model  Additional 
constraints  imposed  by  the  relations  between  stress,  stretching,  and  the  corresponding  plastic 
failure  surface,  however,  may  provide  bounds  on  the  opening  not  seen  in  these  figures. 

CONCLUSIONS 

We  have  developed  a  formulation  for  pack  ice  kinematics  and  an  oriented  ice-thickness 
distribution  that  accounts  for  the  orientation  of  existing  leads  and  ridges.  It  is  concluded  that: 

1.  The  plastic  deformations  of  pack  ice  can  be  expressed  as  the  sum  of  contributions 
from  the  active  leads  and  ridges. 

2.  Two  active  lead/ridge  systems  can  produce  all  possible  deformations,  so  a  third 
cannot  form  until  one  of  the  first  two  is  closed. 

3.  The  opening/closing  of  individual  lead/ridge  systems  can  be  determined  from 
solving  the  plastic  equations  for  pack  ice  but  cannot  be  determined  solely  from  the 
stretching  tensor. 

4.  The  opening/closing  of  individual  lead/ridge  systems  can  be  used  to  directly 
redistribute  the  ice  thickness. 

This  oriented  thickness  distribution  will  be  appropriate  for  use  in  ice  dynamics  models,  in 
climate  dynamics,  in  pollution-transport  studies,  and  in  ice-motion  and  noise  forecasting.  In 
addition,  we  have  concluded  that  it  is  inappropriate  to  scale  down  an  isotropic  large-scale  ice 
dynamics  model  to  a  cell  size  of  a  few  kilometers  without  separately  modeling  the  isotropic  ice 
and  the  oriented  ice. 
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NOTATION 

D  plastic  stretching  tensor  (2nd-order)  (1/s) 

Dj  contribution  to  plastic  stretching  tensor  due  to  j-th  lead  direction  (1/s) 

Dj  components  of  Dj  in  j-th  lead  coordinates  (2nd-order  matrix)  (1/s) 

Di,  Dn  first  invariant  (divergence)  and  second  invariant  (shear)  of  plastic  stretching  (1/s) 

Dxx,  Dxy,  Dyy  components  of  plastic  stretching  tensor  (1/s) 
f0  growth  rate  of  ice  of  thickness  h0  (m/s) 
h  ice  thickness  (m) 

ho  ice  thickness  separating  open  water  from  thin-ice  category  (fixed)  (m) 

hy  greatest  ice  thickness  in  i-th  category  in  j-th  lead  system  (m) 

hjj  average  ice  thickness  in  i-th  category  of  j-th  lead  system  (m) 

Li  contribution  of  velocity  gradient  tensor  due  to  j-th  lead  (1/s) 

N  number  of  lead  directions 

Qj  orthonormal  rotation  tensor  (2nd-order)  for  j-th  lead  direction 
u,  v  velocity  components  in  reference  coordinate  system  (m/s) 

U,  V  velocity  components  in  lead  coordinate  system  (m/s) 

ua  unit  step  function 

x,  y  reference  coordinate  system  (m) 

a0  opening  coefficient  from  Thorndike  (1975) 

AGy  area  fraction  of  i-th  ice-thickness  category  in  the  j-th  lead  system  (m) 

AU,  AV  velocity  jumps  across  Ar|  in  lead  coordinate  system  (m/s) 

Ax,  Ay  cell  dimensions  in  reference  coordinate  system  (m) 

At,,  At|  cell  dimensions  in  lead  coordinate  system  (m) 

0  in  stretching  invariant  space,  the  angle  from  pure  divergence  to  the  invariant  stretching 
vector,  equal  to  the  inverse  tangent  of  the  shear  divided  by  the  divergence  (radians) 

9j  angle  from  x-axis  to  j-th  lead  direction  (radians) 

£,T|  lead  coordinate  system  (m) 

d>0j  rate  of  area  change  due  to  thermal  growth  (1/s) 
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LIST  OF  FIGURES 

Figure  1.  Coordinate  system  and  velocity  components  of  a  rotated  cell  aligned  with  the  lead. 

Figure  2.  Ice-thickness  distribution  includes  the  oriented  ice  in  each  lead  direction  separately 
along  with  the  isotropic  ice. 

Figure  3.  Example  of  redistribution  of  ice  thickness  expressed  as  area  fractions  for  four  ice- 
thickness  categories  resulting  for  opening  and  closing  a  lead. 

Figure  4.  Accounting  for  lead  orientation  results  in  families  of  curves  for  opening  as  a  function  of 
©.  Failure  of  the  isotropic  ice  falls  on  the  curve  for  one  active  lead.  The  three  higher 
loops  result  from  varying  the  opening  of  the  second  lead  for  three  values  of  shear  to 
opening  for  the  first  lead. 
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Figure  5.  Family  of  curves  for  opening  formed  by  holding  shear  in  the  second  lead  at  zero  and 
selecting  specific  values  of  the  ratio  of  opening  in  lead  2  to  lead  1. 

Figure  6.  Family  of  curves  for  opening  formed  by  holding  shear  in  the  first  lead  at  zero  and 
selecting  specific  values  of  the  ratio  of  opening  in  lead  2  to  lead  1. 


Table  1.  Ice-Thickness  Categories 


j  Ice-thickness 
Categories 

Ice  Types 

Area  Fraction 
in  j-th  Lead 

Average 

Thickness 

Open  water 

Open  water,  new 
ice,  and  nilas 

AGoj 

0  <  h  <  h0 

h0 

Thin  ice 

Young  ice  and 
rafted  ice 

AGij 

h0  <  h  <  hXj 

Klj 

Medium  ice 

mmm 

AG2j 

hjj  <  h  <  h2j 

h2j 

Thick  ice 

Ridged  ice  and 
multi-year  ice 

AG3j 

h2j<h 

*3j 
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tem  and  velocity  components  of  a  rotated  cell  aligned  with  the  slip  line. 
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Figure  2.  Ice  thickness  distribution  includes  the  oriented  ice  in  each  slip-line  directions  separately 
along  with  the  isotropic  ice. 
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Figure  3.  Example  of  redistribution  of  ice  thickness  expressed  as  area  fractions  for  four  ice-thickness  categories  resulting  for  opening  and 
closing  a  lead. 
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Figure  4.  Accounting  for  lead  orientation  results  in  families  of  curves  for  opening  as  a  function  of  0.  Failure  of  isotropic  ice  falls  on 
the  curve  for  one  active  lead.  The  three  higher  loops  result  from  varying  the  closing  of  the  second  lead  for  three  values  of  the  shear  to 
opening  ratio  for  the  first  lead. 
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Figure  5.  Family  of  curves  for  opening  formed  by  holding  shear  in  the  second  lead  at  zero  and  selecting  specific  values  of  the  ratio 
of  opening  in  lead  2  to  lead  1. 


o 

oo 


IQI  /  SuiuadQ 


Figure  6.  Family  of  curves  for  opening  formed  by  holding  shear  in  (lie  first  lead  at  zero  and  selecting  specific  values  of  the  ratio 
of  opening  in  lead  2  to  lead  1. 


